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SUMMARY 

The  Method  of  Moments  (MoM)  is  a  powerful  computational  electromagnetics  method  for 
analyzing  antennas  in  situ.  In  a  previous  report,  we  presented  a  new  method  for  evaluating  the 
inner  integral  of  the  impedance  matrix  element  in  the  traditional  Rao-Wilton-Glisson  formulation 
of  MoM  for  the  electric  field  integral  equation.  In  this  report,  we  extend  our  method  to  the 
magnetic  field  integral  equation.  We  show  how,  in  this  case,  the  inner  integral  can  be  written  in 
tenns  of  two  scalar  integrals,  a  surface  integral  over  the  integration  triangle  (IT)  and  a  line 
integral  over  the  boundary  of  the  IT.  Just  as  with  the  original  integral,  neither  of  these  integrals 
can  be  evaluated  analytically.  In  our  method,  we  bypass  this  obstacle  by  replacing  the  original 
integrand  (modified  by  a  constant  phase  factor)  by  its  Taylor  series  and  by  keeping  enough  tenns 
to  guarantee  a  number  of  significant  digits  in  the  integration  outcome.  We  have  accomplished 
this  for  both  the  surface  integral  and  the  line  integral.  We  present  a  systematic  derivation  of  the 
formulas  for  the  two  integrals  and  we  conduct  an  extensive  testing  of  our  results  against 
quadratures  and  cubatures. 
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INTRODUCTION 

The  Method  of  Moments  (MoM)  is  one  of  the  principal  methods  for  electromagnetic  simulations 
in  the  frequency  domain.  In  this  method,  we  express  an  unknown  current  density  as  a  linear 
combination  of  known  functions  and  the  objective  of  the  MoM  is  to  detennine  the  coefficients  in 
this  expansion  by  minimizing  the  square  of  the  modulus  of  the  residual  error.  It  ultimately  leads 
to  a  system  of  linear  algebraic  equations,  the  solution  of  which  yields  the  values  of  the 
coefficients  of  the  current  density  expansion.  The  elements  of  the  matrix  (known  as  the 
impedance  matrix  (IM))  are  expressed  as  double  surface  integrals  over  flat,  triangular  regions  of 
space.  In  a  previous  report  (reference  1),  we  computed  the  inner  of  these  integrals  to  a  prescribed 
precision  for  the  Electric  Field  Integral  Equation  (EFIE).  This  work  appeared  also  as  a  paper  but 
in  a  more  condensed  form  (reference  2).  In  the  present  study,  we  compute  the  inner  integral  for 
the  Magnetic  Field  Integral  Equation  (MFIE).  On  both  occasions,  the  object  under  consideration 
is  a  perfect  conductor. 

Our  motivation  for  including  the  MFIE  in  this  prescribed-precision  analysis  is  that,  usually,  the 
MFIE  is  combined  linearly  with  the  EFIE  into  one  equation  so  as  to  eliminate  the  effect  of  the 
internal  resonances  of  the  target  structure  (reference  3,  Theorem  4.47).  Moreover,  the  MFIE 
operator  appears  also  in  the  pair  of  integral  equations  that  describe  the  behavior  of  penetrable 
materials,  both  lossless  and  lossy.  For  such  materials,  the  wave  number  is  generally  complex. 

The  motivation  for  computing  the  matrix  elements  to  prescribed  precision  is  presented  in 
reference  1  and  repeated  here.  Until  recently,  the  quantities  of  interest  in  electromagnetic 
simulations  (for  both  scattering  and  radiation  problems)  were  the  far  electric  and  magnetic  fields. 
From  these,  we  could  compute  the  radar  cross  section  of  a  target  or  the  radiation  patterns  of  an 
antenna.  This  allowed  for  a  considerable  degree  of  error  in  computing  the  elements  of  the  IM  in 
the  MoM  because  of  the  error  smoothing  effect  of  the  near-  to  far- II eld  transformation 
(integration).  Moreover,  due  to  computer  hardware  limitations,  the  size  of  the  IM  was  small 
enough  so  that  the  round-off  error  did  not  have  a  severe  effect  on  the  accuracy  of  the  solution. 

In  recent  years,  the  MoM  is  being  applied  to  problems  where  near-field  information  is  needed 
(e.g.,  antennas  and  antenna  arrays).  This  necessitates  a  more  accurate  computation  of  the 
elements  of  the  IM  than  when  only  the  far  fields  are  of  interest.  Additionally,  advances  in 
computer  hardware  allow  us  to  address  problems  that  result  in  an  IM  system  with  millions  of 
unknowns.  Thus,  the  effect  of  the  round-off  error  becomes  more  pronounced.  Both  of  these 
reasons  lead  us  to  the  conclusion  that  the  more  accurate  the  representation  of  the  IM  is,  the  better 
the  quality  of  the  solution.  A  third  observation  we  made  recently  is  the  following.  In  computing 
the  interior  fields  due  to  a  slit  on  a  perfectly  conducting  enclosure,  we  noticed  that  the  result  was 
corrupted  by  numerical  noise.  To  verily  this,  we  considered  a  hollow  cube  with  perfectly 
conducting  walls  under  the  exterior  influence  of  a  plane,  electromagnetic  wave.  We  computed 
the  surface  current  density  using  MoM,  and  used  the  result  to  compute  the  fields  inside  the 
cavity.  The  fields  inside  the  cavity  should  ideally  be  equal  to  zero  since  the  walls  of  the  cavity 
are  perfectly  conducting.  We  found  that  the  fields  inside  reached  as  high  as  a  tenth  of  the 
amplitude  of  the  incident  field  when  operating  the  MoM  in  its  default  mode.  We  had  to  increase 
the  number  of  unknowns  seventy-fold  before  we  could  get  the  interior  field’s  amplitude  to  be  a 
hundredth  of  that  of  the  exterior  field  (24  dB  down,  to  be  precise).  The  increase  in  number  of 
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unknowns  was  due  to  the  increase  of  the  grid  density.  This  results  in  a  more  accurate  calculation 
of  the  impedance  matrix  elements  since  the  region  of  integration  becomes  smaller  for  each 
element  and,  hence,  the  variation  of  the  integrand  over  it  is  smaller.  The  alternative  to  an 
extremely  fine  grid  and  the  attendant  expense  in  time  and  hardware  is  to  be  able  to  use  a  standard 
grid  (the  side  of  a  triangle  less  or  equal  to  a  tenth  of  a  wavelength)  and  use  the  present  method  to 
calculate  the  inner  integral  to  a  specified  number  of  significant  digits  (SD).  This  is  the  central 
theme  of  the  present  study:  the  computation  of  the  inner  surface  integral  of  the  IM  element  to 
prescribed  precision  for  the  MFIE. 

In  Section  1,  we  apply  to  the  MFIE  the  MoM  in  the  manner  of  Rao-Wilton-Glisson  (RWG) 
(reference  4).  We  then  proceed  to  show  that  the  inner  surface  integral  of  the  IM  can  be  written  in 
tenns  of  a  surface  integral  over  the  integration  triangle  (IT)  and  a  line  integral  around  the  IT’s 
boundary.  Both  of  these  integrals  are  scalar.  In  the  line  integral,  we  have  as  integrand  the  firee- 
space  Green's  function  while,  in  the  surface  integral,  the  normal  derivative  of  the  same  function. 
The  approach  we  employ  to  calculate  these  integrals  to  a  prescribed  precision  is  analogous  to 
that  in  reference  1:  we  approximate  each  of  the  two  integrands  using  Taylor’s  Theorem  with  a 
Remainder  (reference  5).  This  theorem  allows  us  to  approximate  the  integrand  (that  we  do  not 
know  how  to  integrate  analytically)  by  a  polynomial  that  we  can  integrate  analytically.  The 
degree  of  approximation,  i.e.,  the  number  of  SD  to  which  the  approximation  is  accurate  is  known 
and,  hence,  so  is  that  of  the  integration.  The  result  of  this  evaluation  is  given  in  terms  of 
elementary  transcendental  functions  for  which  there  exist  robust  computational  algorithms. 

In  Section  2,  we  deal  with  the  surface  integral  and  show  that  it  can  be  computed  using  a  three- 
term  iteration  that  has  as  its  starting  point  the  solid-angle  integral.  It  also  involves  a  known  term 
that  is  related  to  the  integral  of  the  Green's  function  over  the  IT  and  which  we  have  computed  in 
reference  1  by  converting  it  to  a  line  integral.  We  summarize  the  results  from  reference  1  in 
Section  5.  In  Section  3,  we  turn  to  the  evaluation  of  the  line  integral  and  analyze  it  in  quite  the 
same  way  we  analyzed  the  integral  of  the  Green’s  function  over  the  IT  in  reference  1 .  In  Section 
4,  we  discuss  two  formulas  for  computing  the  solid  angle  integral.  In  Section  6,  we  validate  our 
formulas  using  the  Gauss-Kronrod  quadrature  (GKQ)  (reference  6)  as  it  appears  in  Matlab 
(reference  7)  for  the  line  integral  and  the  cubatures  from  reference  1  for  the  surface  integral. 

In  conclusion,  we  have  extended  our  work  on  calculating  the  inner  integral  of  the  IM  element  to 
prescribed  precision  to  the  MFIE.  We  mention  that  the  algorithms  we  developed  are  not  for  all 
observation  points  (OP)  but,  rather,  for  those  close  to  the  IT.  We  hope  that,  in  a  subsequent 
report,  we  will  carry  the  same  analysis  as  in  Part  II  of  reference  1  to  determine  the  spherical 
region  about  the  IT’s  centroid  where  it  is  necessary  and/or  advantageous  to  use  the  present 
method,  and,  outside  of  which  it  is  preferable  to  use  the  GKQ  for  the  line  integral  and  cubatures 
for  the  surface  integral. 

We  know  of  no  other  method  that  calculates  the  inner  integral  of  the  IM  element  for  MFIE  to  a 
prescribed  precision. 
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SECTION  1:  THE  MAGNETIC  FIELD  INTEGRAL  EQUATION 

The  MFIE  (also  known  as  Maue’s  integral  equation)  for  a  closed  and  bounded  perfectly 
conducting  surface  S  is  (reference  8,  pp.  354-355) 


—  J'(r')  =  J"'f(r') -n'  x  [  J'(r)  x  Vg(r,  r')dS  ,  r  'eS  (1.1) 

2  s 

where  n  is  the  exterior  unit  normal  on  S  and 

-z£|r-r'|  _ 

=  - 7  ,  k  =  (o^£()Ju0  (1.2) 

4^-r-r 

with  r  the  integration  point  and  r'  the  observation  point 

r  =  xx  +  yy  +  zz ,  r'  =  x'x  +  y'y  +  hz .  (1.3) 

The  unknown  and  known  linear  current  densities  are  defined  by 

J'(r')  =  n'xH'(r')  ,  Jinc(r')  =  fl'xHiK(r')  (1.4) 


respectively,  where  H'  and  H",c  are  the  total  (unknown)  and  incident  (known)  magnetic  fields. 
A  harmonic  time  dependence  exp (iaX)  is  suppressed  throughout. 

To  the  integral  equation  in  (1.1),  we  now  apply  the  RWG  (reference  4)  formulation  of  the  MoM. 
This  is  a  well-known  procedure  and,  hence,  we  provide  only  a  brief  outline.  The  scatterer’s 
surface  is  converted  to  one  made  up  of  flat  triangles.  The  total  number  of  edges  is  N.  We 
consider  the  edge  as  a  directed-line  segment  and  we  define  identical  basis  and  test  functions 


f  (r)  = 


24 

_L_ 

2  a: 


P  n  ’ 

P  n  ’ 


r  eC 

elsewhere 


(1.5) 


where  ln  is  the  length  of  the  edge  and  4)  ( An  )  is  the  area  of  the  triangle  attached  to  the  edge  and 

positively  (negatively)  oriented  with  respect  to  the  normal.  The  vectors  p  and  p  lie  on  the  two 

triangles,  the  first  emanating  from  the  comer  opposite  to  the  side  under  consideration  and  the 
second  terminating  there. 
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We  expand  the  unknown  current  density  in  terms  of  the  basis  functions 

N 

J'(r  )  =  2>A(r> 

n= 1 

and  substitute  in  the  integral  equation  (1.1) 

=  J'”(r  ')-/>'  x  f >  J  f„(r)  x  Vg(r,  r  ’)dS ,  r’  £  S 

^  n= 1  «=1  f 

where 


T  =T+uT\ 

n  n  n 

We  test  this  equation  by  taking  the  dot  product  with  f  n(r')  and  integrating  over  S 

|  f,„(r'K(r')<«'  =  j  r(v')dS' 

^  n=  1  Tm  Tm 

f  W  •  «'x  f  f„(r)  X  Vg(r,  r')dS \dS' ,  m  =  1,2,...,  A  . 

«=1  Tm  {  Tn  \ 


(1.6) 


(1.7) 


(1.8) 


(1.9) 


The  only  values  of  n  in  the  first  integral  in  this  expression  that  contribute  are  those  for  which  f 
has  a  domain  in  common  with  fm .  Using  standard  vector  identities  for  the  last  integral,  we  can 
write  this  expression  in  the  fonn 

J  =  j  f>')- r(r')dS' 

^  /j=1  Tm 

[fffl(r')x«']-ffH(r)xVg(r,r')^ \dS' ,  m  =  1,2,...,  N .  (1.10) 

«=1  Tm  [  Tn  J 

We  are  interested  in  the  accurate  calculation  of  the  integral 

I,(r')  =  |f,(r)xVg(r,rVS.  (1.11) 

7 ; 

For  simplicity,  we  consider  a  single  triangle  T7,  as  shown  in  Figure  1.1,  and  define 

I(,) (r ')  =  J f (,) (r) x  Vg(r , r')dS ,  f(,|(r)  =  r-r,  ,  /  =  1,2,3.  (1.12) 
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The  function  f(/)  is  proportional  to  a  basis  function.  The  superscript  refers  to  a  corner  of  T. 
Coordinates  are  centered  at  the  centroid  of  the  triangle  with  the  z-axis  perpendicular  to  the 
triangle. 


Figure  1.1:  The  Integration  Triangle 

The  integral  in  (1.12)  is  the  one  we  are  interested  in  computing  to  a  prescribed  number  of 
significant  digits.  To  this  end,  we  transform  it  to  a  form  that  we  can  evaluate.  We  write 


A(r,r')  =  (r-r,)xVg(r,r')  =  (r-r')xVg(r,r')  +  (r'-r/)xVg(r,r') 

Qq(y  r') 

=  (r-r')xVg(r,r')  +  (r'-r/)x  V^g(r,r')  +  z  ;  (1.13) 

oz 


where  the  subscript  S  denotes  the  surface  gradient.  From  (1.2) 

r  ro  =  ylr - r'l  =  ag(r’r')  r~r' 

S|r  — r'|  9|r  — r'|  |r-r' 

and,  thus,  the  first  of  the  two  terms  in  (1.15)  is  equal  to  zero.  We  can  then  write 


(1.14) 


(1.15) 
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From  ( 1 . 12)  we  can  then  write 

I(,)  (r ' ')  =  | (r  -  r, )  x  Vg(r, r')dS  =  (r'  -  r, )  x  | J Vsg( r,  r')dS  +  zj^^-  ^  dS 
=  (r  1 ' - r, ) x  J  vg(r, r')ds  - z x (r'  - r, )  J  dg^,r  ^ dS 

dT  T  OZ 


(1.16) 


where  v  the  exterior  unit  normal  to  the  side  of  the  triangle  (lying  in  the  triangle’s  plane)  and  dT 
is  the  triangle’s  boundary.  We  have  also  used  Identity  43  on  p.  503  of  reference  8  in 
transfonning  the  first  surface  integral  into  a  line  integral.  We  deal  with  the  remaining  surface 
integral  in  the  next  section. 

We  have  thus  expressed  the  surface  integral  in  (1.12)  in  terms  of  a  line  integral  and  a  surface 
integral.  This  is  the  principal  result.  We  concentrate  next  on  the  case  where  r'  is  a  point  in  the 
triangle’s  interior.  In  this  case,  the  first  term  in  (1.16)  is  a  vector  that  has  only  a  z-component. 
When  that  vector  is  dotted  into  the  vector  f„,(r')  x  n'  in  (1.10)  the  result  is  zero  since  the  latter  is 

transverse  to  this  direction.  In  order  to  evaluate  the  contribution  of  the  second  tenn  in  (1.16),  we 
move  the  OP  to  the  position 

f'  =  r  '  +  ez,  s>0;  r'eT  (1.17) 


we  evaluate  the  integral  and  then  let  ego  to  zero.  For  the  integral,  we  write 


•dg(r,r') 


dz 


dS  =  — 


4x 


J  —  |r  -ez\  *  ^5*+  J 


dz 


d_ 

dz 


—ik\r—€z\  i 

e  1  1  -1 


r  -sz\ 


dS 


(U8) 


The  first  of  these  integrals,  J\,  is  the  integral  that  measures  the  solid  angle  subtended  by  the 
triangle  at  the  point  V'  (see  Appendix  C).  As  s  — >  0 ,  this  integral  goes  to  In 


J,  (r')  =  2x . 


(1.19) 


This  contribution,  however,  has  already  been  accounted  for  in  the  derivation  of  the  MFIE  (1.1). 
Thus,  we  must  omit  it  here  or  we  will  have  a  double  count. 

For  the  remaining  integral  in  (1.18),  we  let 

R  =  \r-ez\  (1.20) 
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and  write 


A(r')  =  J- 


rA 

e  “  -f 

rfS=j 

z; 

rA 

e  lkR  -\ 

dS  + 

T- 

r  A 

-QikR-\ 

1 

■  i  f~~< 

_  R  \ 

dz 

0 

R 

DUSZ 

R 

dS 


(1.21) 


where  Do  is  a  solid  disk  with  center  at  the  origin  and  radius  p0 .  The  radius  is  small  enough  so 

that  the  entire  disk  lies  within  the  triangle.  The  second  of  these  integrals  is  over  the  region 
bounded  by  the  triangle’s  perimeter  and  the  disk’s  circumference;  thus,  this  integral  is  not 
singular  when  s  —  0  and,  since  its  integrand  is  equal  to  zero  when  s-  0 ,  the  integral  itself  is 
equal  to  zero.  For  the  first  integral,  we  write 


d 

QikR  -  1 

_  d 

a 

l 

£ 

1 

_ i 

2  -f) 

's'' 

dz 

R 

dR 

R 

[r  ) 

[r  ) 

JkQm  ,  e  mi2  sin  (AA/2) 

S  .  l  .  KQ 

R2  R  kR/2 


from  which  we  get  that 


-ikRe~m  -(e~m 
R2 


(1.22) 


d 

"  Q~m  -  1 1 

dz 

R 

(1.23) 


Thus, 


dz 


e  - 1 
R 


dS 


< 


Pa 

4 /rkzj  - 


pdp 
2  2 
o  P  +£ 


=  Inks 


ln(/v+g2)-21n(g)]  ^0->0.  (1.24) 


Thus,  (1.21)  is  equal  to  zero  in  the  limit.  From  this  result  and  (1.21),  we  have  for  the  limit  of 
(1.18) that 


•  <3e(r,r')  1 

SV  ’  dS  =  — ,  ref. 

dz  2 


(1.25) 


As  we  pointed  out  above,  this  contribution  has  already  been  accounted  for  in  the  derivation  of 
(1.1).  For  computational  purposes,  this  integral  should  be  set  equal  to  zero  whenever  the  OP  is 
in  the  interior  of  the  triangle.  Consequently,  the  entire  term  in  (1.16)  should  be  set  equal  to  zero 
in  such  a  case 

I(/)(r')  =  J(r-r/)xVg(r,ryS  =  0,  r'eT.  (1.26) 

T 
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We  close  this  section  by  noting  that  an  alternate  testing  of  (1.7)  leads  to  the  same  inner  integral 
as  the  one  above.  We  show  this  in  Appendix  A. 
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SECTION  2:  REMAINING  SURFACE  INTEGRAL 
For  the  last  integral  in  (1 . 16),  and  for  all  r '  &.T  (in  fact,  not  in  the  triangle’s  plane),  we  write 


FI 


(r-)=fM^rfs  =  __Lf£J£_  rfs  =  _£_f£. 

J  r)v  A  it  J  ‘P\r7  /?  A  it  J 


1  r  8  f  e 


-ikR  \  „-ikr'  a  (  -ik(S-r')  A 


dz 


An  JT  dz  {  R 


e'!V  | 

r a 

cos  [&  (R -/)]-?  sin  [&(R-r')]j 

An  * 

[dz 

R  J 

An  *  dz 


dS 


v  *  j 


dS 


-ikr'  N'~]  (-1)"  k2" 


Z 


An  n=0  (27?)!  '  dz 


f  — J 

\(R-rf'} 

'  ] 

R 

ikr'  Ns-l  ( _\  \"  J-2n+l 

>dS  +  - — £- 


w 


(-i  ye 


An  ^  (2/1  +  1)!  { dz 


f  — J 

\R-r')2n+1' 

r  ] 

R 

dS.  (2.1) 


Here  we  assume  that  r'  g  T  for  then,  by  agreement,  we  set  F  equal  to  zero.  F  is  also 
(legitimately)  zero  whenever  the  OP  is  in  the  triangle’s  plane  and  exterior  to  the  triangle.  We 
nonnalize  all  distances  with  respect  to  wavelength  (see  also  reference  1,  p.  5)  and  write 


-z'2 nr'  N-\ 


(~l)\ln)2n  r  d  \{R-r')ln 


F(  r')  =  -- - 

y  ’  An  to  (2 n)\  Jr 


dz 


R 


dS 


le 


-ilnr'  AG-1 


An 


Z 

n=0 


(-1)'  M  f  B 

(277  +  1)!  \dz 


riL 

(R-rf'+r 

1 

1  <0 

R 

dS 


(2.2) 


with  distance  measured  in  wavelengths. 
We  define 


Ur')  =  J 


rJL 

> 

1 

■! 

f 

1 

— » 

R 

dS ,  77=0,1,2,. 


and 


Q„(r')  =  {n  +  l)fV  r  ’  dS,  77=0,1,2,... 


(2.3) 


(2.4) 


then 


f  R-r'^ 

J 


>dS 


■  d(R-r' )" 

dz 


dS  -  r'Pn  (r()  =  -nz'j  ^  dS  -  r'Pn  (r') 


R 


(2.5) 
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from  which  we  get  that 
/J(r')  =  -r7>(r') 


(2.6) 


and 


Pn+2  (r)  =  ~z'Qn  (r)  _  r'P„ 


n+ 1 


(r'),  n  =  0,1,2,... 


(2.7) 


This  iteration  converges  for  r'  <  1 .  Since,  however,  Qn  is  computed  also,  then  the  OP  must  be 
interior  to  a  sphere  centered  at  the  centroid  and  radius  of  one-half  of  a  wavelength  (reference  1, 
Section  3).  We  also  note  that  (2.4)  has  been  thoroughly  studied  and  evaluated  in  (reference  1, 
Section  2  and  Appendix  A).  It  remains  then  to  evaluate  (2.3)  for  n  —  0 .  We  will  do  that  in 
Section  4. 


While  performing  verification  and  validation  work  on  (2.2),  we  noticed  that,  when  the  point  z'  is 
over  the  triangle  and  tending  to  zero,  the  imaginary  part  of  F  does  not  tend  to  zero  as  it  should 
but  to  a  finite  albeit  small  number.  We  traced  the  cause  of  the  problem  to  the  fact  that,  in  the 
final  stage,  we  subtract  two  almost  equal  numbers.  We  modified  the  algorithm  for  this  case,  and 
we  present  the  details  in  Appendix  B.  We  summarize  the  result  as  follows. 

For  a  point  over  a  triangle,  if  either 


,  V 5  •  icrM 

z  < - 

2  n 


(2.8) 


or 


VlOg 

In 


(2.9) 


then 


N-2 


Im{F(r')}  =  z'X 

n= 0 


(-l)"(n  +  l)(2.r 


(2n  +  3)! 


-\p2"dS. 


(2.10) 


Above,  M  is  the  number  of  significant  digits  desired  and  s  is  the  machine  epsilon.  The  latter  is 
given  in  Table  B.  1  for  different  precisions.  Moreover, 


N  =  max  ( Nc ,  Ns ) 


(2.11) 


with  Nc  and  Ns  obtained  from  (3.14)  and  (3.12),  respectively. 

We  proceed  in  the  next  section  with  the  evaluation  of  the  first  integral  in  (1 . 16). 
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SECTION  3:  NUMERICAL  CONSIDERATIONS  FOR  LINE  INTEGRAL  IN  (1.16) 


In  this  section,  we  evaluate  the  integral  in  (1.16) 


>  (r ' ')  =  (r ' '  -  r, )  x  J  vg(r,r')ds  =  (r'  -  r, )  x  £ v}  |  g(r,  r')ds 


where  v,  is  the  exterior  unit  nonnal  on  the  j- th  side  of  the  triangle,  lying  in  the  triangle’s  plane, 
and  Sj  is  the  j- th  side  and,  also,  its  length.  We  can  also  write  this  result  in  terms  of  its  component 
in  the  triangle’s  plane  and  one  along  the  z-direction 


if 1  (r ' ')  =  (p'  +  z  z  -  r, )  x  Y,  Vj  {  g(X,  r  ')ds 

=  (p'  -  ri ) x  X ?j  {  g(r>  r')ds  +  z'z  x  X v  J  g(r .  r')ds 

j  1  Sj  J  1  Sj 

=  -z(p'  -  r7 )  •  Yjj  J  g(U  r')ds  +  z'Yjj  {  g(r,r')ds 


where  t .  is  the  unit  tangent  to  the  j- side  of  the  triangle,  given  by 


t j  =  z  X  V j 


p'  =  x'x  +  y'y 


is  the  projection  of  the  observation  point  onto  the  xy -plane. 

We  consider  the  integral 

— /Ar|r— r'|  -ik(R-r')  -i2x(R-r') 

f  ? - rds  =  e ~ikr'  f - - ds  =  o~i2nr’  f - - ds  =  e  ,2rcr'  I  (r',j) 

J  r_r'  Jr  Jr  v  J} 


where,  in  the  last  integral  all  distances  are  measured  in  wavelengths.  The  multiplication  and 
division  by  the  exponential  is  justified  below.  We  re-write  the  integral  as 


=  jCOs[2^(R-r')]-/sin[2^(R-r')]^ 
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In  evaluating  this  integral,  we  follow  the  procedure  in  reference  1.  We  can  show  that  the 
argument  of  the  trigonometric  functions  is  bounded.  The  centroid  of  a  triangle  divides  each 
median  in  a  2:1  ratio.  The  longer  of  the  two  parts  is  the  distance  from  the  centroid  to  the  vertex 
of  the  median.  Thus,  the  largest  distance  from  the  centroid  to  a  point  on  the  triangle  is  equal  to 
two  thirds  the  longest  of  the  three  medians.  However,  a  median  is  shorter  than  the  longer  of  the 
two  adjacent  sides;  thus,  the  longest  distance  from  the  centroid  to  the  triangle’s  perimeter  is  less 
than  two  thirds  the  length  of  the  triangle’s  longest  side.  We  combine  this  with  the  fact  that  the 
side  of  a  triangle  (in  this  case  that  formed  by  the  observation,  integration  and  centroid  points)  is 
always  greater  than  the  difference  of  the  other  two  sides  to  obtain  the  bound 

(3.7) 


where  /max  is  the  length  of  the  longest  side  of  the  triangle.  Though  this  is  not  a  strict  bound,  it  is 

one  that  holds  for  all  observation  points;  thus,  no  matter  how  far  away  the  observation  point  is 
from  the  integration  triangle,  the  number  of  terms  in  the  expansion  will  be  the  same  as  for  a  point 
near  the  triangle.  Moreover,  for  an  actual  grid,  we  can  search  among  all  triangles  for  the  longest 
side  and  use  that  value  in  (3.7).  In  this  way,  we  do  not  have  to  test  all  triangles  separately, 
conserving  a  lot  of  compute  time.  Had  the  exponential  factor  not  been  inserted  in  (3.5),  we 
would  need  more  and  more  terms  in  the  expansion  as  the  observation  point  receded  from  the 
triangle. 

We  proceed  now  to  expand  (3.6)  in  a  Maclaurin  series  and  keep  the  first  N terms 


N- 1 


Inc,ns  (r'J)  =  X 


(-l)"(2*f 


(2«) ! 


f 


N-\ 


(-ir(2*r  f(R~r') 


R 


-  ds  -  i 


\2k+1 


n=  0 


(2/7  + 1)! 


S' 


R 


-ds 


(3.8) 


In  detennining  the  number  of  terms  N,  we  use  Taylor’s  Theorem  with  a  Remainder  (reference  5, 
p.  113).  If  we  require  a  number  M  of  correct  significant  digits,  then  we  proceed  as  follows  to 
determine  the  number  of  terms.  The  series  for  the  sine  is 


sin  2/r(R-r') 
2k  (R  -  r') 


h  (2;.  +  !)! 


+  R,(N,2x(R~r')) 


(3.9) 


where  Rs  is  the  remainder  of  the  series  and  is  bounded  by  the  first  omitted  tenn 


R,(N,2x(R-r')) 


< 


\2x\R-r'\  1” 

L  I  Imax  J  ^ 


^  Aft  ^ 

3fma: 


2  N 


(2JV  +  1)!  (2JV  +  1)! 


(3.10) 
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We  can  then  write  for  the  relative  error 


sin[2;r(7?- 

r')]  ^(-l)"[2;r(7? -/)]-" 

2  n(R-r')  fa  (2n  +  l)\ 

sin  2n(yR-r'^) 

2n[R-r ') 

R\U,2n(R-r')) 

/ 

An  i 

rs  ^max 

.  ^  y 

2N 

l 

sin  2n(yR-r'} 
2n{R-r') 

(27V  +  1)! 

sin  2 n 

R-r’  I 

max  J 

2  n\l 

l  —  r'\ 

Imax 

An 


,2  N 


(27V +  l)! 


sin 


An 


v  3 
An 


L 


\ 

J 


(3-11) 


The  number  of  terms  is  determined  by  requiring  that  the  relative  error  is  smaller  than  5  divided 
by  the  number  10  raised  to  the  number  of  significant  digits  Mplus  one,  or 


f 

v 


2  N 


An  j  3 


(27V  +  1)! 


<5-irM"' 


sin 


An 


An 


y 


(3-12) 


Once  the  triangle’s  dimensions  are  known,  we  can  solve  this  expression  for  N. 

In  this  discussion,  we  have  assumed  that  the  argument  of  the  sine  function  is  small;  in  practice, 
the  longest  side  of  a  triangle  does  not  exceed  one  tenth  of  a  wavelength  and,  hence,  the  largest 
argument  is  approximately  equal  to  0.4yr/3  of  a  radian  or  24  deg.  This  means  that  the 
trigonometric  functions  in  the  original  integrals  hardly  exhibit  an  oscillatory  behavior. 

For  the  cosine  we  have 


cos  [27r(R  ^  ^  ^ 

[2  *(R~r')Y 

(4  ,  V 

L  1  fa  (2«)! 

.<  (“)!  < 

cos 


[2n(R-r')\ 


cos 


\ln(R  r')max]  (2A^, 


cos 


-nl 


(3.13) 
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and  the  test  is 


'4  h 

~  x 

W  J 


2  N 


(2A)! 


(4  3 

<5-10  M  'cos  — 7ll„ 


V- 


J 


(3.14) 


As  an  example,  in  Table  3.1  we  give  the  number  of  terms  in  the  sine  and  cosine  expansions  that 
are  required  to  guarantee  a  certain  number  of  significant  digits. 


Table  3.1:  Number  of  terms  that  guarantee  a  given  number  of  SD. 
Largest  side  of  triangle  less  than  or  equal  to  a  tenth  of  a  wavelength. 


No.  of 
SD 

Sine 

Cosine 

4 

3 

3 

5 

3 

3 

6 

4 

3 

7 

4 

4 

8 

4 

4 

9 

5 

5 

10 

5 

5 

11 

5 

5 

12 

6 

6 

13 

6 

6 

14 

6 

6 

15 

7 

7 

We  now  let 


J 


(«) 


M=J 


R 


ds 


n  =  0,1,2,... 


(3.15) 


We  can  now  proceed  as  with  (2.5).  We  do  not  see  the  advantage;  however,  over  the  following 
procedure.  We  expand  the  binomial  in  (3.15)  and  write 


(3.16) 


while 


m= 0 


n+m  I  fl 

\m 


J  R'"'ds 

sj 
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^-lyvv'VhZH) 


«+'«  n  in¬ 
i’ 

m 


'  J R"'lds  ,  n  =  1,2,... 


(3.17) 


We  separate  even  and  odd  powers 


jW(r')  =  (-l)V'Vf(r')  +  (-l)"X 


n  X  1  /  /7  1  fn-lm 


J  tf2'"  'A 


[(»-l)/2] 


(-1)”  I  1 2/n?+ 1 1 r 2”'  1 J  Rlmds 


(3.18) 


In  order  to  evaluate  these  integrals,  we  introduce  the  following  notation  (see  Figures  3.1  and  3.2) 


r  =  rj+l+stj,  j  =  1,2,3 


(3.19) 


with  the  indices  running  cyclically.  We  also  write 


r'  =  x'x  +  y'y  +  hz-p'  +  hz ,  rl  -  x,x  +  y,y 


(3.20) 


ay+1  =r.+1-p',  T  =  s  +  tj- a.+1,  A.  =  |fy  x(f.  xay+1)| ,  cr  =  |r-p'|  =  |a .+1 +^.| .  (3.21) 


From  these  we  can  write 


-  s~  +  2*/  '  a/+ls  +  a  /+l  -  (5  +  0  •  a;+ 1  )  +  a/+i  -  ((,'  •  a,'+i  )  -  (s  +  ij  •  a;'+ 1 )  +  (3  -22) 


We  also  define 


b  =r  -r';  B=AA2+h2,  j  =  1,2,3 . 


(3.23) 
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We  distinguish  three  cases.  For  t .  •  b  ,+2  >  t  ,  ■  b  +1  >  0 ,  we  make  the  transformation 
r  =  Bj  tan  0 ,  dr  =  2F  sec2  OdO 


(3.25) 


and  write 


y-0)  (r'j  =  J  sec  OdO  =  J 


'  'sec<9(sec6,+  tan<9) 


=  In  sec  6  +  tan  0\ 


arctan 


arctan 


tra;+ 1 


*  1  {  =  In 


sec  0  +  tan  9 


dd 


v  y 


1  + 


tj  ~by+2 

v  *2  y 


A2 


+ 


tj  ^  ;+2 


B. 


1  + 


tj  'by+i 
v  y 


+ 


{J  'by+i 

5. 


=  In 


r~. — ■ 
Bj+\ 

|  +ty'by+2 

\ 

Ib2+  1 

(6  by*0 

1  +Vby+1 

=  In 


b./+2 

+  Vb;+ 2 

by+1 

+  VVi 

We  examine  next  the  case  where  tj  ■  b  ,+1  <  t  j  ■  b  /+2  <  0  and  write 


Jl")(r')=[  —  =  [  ,  r^lnlsecff  +  tanffl 

2  V  '  J  P  J  _2  .  n  2  1  ! 


7?  .  J  It2  +  B  2 


arctan 

h  *  by'+l 

J 

arctan 

h  ‘  b/+2 

5 

J 

=  In 


1  + 


h  -b./+1 

v  B,  y 


+ 


(J  "by+l 


5. 


ti '  by+2 

v  y 


1  + 


+ 


fj  'by+2 


5 


=  In 


b/+i 

+ 

0  '  b/+, 

by+2 

+ 

tj '  b  /+2 

For  the  case 


tj  'b/+i  <0,  t.  •  b/+2  >  0 


(3.26) 


(3.27) 


(3.28) 
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we  write 


rds  '  | 

'y+i| 

r  _ 

dr 

•IT” 

si  ( 

i 

1  V 

It2+B2 

dr 


=  ln|sec#  +  tan<9| 


arctan 

h  ‘V 

arctan 

B 

+  ln  sec  #  + tan# 

;  5 

j 

1  1 

V  v  9 

0 

0 

=  In 


=  In 


1  + 

ft  b  > 

lj  uj+ 1 

2 

+ 

0  ‘by+i 

+  In 

J 

c 

ff  -b  "1 

lj  Vj+2 

2 

,  Vb2+2 

i  B.  1 

B 

V 

1  5.  I 

B 

V  J  9 

J 

V 

V  29 

J 

bi+i  +  h  -b/+l 

+  ln 

by+2  +Vb/+ 2 

5 

5 

j 

j 

(3.29) 


We  next  evaluate  the  integral 
(r')  =  J  Rds 


(3.30) 


in  each  of  the  above  three  regions.  In  the  first  region  (f;.  -b,+2  >  f.  -b,+1  >  0)  we  use  (3.25)  and 
write 


- yj  -j+Zj 

K(j\r')  =  Bj2  J  sec3  Odd  =  B  f  tan 0 sect? 


arctan 


arctan 


>;B’  <~B,2 

tl  a/+ 1 

V  BJ  J 


J  tan2  0  sec  Od 0 


-Bf\ 


ft  ,a  A 

tj  a./+ 2 
v  BJ  J 


(7  n  V 


1+ 


t  a  , 

j  j+- 

v  B,  , 


Y  -a  ^ 

f  a,/+i 

v  Bj  J 


1  + 


ft  -a  A 

h  a./+i 
v  bj  j 


' i'*i +2 


w|1,(r')  +  5/  j  sec  Odd 


/V 1 


(0  ‘  a,/+2 )  +  (tj  •  a,+2 )  (/ 


,  ViJV5/ 


+  (7-  a;hl  j  -  Yr>  (r')  +  ZY  J  sec  Odd 


(3.31) 


or 


(3.32) 
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with  J{p  (r')  given  by  (3.26). 


In  the  second  region  ( t  .  ■  b  /+l  <  £.  ■  b  .+2  <  0 ) 


arctan  I  A- -a 


r*j+2\) 


Kf( r')=  J  ^t2  +  B2dr  =  B2  J  sec3  Odd 


\'rai+ 2 


r/ay+i|J 


1  riA  ii  i  i/v  ii  i“i  B 


2l 


l^-b J+1  b/+1  -  tj- b.+2  b.+2| 


+^A0)  (>•')• 


In  the  third  region  ( •  b  ;+1  <  0  ,  tJ  ■  b  /+2  >  0  ) 


1 2  y+1i  _  y  y+2  _ 

(r')  =  J  >/r2+£/.2rfr  +  J  ^ r+B/dr 


=  B,  tan 6 sec 6 


arctan 


B. 


y^y+i 


-B, 


|  tan2#  sec  6d6 


+B  r  tan  6  sec  6 


arctan 


■  a  ^ 

a;+2 

v  B,  j 


-B2  J  tan2  9 sec 6d6 


=  B2  \ 


Va/+> 


B 


1  + 


-a  A 

l)  a/+i 

v  B,  , 


+ 


n  .  a  y 

a./+ 2 

v  BJ  J 


1  + 


/  -  \2 

v  , 


-KP(r>) 


+B2  J  sec  OdO  +  B2  J  sec  Odd 


or 


lr 


A1  (r') = jL  o  •  v  I  I + (t  ■  b/+- )  lb/« 


+^-yoV) 


(3.33) 


(3.34) 


(3.35) 
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We  next  evaluate  the  integral  in  the  first  sum  in  (3.18).  In  the  first  region 


2m-l 


iffm+1)  (V)  =  J R2m+lds  =  J  R2m-l(r2+Bi2)dT=  J  (t2+B2)  2  T2dT  +  B2Kfm~1\ r') 


=  — - — (r2  +B2) 
2m  +  P  J  ’ 


2m+l 


tj-Zj+2 


1 


V'a7+ 1 

K{2m+1)(r')  +  B2K{2m-l](r') 
fi-a/+I  2m  + 1  7  '  K  ’ 


(3.36) 


from  which 


K(lm+l)  (r')  =  — - — 
1  v  '  2m +2 


ft  '  a;+2  )|b.,+2  |2,"+1  -  ft  ‘  a;+l  )|b 


+  2m+j_  2  ,2m  1)(r  m  =  12  (3.37) 

2m +2  33  v  ’ 


Similarly,  in  the  second  region 


1  (j  .a  b  2m+1 

2m +  2  v 7+1  7+1 


iff  m+1)  (r')  = 


+  m  =  1,2,...  (3.38) 


while  in  the  third  region 


if 


(2/n+l) 


= 2^rv2  (6- '  )lb>- 1’”"  -  lfj  ■ a^*  llb>*'  r""‘ )  -  <■")  • 


(3.39) 


For  the  integral  in  the  second  summation  in  (3.18),  we  have 
K{p{y')  =  \ds  =  sJ 


(3.40) 


for  all  regions.  In  the  first  region 


iff ")  (r')  =  \R2mds  =  J  R2m-2[z2  +B2)dT=  J  (z2  +  B  2J~\2dz  +  B2K(2m-2)  {r') 

sj  0'ay+ 1  f/‘a/+i 

=  —  (r2  +B2)m  t  h'*J+2  — —K{2m)(r')  +  B2K(2m~2)(r') 

2m  \  J  /  /  • «»  J  \  f  J  J  V/ 


y  j+ 2 , 

ij '  a./+ 1  2/h  7 


(3.41) 


so  that 


jtf2")(r')  =  — L_ 

7  '  ;  2m +  1 


ft  -a,+2)K+2r"'  -ft  ’  ay+i )  |by+i  r  +-^-B2Kf^(r'),  m  =  1,2, 


J  2m  + 1 


(3.42) 
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In  the  second  region 


K{2m)(  r')  =  — - — 
1  v  ’  2m  +  \ 


/v  2m  /\  2m 

tj  '  aj+l  h/.l  /,  'aM  b./+2 


and  in  the  third  region 


This  concludes  the  evaluation  of  the  integral  in  (3.6)  by  its  approximation  in  (3.8). 


(3.43) 


(3.44) 
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SECTION  4:  NUMERICAL  CONSIDERATIONS:  THE  SOLID  ANGLE  INTEGRAL 


We  turn  to  (2.3)  for  n  —  0 


(4.1) 


As  we  point  out  in  Appendix  C,  this  expression  represents  the  solid  angle  subtended  by  the 
triangle  at  an  observation  point  that  lies  in  the  half  space  in  which  the  z-axis  points.  If  it  lies  in 
the  other  half  of  the  space,  then  (4.1)  yields  the  negative  of  the  solid  angle.  We  thus  can  write 


r' -z  >  0 
r'  -z <  0 


(4.2) 


where  Q(r')  >  0  is  the  measure  of  solid  angle. 


For  z'  =  0,  and  for  r'  entirely  outside  the  closure  of  T,  (4.1)  is  equal  to  zero,  as  a  simple 
evaluation  shows.  If  r'  is  in  the  interior  of  the  triangle,  then  the  integral  is  equal  to  +2/T  as  r' 
approaches  the  triangle.  This  contribution,  however,  has  been  accounted  for  in  the  derivation  of 
the  MFIE  (2.1);  in  fact,  it  is  the  reason  for  the  Vi  factor  in  the  left-hand  side  of  it.  Thus,  for  r'  in 
the  interior  of  the  triangle,  we  must  treat  the  integral  in  (4. 1 )  as  an  improper  integral,  in  which 
case  it  is  equal  to  zero.  This  is  so  because  the  solid  angle  of  the  triangle  is  2n  and  so  is  the  solid 
angle  subtended  by  the  small  circle  with  center  the  OP.  The  net  solid  angle  is  the  difference  of 
the  two,  which  is  zero.  The  case  where  r'  is  a  point  on  the  boundary  of  T  does  not  occur  here. 
We  also  point  out  that,  for  the  same  reasons,  the  function  F  in  (2.1)  is  equal  to  zero  for  z'  =  0 .  In 
fact,  if  we  examine  the  first  vector  in  (1.16),  we  see  that  it  points  in  the  z-direction  so  that  its  dot 
product  with  the  test  function  (as  in  (1.10))  is  equal  to  zero  since,  in  this  case,  the  test  function 
lies  on  the  integration  triangle  and  is  orthogonal  to  the  z-axis. 

One  formula  for  computing  the  solid  angle  subtended  by  a  triangle  is  reference  9 


Q(r')  =  2  tan  1 


_ MvM _ 

b\b2b3  +(bj  -b2)h3  +(b2  •  b3  +(b3  -bj)b2 


(4.3) 


where  the  vectors  are  given  by  (3.23)  and  the  scalars  indicate  the  magnitude  of  the  vectors. 
Another  formula  is  given  in  reference  10  but  it  requires  more  operations  and  is  to  be  used  only  if 
we  notice  numerical  instabilities  in  the  above.  This  fonnula  is 
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Q(r')  =  2^  tan  1 


<r(4x4+i)fl ' 


1/2 


bxb 


n+ 1 


_ yn+ 1 

v 1 + 4  ■  4+i  / 


i-r^  + 


<y 


4  x(4  x4+i) 

i 

A-> 

+ 

b  xb , 

n  n+ 1 

U+4-4+J 

,  <7z  -  0 ,  4=^7  (4-4> 

b„ 


where  g  is  a  unit  vector  with  its  tail  at  the  OP  and  such  that  its  extension  does  not  cross  the 
triangle.  If  it  does,  then  4 /r  must  be  added  to  the  result  (reference  11). 

We  point  out  that  the  surface  integral  in  (4.1)  can  be  converted  to  a  line  integral  and  the  same  is 
true  for  the  integral  of  the  normal  derivative  of  the  free-space  Green’s  function  (references  11 
and  12). 

Since  (4.3)  returns  a  value  of  2n  when  the  OP  is  a  point  on  the  triangle’s  interior,  and  this  has 
already  been  accounted  for  as  explained  above,  we  must  also  subtract  2n  from  it  to  bring  it  back 
to  zero.  From  the  discussion  above,  this  amounts  to  setting  (4.1)  equal  to  zero  whenever  the  OP 
is  in  the  same  plane  as  the  triangle. 
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SECTION  5:  NUMERICAL  CONSIDERATIONS:  Qn(r')  IN  (2.4) 


As  mentioned  at  the  end  of  Section  2,  (2n(r')  in  (2.4)  has  been  evaluated  in  (reference  1,  Section 

2  and  Appendix  A).  We  quote  here  results  from  that  document.  The  relevant  integral  appears  in 
(1.21)  of  reference  1  and,  as  a  line  integral,  in  (1.22).  We  quote  the  latter  here 


0„(r')  =  Xv(r/+i 


j= i 


(5.1) 


where  the  notation  is  that  of  Section  3  (Equations  (3.19)  -  (3.23)  and  Figures  3.1  and  3.2.  We 
further  write  this  as 


Q,(r')  =  Zv(r,-+i  -?')Kn,j(r') 

7=1 


(5.2) 


where 


K. 


».y(r')  =  j 


n  =  0,1,2,... 


(5.3) 


From  this 


=  I 


yj p~  +  h 1  - 1  h 


ds 


P 


K\j  =  Sj  -  2r'K0  j 


Kn+lJ=V„_UJ-2r’KnJ-p'2Kn_UJ,  n  =  1,2, 


where 


V„ 


,(r  ')  =  \Up2+h2 -r’\+l  ds,  n  =  0,1,2 


(5.4) 

(5.5) 

(5.6) 


(5.7) 


5.1.  RESULTS  FOR  Vn  j 

We  have  three  distinct  regions  of  space.  In  the  region  t.  ■  b  ,+2  >  t  •  b/+1  >  0 ,  we  have  from  (2.21) 
of  reference  1  that 
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n+ 1 

v . =y 

n,J 


m= 0 


(«  +  l)! 


m\(n  +  \  —  my. 


r') 


n+ 1- 


vm  ;. ,  n  =  0,1,2,. 


(5.8) 


where 


V0  J=SJ 


vu  2  < 


bi+2|(0  -by+2)-|by+I|(f,  b J+l)  +  Br  In 


b/+2 

+  *j'hJ+2 

Vi 

+  *j'hJ+l 

(5.9) 

(5.10) 


and 


"uj  m  + 1 


77 {K4"  (tj -b,+2)-|by+ir (ij  -by+1)}  + 


m  + 1 


Bj\-2j »  m=2,3„ 


(5.11) 


In  the  region  f .  •  b /+1  <  ?.  •  b /+2  <  0 ,  we  still  have  the  expansion  (5.8)  with  (5.9)  holding,  while 


vu  2  ^ 


b/+i  tj  •  b/+l  -b7+2  tj  •  b y+2  +7  In 


b  /+i  +  tj  ■  b  /+1 


b/+2  +  tj  -  b/+2 


(5.12) 


and 


'”'i  111  +  1 


b/+i  *,•  •  b/+i  -  b/+2  tj  -by+2l!-  + 


/H  +  1 


>  m  =  2,3,. 


(5.13) 


In  the  region  f .  •  b  /+1  <  0 ,  1 1  •  b  /+2  >  0 ,  we  still  have  (5.8)  and  (5.9),  while 


*7  j+2 


vu  2  < 


by+i  1 1  tj  •  b  /+,  |  +  |b  J+2 1  (ij  •  b  J+2 )  +  Br  In 


by+l  +  h 


•by+1|)(|by,2|  +  fy-b,+2) 


B 


(5.14) 


and 


v  .  = - 

m,j  ,  i 

m  + 1 


V+i 


tj  ‘hj+i 


+  b 


j+2 


t ,  b 


j+2 


+ 


m 


m  + 1 


B  2v  ,  . 

j  m-2,j 


m  =  2,3,... 


(5.15) 


This  concludes  the  evaluation  of  (5.7).  The  iteration  formulas  (5.11),  (5.13),  and  (5.15)  do  not 
converge  everywhere.  In  Appendix  B  of  reference  1,  we  show  that,  for  the  iteration  to  be  stable, 
we  need  0  <  B t  <  1 . 
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5.2.  RESULTS  OF  K0  J 


In  the  first  region,  we  have  from  (2.57)  in  reference  1  that 


K0J  =  In 


by +2 

+ 

tj  *  by +2 

> 

h 

4+i 

+ 

44+1 

7 

1  4 

tan 


h  1 

'  h  + 

V 

by +2 

)  +  +/ 

4 

tj  -4+2 

-tan 


4 

4  + 

4+1 

)+A 

4 

tj  4+ 1 

>•  tj  '4+2  >tj  '4+1  >°- 


(5.16) 


If 


A.  <  a/3  10 


4 

4  + 

4+2  ) 

tj 

b 

/'+ 2 

(5.17) 


with  M  the  number  of  significant  digits  desired,  then  we  replace  (5.16)  by 


K0J  =  In 


by+2 

+ 

tj  *  4+2 

p 

_  |/7L 

44+2 

1 

o* 

' — 

+ 

4+1 

+ 

tj  '4+1 

7 

rl 

h(  h  +  bj+2  )  +  A72 

h  (  h  +  by +,  ) 

i+42 

tj  'b/+2  >  tj  •  b  /+1  >0 


(5.18) 


while,  if  h  -  0, 


K0J  =  In 


p 

by+2 

+ 

tj  'by  +  2 

> 

V 

4+1 

+ 

h  ’  4+i 

7 

,  h  0  ,  tj  -b;.+2  >  tj  ■  b  /+t  >  0  . 


(5.19) 


Similarly, 


K0J=-fo 


by+2 

+ 

fy  '  4+2 

p 

1  N 

4+i 

+ 

tj  '  4+1 

7 

4 

tan 


4 

4  + 

by+2 

)+A 

4 

tj  '4+2 

-tan 


4 

4  + 

4+1 

)+A 

4 

tj  ’  4+i 

>-  tj  '  b/+!  <tj  ■  b/+2  <0 


(5.20) 


and,  if 
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*1 

( h  + 

*V  ) 

b 

;+i 

then 


K0J  =-ln 


4+2 

+ 

h '  b/+2 

4+1 

+ 

'4+i 

+  h 


tj  '  4+2 

<J  '4+1 

4(4+ 4+2  ^ 

-g 

<N 

'gf 

+ 

(4+4+1  , 

+ 

v> 

while,  if  h  =  0, 


Koj  =-ln 


4+2 

+ 

'4+2 

4+i 

+ 

tj  '4+ 1 

,  h  0  ,  tj  -b^+1  <  tj  -b  j+2  <  0  . 


Finally, 


K0J=fo 


(|b;+l|  +  |^/'  ‘  b  y+l  I )  ( |b  y+2  I  +  1^/  'b/+2 
B  2 


-  — tan- 
A 


tj '  b7-+1  4 


A;  +  h  b  .+1  + 


-—tan 

A. 


tj  ‘ b /+2  4 


Ar  +|A|(  b;+2  +|/?| 


t.-b.+i<°,  tj  -b/+2  > 0 


and,  if 


tj-  V.4 


<io-8 


we  replace  (5.24)  by 


K0J  =  In 


tan- 

4 


K, 

+ 

^  4+ 1 ) 

+ 

tj  '  4+2 

f 

h  -4+i 

A 

1 

42 

42  + h  \ 

+  /z) 

t.-b.+24 


42+4(b2+2  +4, 


,  tj  -b/+1  <  0 ,  ti  ■  b  (+2  >  0 


j  "j+ 2 


(5.21) 


(5.22) 


(5.23) 


(5.24) 


(5.25) 


(5.26) 
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while,  if 


<10 


then  we  replace  (5.24)  by 


K0J=  In 


( b;+i  +  tj  •  by+j )( by.+2  +  tj  •  by.+2 ) 

\h\ 

_  1  1  ton  1 

1 

+ 

_ 1 

1 

<N 

_ i 

g 

H 

42  + h  ( b7+i  + h ) 

h  *  b./+2 

h 

42  + h  1 

!b7+2 

+  h) 

0  ’  b7+l  -  0  ’  tj'bJ+ 2^0- 

If  both  (5.25)  and  (5.27)  obtain,  then  we  replace  (5.24)  by 


K0J  =ln 


(■V. 

+ 

h  *  b  7  +1 

)( 

b7'+2 

+ 

h  'b7+2 

f 

h  "b7+l 

|  h 

j 

B 

J 

l 

Aj2  +  h  | 

+  /*) 

h  *  b7+2 

h 

Ar  +  h  | 

!b7-+2 

+  h) 

Vb/  +  I  -  °,  Vb/+2^0- 


7  "j+ 2 


If  h  -  0 ,  then  A  .  *  0  and 


K-oj  =  hi 


(Ki|+Ivv.|)(K>2I+Ivv2|) 


B 


,  h  =  0  ,  tj-  b/+l  <  0  ,  f .  b  ;+2  >  0 


7  7+2 


(5.27) 


(5.28) 


(5.29) 


(5.30) 


The  rest  of  the  Kn  can  be  computed  from  the  iteration  formulas  (5.5)  and  (5.6).  In  Appendix  B 

of  reference  1,  we  show  that  the  iteration  fonnula  does  not  converge  for  all  positions  of  the 
observation  point.  We  discuss  this  issue  in  the  concluding  section. 
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SECTION  6:  VALIDATION 


In  this  section,  we  provide  validation  for  our  numerical  approach  to  evaluating  the  surface  and 
line  integrals  in  (1.16).  We  begin  with  the  surface  integral,  given  by  (2.1) 


F(r') 


1  r  <3 

4ft  j  dz 


(6.1) 


We  use  the  Khayat  and  Wilton  test  triangle  (reference  13),  shown  in  Figure  6.1,  and  a  10-m 
wavelength.  We  choose  the  35  OPs  displayed  in  Tables  6.1  -  6.11.  All  distances  are  in  meters 
and  all  calculations  are  perfonned  in  double  precision  (DP).  The  number  of  terms  in  the 
expansion  of  sine  and  cosine  is  8  each,  providing  15  SDs  of  accuracy.  The  reason  for  having 
eight  terms  rather  than  seven  (as  suggested  in  Table  3.1)  is  that  the  hypotenuse  of  the  triangle  is 
longer  than  a  tenth  of  a  wavelength.  Of  the  35  OPs,  3  are  off  the  triangle's  plane.  Of  the  rest,  16 
are  inside  the  triangle  and  16  outside.  From  Sections  2  and  4,  the  result  for  the  surface  integral 
for  a  point  on  the  plane  of  the  triangle  and  inside  the  triangle  should  be  (-0.5  +  /0);  for  a  point 
outside  the  triangle,  it  should  be  (0  +  /0).  For  all  the  exterior  points,  we  obtained  the  correct 
result  to  15  SDs,  as  shown  in  Table  6.12. 

For  the  interior  points,  shown  in  Table  6.13,  the  results  are  mixed.  For  the  real  part,  we  get 
agreement  to  15  SDs  for  9  out  of  the  16  outcomes;  for  the  rest,  we  get  agreement  to  14  SDs.  We 
attribute  this  to  round-off  error  rather  than  failure  of  the  algorithm.  For  the  imaginary  part,  the 
agreement  is  not  as  good.  Although  we  get  numbers  of  the  order  of  10  to  the  negative  17  or  18, 
these  numbers  are  not  close  enough  to  zero,  the  correct  answer.  This  is  because  the  final  stage  of 
the  calculation  involves  the  subtraction  of  two  almost  identical  numbers.  For  this  case,  we 
developed  a  special  formula  that  is  valid  when  z'  is  appropriately  close  to  the  triangle's  surface.  It 
is  given  by  (2.10)  and  we  repeat  it  here 


N-2 


Im{F(r')}  =  z'X 


(-1)>  +  I)(2^r 

(2n  +  3)\ 


plndS 


(6.2) 


We  present  the  derivation  of  this  formula  in  Appendix  B.  The  integral  in  (6.2)  can  be  evaluated 
either  analytically  or  by  the  GKQ  (as  implemented  in  Matlab  (reference  7))  when  converted  to  a 
line  integral  (see  Appendix  B)  or  by  cubatures.  We  provide  a  typical  example  of  the  use  of  (6.2) 
in  Figure  6.2.  We  note  that  (6.2)  takes  over  when  the  distance  from  the  triangle  is  less  than  about 
10  to  the  -8.  That  two  triangles  will  be  so  near  each  other  is  unlikely  to  happen  in  practice.  It  is 
prudent,  however,  to  know  the  limitations  of  an  algorithm  and  to  provide  alternatives. 
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3:  (-1/3,  2/3,  0) 


F  igure  6.1:  TestTriangle 


Table  6. 1 :  Points  on  the  z  =  0  plane  and  inside  the  triangle.  They  lie  on  the  line  y  =  x . 


Title 

Centroid  Coordinates 

Case  1 

(-7/30,  -7/30,  0)  | 

Case  2 

(-2/15,  -2/15,  0) 

Case  3 

(-1/30,  -1/30,  0) 

Case  4 

(1/15,1/15,0)  ; 

Table  6.2:  Points  on  the  z  =  0  plane  and  outside  the  triangle.  They  lie  on  the  line  y  =  x . 


Title 

Centroid  Coordinates 

Case  5 

(4/15,  4/15,  0)  ! 

Case  6 

(11/30,11/30,0)  | 

Case  7 

(7/15,7/15,0)  | 

Case  8 

(17/30,17/30,0)  ! 

Table  6.3:  Points  above  the  triangle. 


Title 

Centroid  Coordinates 

Case  9 

(-7/30,-7/30,0.1) 

Case  10 

(-7/30,-7/30,0.01)  i 

Case  1 1 

(-7/30,-7/30,0.001)  ; 
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Table  6.4:  Points  on  the  z-0  plane,  inside  the  triangle  and  along  the  median  of  the  upper  left 

angle:  (- 1/3+1 0'*/2,  2/3-10'*,  0),  N  =  1,  2,  3. 


Title 

Centroid  Coordinates 

Case  12 

(-1/3+0. 1/2,  2/3-0. 1,0) 

Case  13 

(-1/3+0.01/2,2/3-0.01,0) 

Case  14 

(-1/3+0.001/2,2/3-0.001,0) 

Table  6.5:  Points  on  the  z  =  0  plane,  outside  the  triangle  and  along  the  median  of  the  upper  left 

angle:  (- 1/3-1 0'*/2,  2/3+10'*,  0),  N =  1,  2,  3. 


j  Title 

Centroid  Coordinates 

Case  15 

(-1/3-0. 1/2,  2/3+0. 1,0) 

!  Case  16 

(-1/3-0.01/2,2/3+0.01,0) 

i  Case  17 

(-1/3-0.001/2,2/3+0.001,0) 

Table  6.6:  Points  on  the  z  =  0  plane,  inside  the  triangle  and  along  the  bisector  of  the  upper  left 

(  i  2  1 0~*  ^ 

--  +  10-*,  — — -j= — ,  0  ,  /V  =  1,2,3 
^  3  3  V2-1 


angle: 


;  Title 

Centroid  Coordinates 

Case  1 8 

f  1  m  2  01  n) 
l  3  3  V2-1  J 

Case  19 

f  1  nni  2  0.01 
l  3  3  V2-I  J 

Case  20 

/ 

V 

1  nnni  2  0.001  A 

--  +  0.001, - -= — ,  0 

3  3  V2-1 

\ 

) 
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Table  6.7:  Points  on  the  z  =  0  plane,  outside  the  triangle  and  along  the  bisector  of  the  upper  left 

r  1  2  10_JV  ^ 

angle:  ---10",-  +  ^— ,0  ,  JV  =  1,2,3 
3  3  V2-1 


V 


J 


V 

Centroid  Coordinates 

Case  2 1 

f  1  n  ,  2  0.1  n) 
L  3  3  V2-1  J 

Case  22 

f  1  2  0.01 

l  3  3  V2-1  J 

i  Case  23 

/ 

V 

3  3  V2-1 

\ 

) 

Table  6.8:  Points  on  the  z  =  0  plane,  inside  the  triangle  and  along  the  median  of  the  lower  right 

angle:  (2/3-  10'w,- 1/3+ 10^/2,  0),  A=l,2,  3. 


Title 

Centroid  Coordinates 

Case  24 

(2/3-0. 1,  - 1/3+0. 1/2,  0) 

Case  25 

(2/3-0.01,-1/3+0.01/2,  0) 

Case  26 

(2/3-0.001,-1/3+0.001/2,  0) 

Table  6.9:  Points  on  the  z  =  0  plane,  outside  the  triangle  and  along  the  median  of  the  lower  right 

angle:  (2/3+ 10'*,- 1/3- 10^/2,  0),  A  =  1,  2,  3. 


i  Title 

Centroid  Coordinates 

,  Case  27 

(2/3+0. 1,  -1/3-0. 1/2,  0) 

1  Case  28 

(2/3+0.01,-1/3-0.01/2,  0) 

Case  29 

(2/3+0.001,-1/3-0.001/2,  0) 
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Table  6.10:  Points  on  the  z  =  0  plane,  inside  the  triangle  and  along  the  bisector  of  the  lower  right 


if210 
an§le:  A — A 


--^^,--+io-w,o 

3  V2-1  3 


N  =  1,2,3. 


Table  6.11:  Points  on  the  z  =  0  plane,  outside  the  triangle  and  along  the  bisector  of  the  lower 

( ?  i(rw  1  ^ 

right  angle:  -  +  -j= — ,—  —  — 10  v  ,  0  ,  A  =  1,2,3. 

13  V2-1  3 
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Table  6.12:  Results  for  the  surface  integral  (6.1)  for  points  on  the  triangle's  plane  but  outside  the 

triangle.  Algorithms  yield  exact  answer. 


RE{INTG_2> 

IM{INTG_2> 

!  Case  5 

0.000000000000000E+00 

0.000000000000000E+00 

Case  6 

0.000000000000000E+00 

0.000000000000000E+00 

Case  7 

0.000000000000000E+00 

0.000000000000000E+00 

Case  8 

0.000000000000000E+00 

0.000000000000000E+00 

|  Case  15 

0.000000000000000E+00 

0.000000000000000E+00 

Case  16 

0.000000000000000E+00 

0.000000000000000E+00 

Case  17 

0.000000000000000E+00 

0.000000000000000E+00 

Case  21 

0.000000000000000E+00 

0.000000000000000E+00 

Case  22 

0.000000000000000E+00 

0.000000000000000E+00 

Case  23 

0.000000000000000E+00 

0.000000000000000E+00 

i  Case  27 

0.000000000000000E+00 

0.000000000000000E+00 

J  Case  28 

0.000000000000000E+00 

0.000000000000000E+00 

j  Case  29 

0.000000000000000E+00 

0.000000000000000E+00 

j  Case  33 

0.000000000000000E+00 

0.000000000000000E+00 

Case  34 

0.000000000000000E+00 

0.000000000000000E+00 

Case  35 

0.000000000000000E+00 

0.000000000000000E+00 

Table  6.13:  Results  for  the  surface  integral  (6.1)  for  points  on  the  triangle's  plane  and  inside  the 

triangle. 


RE{INTG_2> 

IM{INTG_2> 

Case  1 

-5. 00000000000000 E-01 

2. 7 1 5926442076 18E- 17  ! 

Case  2 

-4. 99999999999999 E-01 

-1. 124 18212759  590  E- 17 

Case  3 

-5.00000000000000E-01 

-1. 075731 84301 296 E- 18 

Case  4 

-5.00000000000000E-01 

-5. 6209 1063797953 E- 18 

Case  12 

-4. 99999999999999 E-01 

-2.31341638554094E-17 

Case  13 

-4. 99999999999999 E-01 

-1.01101852584273E-17 

Case  14 

-5.00000000000000E-01 

- 1. 09368894 149475 E- 17 

Case  18 

-5.00000000000000E-01 

-2. 60750622482763 E- 17 

Case  19 

-5.00000000000000E-01 

-6.23551774450725E-17 

Case  20 

-4. 99999999999999 E-01 

-7. 382061 54 191 067 E- 17 

Case  24 

-4. 99999999999999 E-01 

-2.31341638554094E-17  ! 

Case  25 

-4. 99999999999999 E-01 

-1.01101852584273E-17  1 

Case  26 

-5.00000000000000E-01 

- 1. 09368894 149475 E- 17  ! 

Case  30 

-5.00000000000000E-01 

-2. 60750622482763 E- 17 

Case  31 

-5.00000000000000E-01 

-6.23551774450725E-17  i 

Case  32 

-4. 99999999999999 E-01 

-7. 382061 54 191 067 E- 17 
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Comparison  of  imaginary  part  of  surface  integral 
for  old  (in  black)  and  new  (in  blue)  algorithm 


Distance  from  triangle  (A) 

l.E-26  l.E-23  l.E-20  l.E-17  l.E-14  l.E-11  l.E-08  l.E-05  l.E-02 


l.E-01 


l.E-04 
l.E-07  1 

03. 

5* 

QJ 

l.E-10  ^ 

T3 

0) 

l.E-13  * 

-h 

V) 

l.E-16  | 

QJ 

r> 

ID 

l.E-19  s' 

r-* 

n> 

era 

l.E-22  3 


l.E-25 


l.E-28 


Figure  6.2:  Comparison  of  imaginary  part  of  surface  integral  from  old  (6.1)  and  new  (6.2) 
algorithm.  Distance  is  measured  in  wavelengths. 

We  proceed  next  to  a  closer  examination  of  the  surface  integral.  Since  its  calculation  depends  on 
an  iteration  formula  that  has  as  a  starting  point  at  the  solid-angle  integral,  we  test  the  sensitivity 
of  the  initial  term.  To  this  end,  we  use  (4.2)  with  each  of  expressions  (4.3)  and  (4.4).  We  refer  to 
them  as  Formula  1  and  Formula  2.  Clearly,  Formula  1  is  more  economical  than  Formula  2.  We 
used  this  formula  to  compute  (6.1)  above;  thus,  it  is  well  behaved  for  OPs  in  the  triangle's  plane, 
both  inside  and  outside  the  triangle.  We  proceed  to  check  this  formula  for  OPs  near  one  of  the 
triangle's  sides.  We  include  here  Formula  2  as  a  check  on  Formula  1.  We  test  these  formulas 
using  the  points  of  Table  6.14.  These  points  lie  on  the  triangle's  plane  and  are  approached 
through  values  of  z'  that  are  also  given  in  Table  6.14.  The  first  point  (Case  36)  is  the  same  as  for 
Case  1 .  The  next  two  are  the  vertices  of  the  base  of  the  triangle.  Case  39  is  the  middle  point  of 
the  hypotenuse  while  Case  40  is  a  point  along  the  bisector  of  the  right  angle,  and  just  outside  the 
triangle.  As  a  reminder,  the  OP  is  not  allowed  to  lie  on  the  sides  or  vertices  of  the  IT;  thus,  for  z' 
equal  to  zero,  cases  36-38  are  here  to  test  the  solid-angle  formulas  since  we  know  the  outcome  at 
these  points  when  z'  is  equal  to  zero.  The  question  of  how  close  the  OP  can  be  to  a  side  is  one 
that  depends  on  wavelength.  We  do  not  anticipate  situations  where  the  distance  will  be  less  than 
a  millionth  of  a  wavelength  (1  p/t)  and  we  set  this  number  as  the  lower  limit.  For  the  10-m 
wavelength  we  are  using,  this  distance  is  10  pm. 
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Table  6.14:  Points  in  triangle's  plane  for  testing  Formulas  1  and  2. 
The  values  of z' are  5,  1,  10'2,  10'4,  10'6,  and  0. 


i  Title 

Centroid  Coordinates 

Case  36 

(-7/30,  -7/30) 

Case  37 

(-1/3, -1/3) 

Case  38 

(2/3, -1/3) 

Case  39 

(1/6,  1/6) 

Case  40 

(1/5,  1/5) 

We  perform  the  calculations  for  the  points  in  Table  6.14  in  single  and  double  precision  and 
exhibit  the  results  in  Table  6.15.  We  discuss  each  case  separately.  In  Case  36,  the  OP  approaches 
a  point  in  the  interior  of  the  triangle.  All  DP  results  but  the  first  agree  to  15  SDs.  The  first  pair 
agrees  to  14  SDs.  The  single  precision  (SP)  results  for  Fonnula  1  are  better  than  those  of 
Formula  2.  All  four  results  are  exact  when  z’  is  zero. 

In  Case  37,  the  OP  approaches  the  90-deg  corner  of  the  triangle.  The  correct  result  at  the  corner 
point  is  0.25.  None  of  formulas  gives  this  result.  In  fact,  Fonnula  1  sees  the  comer  point  as  being 
outside  the  triangle  while  Fonnula  2  sees  it  inside  the  triangle.  This  may  explain  the  diverging 
results  in  DP  as  that  point  is  approached.  We  also  see  that  Formula  1  is  much  better  behaved  in 
SP  than  Formula  2. 

In  Case  38,  the  OP  approaches  the  lower-right  corner  of  the  triangle.  The  correct  answer  at  the 
comer  point  is  0.125.  Both  formulas  see  this  point  as  being  outside  the  triangle  and  give  the 
wrong  result  (zero).  For  the  rest  of  the  OPs,  there  is  perfect  agreement  between  the  two  formulas 
in  DP  except  when  z’  is  equal  to  5  m,  in  which  case  the  agreement  is  to  14  SDs.  The  SP  results 
are  good  except  for  one  occasion  for  Formula  2. 

In  Case  39,  the  OP  approaches  the  midpoint  of  the  hypotenuse.  At  that  point,  the  correct  answer 
is  0.5.  Clearly,  in  DP  both  formulas  see  this  point  outside  the  triangle.  In  SP,  Formula  1  sees  it 
inside  while  Formula  2  sees  it  outside.  The  SP  results  are  not  very  good  in  general,  while  the  DP 
results  provide  from  10  to  15  SDs. 

In  Case  40,  the  OP  approaches  a  point  in  the  triangle's  plane  but  slightly  outside  the  triangle. 
Both  formulas  give  the  exact  result  when  z’  becomes  zero.  In  DP,  there  is  agreement  to  at  least 
14  SDs.  In  SP,  they  mostly  provide  correct  answers  to  six  SDs. 

The  conclusion  we  can  draw  from  these  runs  is  that  DP  should  be  used  when  possible  and  that 
Formula  1  performs  respectably  provide  the  OP  is  not  very  close  over  a  side  of  the  triangle.  In 
Appendix  D,  we  present  a  more  detailed  comparison  when  z'  is  zero  or  very  close  to  it. 
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Table  6.15:  Results  for  Formulas  1  and  2  for  the  OPs  in  Table  6.14  in  SP  and  DP. 

Numbers  in  DP  occupy  two  rows.  Numbers  in  grey  background  differ  from 
one  another.  Numbers  in  SP  in  red  do  not  provide  7  SD  according  to  the  DP  results. 
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z' 

5.00E+00 

1.00E+00 

1.00E-02 

1.00E-04 

1.00E-06 

0.00E+00 

Formula  1 

3.141805E-03 

6.106466E-02 

9.416611E-01 

9.994146E-01 

9.999941E-01 

1.000000E+00 

Formula  2 

3.141849E-03 

6.106465E-02 

9.416611E-01 

9.994146E-01 

9.999941E-01 

1.000000E+00 

3.1418053817 

6.106466233674 

9.416610540466 

9.994146187948 

9.999941461859 

1.000000000000 

DP 

Formula  1 

4860E-03 

18E-02 

21E-01 

04E-01 

44E-01 

00E+00 

3.1418053817 

6.106466233674 

9.416610540466 

9.994146187948 

9.999941461859 

1.000000000000 

Formula  2 

4867E-03 

18E-02 

21E-01 

04E-01 

44E-01 

00E+00 

Case 

37 

z' 

5.00E+00 

1.00E+00 

1.00E-02 

1.00E-04 

1.00E-06 

0.00E+00 

SP 

Formula  1 

3.120885E-03 

5.408673E-02 

2.468172E-01 

2.499682E-01 

2.499997E-01 

0.000000E+00 

Formula  2 

3.120894E-03 

5.408673E-02 

2.468187E-01 

2.499841E-01 

2.499998E-01 

1.000000E+00 

3.1208851602 

5.408672398469 

2.468171663621 

2.499681690116 

2.499996816901 

0.000000000000 

DP 

Formula  1 

7894E-03 

63E-02 

87E-01 

46E-01 

13E-01 

00E+00 

3.1208851602 

5.408672398469 

2.468171663621 

2.499681690116 

2.499996816759 

1.000000000000 

Formula  2 

7894E-03 

63E-02 

92E-01 

24E-01 

64E-01 

00E+00 

Case 


Case 
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z' 

5.00E+00 

1.00E+00 

1.00E-02 

1.00E-04 

1.00E-06 

0.00E+00 

SP 

Formula  1 

3.061427E-03 

4.166667E-02 

1.238747E-01 

1.249887E-01 

1.249999E-01 

0.000000E+00 

Formula  2 

3.061440E-03 

4.166667E-02 

1.238747E-01 

1.249887E-01 

1.249999E-01 

0.000000E+00 

3.0614270888 

4.166666666666 

1.238746514932 

1.249887460460 

1.249998874604 

0.000000000000 

DP 

Formula  1 

5004E-03 

66E-02 

54E-01 

94E-01 

60E-01 

00E+00 

3.0614270888 

4.166666666666 

1.238746514932 

1.249887460460 

1.249998874604 

0.000000000000 

Formula  2 

5008E-03 

66E-02 

54E-01 

94E-01 

60E-01 

00E+00 

Case 

39 

z' 

5.00E+00 

1.00E+00 

1.00E-02 

1.00E-04 

1.00E-06 

0.00E+00 

SP 

Formula  1 

3.151634E-03 

6.409422E-02 

4.909984E-01 

5.000723E-01 

5.094833E-01 

1.000000E+00 

Formula  2 

3.151670E-03 

6.409420E-02 

4.909987E-01 

4.999069E-01 

4.999820E-01 

0.000000E+00 

3.1516345252 

6.409421684897 

4.909983369786 

4.999099683697 

4.999990997276 

0.000000000000 

DP 

Formula  1 

0085E-03 

49E-02 

01E-01 

99E-01 

50E-01 

00E+00 

3.1516345252 

6.409421684897 

4.909983369786 

4.999099683698 

4.999990996872 

0.000000000000 

Formula  2 

0092E-03 

49E-02 

00E-01 

99E-01 

48E-01 

00E+00 

Case 

40 

z' 

5.00E+00 

1.00E+00 

1.00E-02 

1.00E-04 

1.00E-06 

0.00E+00 

SP 

Formula  1 

3.147065E-03 

6.223706E-02 

5.823865E-02 

5.922366E-04 

5.922377E-06 

0.000000E+00 

Formula  2 

3.147028E-03 

6.223706E-02 

5.823865E-02 

5.922369E-04 

5.922379E-06 

0.000000E+00 

3.1470641589 

6.223706005980 

5.823863498580 

5.922367555002 

5.922377670355 

0.000000000000 

DP 

Formula  1 

7432E-03 

60E-02 

16E-02 

76E-04 

95E-06 

00E+00 

3.1470641589 

6.223706005980 

5.823863498580 

5.922367555002 

5.922377670355 

0.000000000000 

Formula  2 

7428E-03 

60E-02 

16E-02 

75E-04 

95E-06 

00E+00 

Now  that  we  chose  Formula  1  (Equations  (4.2)  and  (4.3))  as  the  starting  term  in  the  iteration 
scheme  of  Section  4,  we  proceed  to  compute  (6.1),  the  surface  integral  of  interest,  for  the  points 
of  Table  6.14  except  for  z'  equal  to  zero.  The  reason  for  excluding  this  point  is  that  we  have 
tested  interior  and  exterior  points  in  the  triangle's  plane  above  and  that  points  on  the  triangle's 
boundary  are  not  allowed.  In  its  place,  we  display  the  value  4.5.  In  Table  6.16,  we  display  the 
convergence  of  our  algorithm  for  Case  36.  The  tables  for  the  rest  of  the  cases  are  presented  in 
Appendix  E  so  as  not  to  clutter  the  main  body  of  the  report. 


37 


NAWCADPAX/TR-20 1 2/200 


We  note  that  on  1  occasion  in  this  table  (and  1  each  in  Tables  E.2  and  E.3,  a  total  of  3  occasions 
out  of  60)  our  formula  fails  to  predict  the  number  of  terms  in  the  expansions  of  Section  4 
correctly.  In  these  3  instances,  and  for  the  real  part  of  the  integral,  9  terms  are  required;  still,  the 
predicted  8  terms  yield  14  SD.  We  hesitate  to  change  our  formula  since,  as  can  be  seen  from  the 
tables,  it  is  a  rather  conservative  one. 

Table  6.16:  Convergence  of  our  algorithm  for  the  surface  integral  for  the  OP  of  Case  36, 
Table  6. 14.  Convergence  is  to  15  SD  and  as  a  function  of  the  number  of  terms  in  the 
trigonometric  sums  in  (2.2).  The  entries  in  light  green  are  the  first  one  for  which  our 

formula  has  stabilized  to  15  SD. 


z'  = 

5.000 

4.500 

REAL{F(r')} 

IMAG{F(r’)} 

REAL{F(r')} 

IMAG{F(r')} 

15  SD,  Nc=Ns=4 

1. 63858509704585  E-03 

4. 9343951 1546236E-03 

2.28919850402077E-04 

5.82159533394148E-03 

15  SD,  Nc=Ns=5 

1. 63858509704592  E-03 

4. 9343951 1546241E-03 

2. 289 1985040209  IE-04 

5. 82 15953339415  IE-03 

15  SD,  Nc=Ns=6 

1.63858509704589E-03 

4. 9343951 1546239E-03 

2.28919850402082E-04 

5.82159533394150E-03 

15  SD,  Nc=Ns=7 

1.63858509704590E-03 

4. 9343951 1546240E-03 

2.28919850402084E-04 

5.82159533394150E-03 

15  SD,  Nc=Ns=8 

1.63858509704589E-03 

4. 9343951 1546239E-03 

2.28919850402084E-04 

5.82159533394150E-03 

15  SD,  Nc=Ns=9 

1.63858509704590E-03 

4. 9343951 1546239E-03 

2.28919850402084E-04 

5.82159533394150E-03 

15  SD,  Nc=Ns=10 

1.63858509704590E-03 

4. 9343951 1546239E-03 

2.28919850402084E-04 

5.82159533394150E-03 

15  SD,  Nc=Ns=ll 

1.63858509704590E-03 

4. 9343951 1546239E-03 

2.28919850402084E-04 

5.82159533394150E-03 

15  SD,  Nc=Ns=12 

1.63858509704590E-03 

4. 9343951 1546239E-03 

2.28919850402084E-04 

5.82159533394150E-03 

z'  = 

1.000 

0.010 

REAL{F(r')} 

IMAG{F(r’)} 

REAL{F(r')} 

IMAG{F(r’)} 

15  SD,  Nc=Ns=4 

-3. 686050508 11883E-02 

3.13416887059356E-03 

-4.71122428585698E-01 

3.26142239362246E-05 

15  SD,  Nc=Ns=5 

-3.68605050809884E-02 

3. 13416887043 198E-03 

-4.71 122428617225  E-01 

3.26142297742437E-05 

15  SD,  Nc=Ns=6 

-3.68605050809884E-02 

3.13416887043203E-03 

-4.71 1224286172 12  E-01 

3. 261422977 17 19  IE-05 

15  SD,  Nc=Ns=7 

-3.68605050809884E-02 

3.13416887043203E-03 

-4.7 11224286172 12  E-01 

3. 26 1422977 1719  IE-05 

15  SD,  Nc=Ns=8 

-3.68605050809884E-02 

3.13416887043203E-03 

-4.7 11224286172 12  E-01 

3. 26 1422977 1719  IE-05 

15  SD,  Nc=Ns=9 

-3.68605050809884E-02 

3.13416887043203E-03 

-4.7 11224286172 12  E-01 

3. 26 1422977 1719  IE-05 

15  SD,  Nc=Ns=10 

-3.68605050809884E-02 

3.13416887043203E-03 

-4.7 11224286172 12  E-01 

3. 26 1422977 1719  IE-05 

15  SD,  Nc=Ns=ll 

-3.68605050809884E-02 

3.13416887043203E-03 

-4.7 11224286172 12  E-01 

3. 26 1422977 1719  IE-05 

15  SD,  Nc=Ns=12 

-3.68605050809884E-02 

3.13416887043203E-03 

-4.7 11224286172 12  E-01 

3.26142297717191E-05 

z'  = 

0.0001 

0.000001 

REAL{F(r')} 

IMAG{F(r’)} 

REAL{F(r')} 

IMAG{F(r’)} 

15  SD,  Nc=Ns=4 

-4.997103232 12000E-01 

3.26135850866057E-07 

-4. 99997 103199681  E-01 

3. 2536773246842  IE-09 

15  SD,  Nc=Ns=5 

-4.99710323253529E-01 

3. 26 143592 14025  IE-07 

-4.99997 103241321  E-01 

3.26143966793828E-09 

15  SD,  Nc=Ns=6 

-4.99710323253509E-01 

3.26143588367667E-07 

-4. 99997 103241301  E-01 

3.26143589320174E-09 

15  SD,  Nc=Ns=7 

-4.99710323253509E-01 

3.26143588367667E-07 

-4. 99997 103241301  E-01 

3.26143589320174E-09 

15  SD,  Nc=Ns=8 

-4.99710323253509E-01 

3.26143588367667E-07 

-4. 99997 103241301  E-01 

3.26143589320174E-09 

15  SD,  Nc=Ns=9 

-4.99710323253509E-01 

3.26143588367667E-07 

-4. 99997 103241301  E-01 

3.26143589320174E-09 

15  SD,  Nc=Ns=10 

-4.99710323253509E-01 

3.26143588367667E-07 

-4. 99997 103241301  E-01 

3.26143589320174E-09 

15  SD,  Nc=Ns=ll 

-4.99710323253509E-01 

3.26143588367667E-07 

-4. 99997 103241301  E-01 

3.26143589320174E-09 

15  SD,  Nc=Ns=12 

-4.99710323253509E-01 

3.26143588367667E-07 

-4. 99997 103241301  E-01 

3.26143589320174E-09 

38 


NAWCADPAX/TR-20 1 2/200 


We  compare  next  the  values  of  the  surface  integral  for  the  five  points  of  Table  6.14.  In  Tables 
6. 17. a  -  6.17.C,  we  display  the  results  for  Case  36.  In  the  first  row  of  each  table  we  display  the 
values  of  z'  in  meters.  In  the  second  row,  we  present  the  required  number  of  terms  in  (2.2)  of 
Section  2  to  obtain  15  SD  (8  terms  each  for  the  sine  and  the  cosine).  We  follow  that  by  a 
Riemann  sum  of  the  integral  that  involves  10,000  points,  not  a  very  economical  method 
especially  since  it  does  not  provide  the  required  15  SD.  We  present  the  cubature  results  after  the 
Riemann  sum,  first  for  the  real  part  of  the  integral  and,  next,  for  the  imaginary  part.  In  Appendix 
F,  we  present  the  methodology  we  used  to  compute  the  integral  using  cubatures.  Although  this 
procedure  is  not  necessary  when  the  OP  is  as  far  as  half  a  wavelength  from  the  IT,  it  is  absolutely 
necessary  when  the  OP  is  near  the  IT.  The  methodology  is  not  new  (see,  for  example,  reference 
14).  The  way,  however,  we  computed  the  integral  of  the  static  Green's  function  is  ours.  This  is  a 
time-consuming  way  to  compute  using  cubatures  but  it  is  used  here  for  the  single  purpose  of 
comparing  with  our  method  since  we  do  not  have  any  other  data  available.  We  see  that,  with  a 
28-point  cubature,  we  get  agreement  to  14  SD  for  the  real  part  of  the  integral.  For  the  imaginary 
part,  we  get  15  SD  with  a  54-point  cubature  for  the  5.0  case  and  as  many  SD  for  the  4.5  case  but 
with  a  28-point  cubature.  We  note  that  we  get  14  SD  for  either  case  with  a  21-point  cubature. 
The  terms  highlighted  in  green  are  the  lowest-point  cubature  terms  for  which  we  get  the 
maximum  number  of  SD  to  agree  with  our  result. 

We  remark  that,  for  both  cases  in  Table  6. 17. a,  we  could  have  gotten  the  same  number  of  SD  by 
applying  the  cubatures  directly  to  the  integrals.  This  is  because  the  OP  is  sufficiently  far  away 
from  the  IT.  This  is  no  longer  true  for  the  next  four  cases  shown  in  Tables  6.18  and  6.19.  Here, 
we  have  to  use  the  methodology  of  Appendix  F  to  increase  the  number  of  SD. 

The  three  tables  are  summarized  in  Figure  6.3  where  we  display  the  maximum  number  of  SD 
obtained  using  cubatures  as  a  function  of  distance  from  the  triangle.  The  remaining  four  cases  in 
Table  6.14  are  similar  to  Case  36.  For  completeness,  the  data  are  provided  in  Appendix  G  and 
the  relevant  graphs  in  Figures  6.4  to  6.7.  From  Table  6.14  (and  the  discussion  following  it),  Case 
36  involves  an  OP  whose  projection  is  an  interior  point  of  the  triangle.  In  the  next  two  cases,  the 
OPs  project  to  the  lower  comers  of  the  triangle,  while,  in  Case  39,  to  the  midpoint  of  the 
hypotenuse.  In  Case  40,  the  OP  projects  to  a  point  just  outside  the  triangle.  Thus,  in  the  first  four 
cases,  the  OP  approaches  a  singular  point  while  in  the  fifth  it  approaches  a  regular  point  that  is 
near  a  singular  point. 

For  the  first  four  cases,  we  see  that  the  imaginary  part  (which  is  computed  using  a  cubature  only) 
has  a  quasi-linear  behavior  and  that  it  loses  SD  as  the  OP  approaches  the  IT.  At  the  same  time, 
the  order  of  the  cubature  that  provides  the  greatest  number  of  SD  increases.  In  the  case  of  the 
real  part  (involving  two  integrals  that  are  evaluated  analytically  and  one  through  cubature),  the 
result  is  a  U-shaped  graph;  thus,  in  the  approach  of  the  OP  to  the  IT,  we  have  a  good  number  of 
SD  when  far  away,  we  lose  a  few  along  the  way,  and  we  recover  them  as  we  get  close  to  the  IT. 
For  the  last  case  (where  the  OP  projects  to  a  point  outside  the  triangle),  the  behavior  of  the  real 
part  is  essentially  the  same  as  before.  The  imaginary  part,  however,  quickly  settles  to  13  SD  and 
does  not  fall  below  that  level. 
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It  is  clear  from  the  above  discussion  that  we  put  more  trust  in  the  results  of  our  method  and  much 
less  in  those  provided  by  cubatures.  The  reasons  for  it  are  the  following.  When  z'  is  large  (greater 
than  a  tenth  of  a  wavelength),  we  get  15-SD  agreement  between  our  results  and  those  provided 
by  cubatures.  When  z'  is  less  than  that,  the  cubature  results  fail  to  settle  down  as  the  cubature 
order  increases.  Our  results  derive  from  analytical  formulas.  The  only  way  that  they  can  be 
correct  on  one  occasion  and  wrong  on  another  is  for  them  to  be  unstable  for  certain  values  of  the 
parameters/arguments.  As  we  discussed  earlier,  we  detected  such  instability  in  the  imaginary  part 
and  corrected  for  it.  We  note,  however,  that  this  instability  occurs  for  values  of  z'  smaller  than 
the  ones  we  have  used.  Until  fresh  numbers  appear  from  another  source,  we  will  consider  our 
expressions  as  providing  correct  answers  to  the  requested  number  of  SD. 
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Table  6. 17. a:  Comparison  of  our  results  for  the  surface  integral  of  Case  36  with  those 
obtained  through  a  Riemann  sum  and  cubatures. 


z'  = 

15  SD,  Nc=Ns=8 

MFIE  METHOD 

5.0 

4.5 

REAL{F(r')} 

IMAG{F(r')} 

1.63858509704589  E-03 

4.9343951 1546239  E-03 

2.28919850402084E-04 

5.82159533394150E-03 

RIEMANN  SUM 

REAL{F(r')} 

1.63858509715993  E-03 

2.28919850599920E-04 

N  =  10,000 

IMAG{F(r')} 

4. 934395 115 16726E-03 

5. 8215953336272  IE-03 

CUBATURE  DIM 

SD 

RE{F(r')} 

21 

1.63858509704539E-03 

2.28919850400982E-04 

28 

14,  14 

1.63858509704591E-03 

2. 28919850402085  E-04 

36 

1.63858509704591E-03 

2.28919850402083E-04 

45 

1.63858509704592  E-03 

2.28919850402088E-04 

46 

1.63858509704591E-03 

2.28919850402085E-04 

54 

1.63858509704591E-03 

2.28919850402080E-04 

55 

1.63858509704591E-03 

2.28919850402085E-04 

66 

1.63858509704591E-03 

2.28919850402083E-04 

78 

1.63858509704591E-03 

2.28919850402088E-04 

85 

1.63858509704591E-03 

2.28919850402088E-04 

91 

1.63858509704591E-03 

2.28919850402090E-04 

105 

1.63858509704591E-03 

2.28919850402083E-04 

120 

1.63858509704591E-03 

2.28919850402087E-04 

126 

1.63858509704591E-03 

2.28919850402087E-04 

CUBATURE  DIM 

SD 

IM{F(r')} 

21 

4. 934395 11546240E-03 

5. 8215953339415  IE-03 

28 

15 

4. 934395 11546240E-03 

5.82159533394150E-03 

36 

4. 934395 11546240E-03 

5.82159533394150E-03 

45 

4. 934395 11546241E-03 

5.82159533394152E-03 

46 

4. 934395 11546240E-03 

5.82159533394150E-03 

54 

15 

4. 934395 11546239E-03 

5.82159533394150E-03 

55 

4. 934395 11546240E-03 

5.82159533394150E-03 

66 

4. 934395 11546239E-03 

5.82159533394150E-03 

78 

4. 934395 11546240E-03 

5.82159533394150E-03 

85 

4. 934395 11546240E-03 

5. 8215953339415  IE-03 

91 

4. 934395 11546240E-03 

5.82159533394150E-03 

105 

4. 934395 11546239E-03 

5.82159533394150E-03 

120 

4. 934395 11546240E-03 

5.82159533394150E-03 

126 

4. 934395 11546240E-03 

5.82159533394150E-03 
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Table  6.17.b:  Comparison  of  our  results  for  the  surface  integral  of  Case  36  with  those 
obtained  through  a  Riemann  sum  and  cubatures. 


z'  = 

15  SD,  Nc=Ns=8 

MFIE  METHOD 

1.0E+00 

1.0E-02 

REAL{F(r')} 

IMAG{F(r')} 

-3.68605050809884E-02 

3. 13416887043203  E-03 

-4.71122428617212E-01 

3.26142297717191E-05 

RIEMANN  SUM 

REAL{F(r')} 

-3.68605050787650E-02 

-4. 7 1 12242861923  IE-01 

N  =  10,000 

IMAG{F(r')} 

3.13416887030680E-03 

3.26142297704777E-05 

CUBATURE  DIM 

SD 

RE{F(r')} 

21 

-3. 6860505 1035008E-02 

-4.71122424239280E-01 

28 

-3.68605050827459E-02 

-4.71122429698403E-01 

36 

-3. 6860505081043  IE-02 

-4.71122430974795E-01 

45 

-3.68605050809920E-02 

-4.71122429253321E-01 

46 

-3.68605050809830E-02 

-4.71122428783164E-01 

54 

-3. 68605050809892  E-02 

-4.71122428079805E-01 

55 

-3. 68605050809885  E-02 

-4.71122427827648E-01 

66 

-3. 68605050809885  E-02 

-4.71122427958889E-01 

78 

15 

-3.68605050809884E-02 

-4.71122429539633E-01 

85 

-3.68605050809884E-02 

-4.71122429631584E-01 

91 

-3.68605050809884E-02 

-4.71122429294459E-01 

105 

10 

-3.68605050809884E-02 

-4.71122428577606E-01 

120 

-3.68605050809884E-02 

-4.71122428452062E-01 

126 

-3.68605050809884E-02 

-4.71122428824172E-01 

CUBATURE  DIM 

SD 

IM{F(r')} 

21 

3.13416887043204E-03 

3.26142297717566E-05 

28 

3.13416887043204E-03 

3.26142297717564E-05 

36 

15 

3. 13416887043203  E-03 

3.26142297717573E-05 

45 

11 

3.13416887043204E-03 

3.26142297717561E-05 

46 

3. 13416887043203  E-03 

3.26142297717565E-05 

54 

3. 13416887043203  E-03 

3.26142297717563E-05 

55 

3.13416887043204E-03 

3.26142297717566E-05 

66 

3. 13416887043203  E-03 

3.26142297717565E-05 

78 

3.13416887043204E-03 

3.26142297717561E-05 

85 

3.13416887043204E-03 

3.26142297717577E-05 

91 

3.13416887043204E-03 

3.26142297717573E-05 

105 

3. 13416887043203  E-03 

3.26142297717570E-05 

120 

3.13416887043204E-03 

3.26142297717566E-05 

126 

3. 13416887043203  E-03 

3.26142297717564E-05 
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Table  6.17.c:  Comparison  of  our  results  for  the  surface  integral  of  Case  36  with  those 
obtained  through  a  Riemann  sum  and  cubatures. 


z'  = 

15  SD,  Nc=Ns=8 

MFIE  METHOD 

1.0E-04 

1.0E-06 

REAL{F(r')} 

IMAG{F(r')} 

-4. 99710323253509  E-01 

3.26143588367667E-07 

-4. 99997 103241301  E-01 

3.26143589320174E-09 

RIEMANN  SUM 

REAL{F(r')} 

-4.99710323252418E-01 

-4. 99997 103238861  E-01 

N  =  10,000 

IMAG{F(r')} 

3. 2614358833523  IE-07 

3. 261435884642 11 E-09 

CUBATURE  DIM 

SD 

RE{F(r')} 

21 

-4. 99710323208003  E-01 

-4.99997103240846E-01 

28 

-4. 997 10323263749  E-01 

-4.99997103241403E-01 

36 

-4.99710323279561E-01 

-4.99997103241561E-01 

45 

-4.997 10323261 100E-01 

-4.99997103241377E-01 

46 

-4. 997 10323254762  E-01 

-4.99997103241313E-01 

54 

-4. 997 10323247283  E-01 

-4.99997103241238E-01 

55 

-4. 997 10323244832  E-01 

-4.99997103241214E-01 

66 

-4.99710323246108E-01 

-4.99997103241227E-01 

78 

-4.99710323264184E-01 

-4.99997103241407E-01 

85 

-4. 997 10323265882  E-01 

-4.99997103241424E-01 

91 

-4.99710323263249E-01 

-4.99997103241398E-01 

105 

12,  12 

-4. 99710323253062  E-01 

-4.99997103241296E-01 

120 

-4. 997 1032325 1477E-01 

-4.99997103241280E-01 

126 

-4.99710323255807E-01 

-4.99997103241324E-01 

CUBATURE  DIM 

SD 

IM{F(r')} 

21 

3. 2614358834815  IE-07 

3.26143588477216E-09 

28 

3. 26143588348152  E-07 

3.26143588477219E-09 

36 

3.26143588348140E-07 

3.26143588477203E-09 

45 

3.26143588348144E-07 

3.26143588477210E-09 

46 

3.26143588348154E-07 

3.26143588477214E-09 

54 

3.26143588348150E-07 

3.26143588477214E-09 

55 

3. 26143588348151  E-07 

3.26143588477214E-09 

66 

3. 26143588348151  E-07 

3.26143588477215E-09 

78 

10 

3. 26143588348163  E-07 

3.26143588477196E-09 

85 

3.26143588348147E-07 

3.26143588477194E-09 

91 

8 

3.26143588348247E-07 

3.26143588477247E-09 

105 

3. 26143588348149  E-07 

3.26143588477214E-09 

120 

3.26143588348154E-07 

3.26143588477216E-09 

126 

3.26143588348148E-07 

3.26143588477210E-09 
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Figure  6.3:  Case  36  for  surface  integral:  maximum  number  of  SD  obtained  using  cubatures 
as  a  function  of  distance  from  IT.  A  number  attached  to  a  calculation  point  (z',  vv)  gives 
the  size  of  the  smallest  cubature  that  provides  w  SD  at  the  point  z'. 
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Figure  6.4:  Case  37  for  surface  integral:  maximum  number  of  SD  obtained  using  cubatures 
as  a  function  of  distance  from  IT.  A  number  attached  to  a  calculation  point  (z',  w)  gives 
the  size  of  the  smallest  cubature  that  provides  w  SD  at  the  point  z'. 
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Figure  6.5:  Case  38  for  surface  integral:  maximum  number  of  SD  obtained  using  cubatures 
as  a  function  of  distance  from  IT.  A  number  attached  to  a  calculation  point  (z',  vv)  gives 
the  size  of  the  smallest  cubature  that  provides  w  SD  at  the  point  z'. 
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Figure  6.6:  Case  39  for  surface  integral:  maximum  number  of  SD  obtained  using  cubatures 
as  a  function  of  distance  from  IT.  A  number  attached  to  a  calculation  point  (zr,  w)  gives 
the  size  of  the  smallest  cubature  that  provides  w  SD  at  the  point  z'. 
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Figure  6.7:  Case  40  for  surface  integral:  maximum  number  of  SD  obtained  using  cubatures 
as  a  function  of  distance  from  IT.  A  number  attached  to  a  calculation  point  (zr,  w)  gives 
the  size  of  the  smallest  cubature  that  provides  w  SD  at  the  point  z'. 

We  turn  now  to  the  evaluation  of  the  line  integral.  The  exact  form  of  this  integral  is  given  by  the 
left-hand  side  of  (3.5)  and  repeated  here 

-/&|r-r'| 

=  Jt - jds,  j  =  1,2,3  (6.3) 

J  r  —  r 

Sj  i  i 

where  Sj  stands  for  one  of  the  sides  of  the  triangle.  We  compare  each  side  separately.  The  number 
of  a  side  is  the  same  as  the  number  of  the  vertex  opposite  that  side.  The  comparison  in  this  case 
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is  with  Matlab's  GKQ  algorithm.  We  set  the  relative  tolerance  in  this  subroutine  to  its  minimum, 
2.22*10A(-14).  The  numbers  we  obtained  using  our  method  (and  requiring  15  SD)  agree  with 
those  of  the  GKQ  algorithm  to  at  least  14  SD.  We  present  part  of  these  results  in  Table  6. 18. 


Table  6.18:  Evaluation  of  line  integral  in  (6.3)  and  comparison  with  Matlab's  GKQ. 


!  RE-fINTG  1} 

IM-fINTG  1}  1  MATLAB  QUADGK  ALGORITHM 

CASE  1  SIDE  1  1.90582350126434E+00 

-8. 60420498 19 1187E-01  1.90582350126434E+00 

-8. 60420498191 187E-01 

SIDE  2  3.69035907422533E+00 

-6. 17945 178406258E-01  3.69035907422534E+00 

-6. 17945178406258E-01 

SIDE  3  3.69035907422533E+00 

-6. 17945 178406258E-01  3.69035907422534E+00 

-6. 17945178406258E-01 

|  RE{INTG 1} 

IM{INTG 1>  |  MATLAB  QUADGK  ALGORITHM  | 

CASE  5  SIDE  1  4.51600791903947E+00 

-8. 77728566109 124E-01  4.51600791903947E+00 

-8.77728566109124E-01 

SIDE  2  1.37627563290081E+00 

-6. 09751 10850262  IE-01  1.37627563290082E+00 

-6.09751108502621E-01 

SIDE  3  1.37627563290081E+00 

-6. 09751 10850262  IE-01  1.37627563290082E+00 

-6.09751108502621E-01 

!  RE-fINTGl} 

IM-fINTG l>  |  MATLAB  QUADGK  ALGORITHM  | 

CASE  11  SIDE  1  1.90582086366456E+00 

-8.60420440840094E-01  1.90582086366456E+00 

-8.60420440840094E-01 

SIDE  2  3.69027366103930E+00 

-6. 17945137475203E-01  3.69027366103930E+00 

-6. 17945137475203E-01 

SIDE  3  3.69027366103930E+00 

-6. 17945137475203E-01  3.69027366103930E+00 

-6. 17945137475203E-01 

I!  RE-fINTG  1} 

IM-fINTG  1>  |  MATLAB  QUADGK  ALGORITHM  | 

CASE  14  SIDE  1  1.06141588860891E+01 

-8.50594455858452E-01  1.06141588860895E+01 

-8. 50594455858453E-01 

SIDE  2  9.63978848276054E+00 

-6. 14740566007970E-01  9.63978848276057E+00 

-6. 14740566007970E-01 

SIDE  3  6. 6584199523553 IE-01 

-5.74796541407303E-01  6.65841995235532E-01 

-  5. 7479654 1407304E-01 

|  RE-fINTG  1} 

IM-fINTG  1}  |  MATLAB  QUADGK  ALGORITHM  | 

CASE  17  SIDE  1  6.97821202303442E+00 

-8.50425853556906E-01  6.97821202303448E+00 

-8.50425853556907E-01 

SIDE  2  6.75413585870779E+00 

-6. 14659499455 112E-01  6.75413585870781E+00 

-6. 14659499455112E-01 

SIDE  3  6.63760055152963E-01 

-5.74600762567182E-01  6.63760055152964E-01 

-5.74600762567182E-01 

i  RE{INTG 1} 

IM{INTG 1>  |  MATLAB  QUADGK  ALGORITHM  | 

CASE  20  SIDE  1  9. 37027630463920E+00 

-8.50701732371761E-01  9.37027630463922E+00 

-8.50701732371761E-01 

SIDE  2  9.11674874498877E+00 

-6. 14797666248879E-01  9. 11674874498878E+00 

-6. 14797666248879E-01 

SIDE  3  6.67242696421228E-01 

-5.74926715022366E-01  6.67242696421228E-01 

-5.74926715022366E-01 

!  RE{INTG 1> 

IM{INTG 1>  |  MATLAB  QUADGK  ALGORITHM  | 

CASE  23  SIDE  1  6. 14264876010087E+00 

-8. 503179696 12828E-01  6. 14264876010089E+00 

-8.50317969612829E-01 

SIDE  2  5.89087476299632E+00 

-6. 14601954320508E-01  5.89087476299633E+00 

-6. 14601954320508E-01 

SIDE  3  6.62366864094923E-01 

-5.74470195405343E-01  6.62366864094923E-01 

-5. 74470 195405343 E-01 

i  RE{INTG 1> 

IM-fINTG l>  |  MATLAB  QUADGK  ALGORITHM  | 

CASE  26  SIDE  1  1.06141588860894E+01 

-8.50594455858452E-01  1.06141588860894E+01 

-8.50594455858453E-01 

SIDE  2  6. 6584199523553 IE-01 

-5. 7479654 1407303E-01  6.65841995235532E-01 

-5. 7479654 1407304E-01 

SIDE  3  9. 63978848276054E+00 

-6. 14740566007970E-01  9.63978848276057E+00 

-6. 14740566007970E-01 

S  RE-fINTG  1} 

IM-fINTG  1>  |  MATLAB  QUADGK  ALGORITHM  | 

CASE  29  SIDE  1  6.97821202303442E+00 

-8.50425853556906E-01  6.97821202303446E+00 

-8.50425853556907E-01 

SIDE  2  6. 63760055 152963E-01 

-5.74600762567182E-01  6.63760055152964E-01 

-5.74600762567182E-01 

SIDE  3  6.75413585870779E+00 

-6. 14659499455 112E-01  6.75413585870781E  +  00 

-6. 1465949945 5 112E-01 

!  RE-fINTG  1} 

IM-fINTG  1>  |  MATLAB  QUADGK  ALGORITHM  | 

CASE  32  SIDE  1  9.37027630463920E+00 

-8.50701732371761E-01  9.37027630463922E+00 

-8.50701732371761E-01 

SIDE  2  6.67242696421228E-01 

-5. 749267 15022366E-01  6.67242696421228E-01 

-5.74926715022366E-01 

SIDE  3  9.11674874498877E+00 

-6. 14797666248879E-01  9. 11674874498878E+00 

-6. 14797666248879E-01 

i  RE-fINTG_l> 

IM-fINTG_l>  |  MATLAB  QUADGK  ALGORITHM  | 

CASE  35  SIDE  1  6. 14264876010087E+00 

-8.50317969612828E-01  6. 14264876010089E+00 

-8.50317969612829E-01 

SIDE  2  6.62366864094923E-01 

-5.74470195405343E-01  6.62366864094923E-01 

-5. 74470 195405343 E-01 

SIDE  3  5.89087476299632E+00 

-6. 14601954320508E-01  5.89087476299633E+00 

-6. 14601954320508E-01 

We  also  compute  according  to  Table  6.14  but  for  the  z'  values  we  use  for  the  surface  integral,  as 
in  Table  6.16.  We  note  here  that  if  the  OP  were  exactly  on  the  boundary,  i.e.,  z'  were  equal  to 
zero,  the  line  integral  in  question  would  not  exist.  Fortunately,  the  OP  cannot  be  on  the 
boundary;  however,  if  it  is  near  enough,  it  may  create  problems  for  any  quadrature  used  to 
evaluate  the  integral.  We  demonstrate  this  in  Table  6.19  with  Case  38.  In  this  case,  the  OP  is 
right  above  the  right-lower  corner.  We  see  that,  in  the  approach  to  the  comer,  the  GKQ  starts 
losing  SD.  This  is  true  for  the  real  part  of  the  integral,  the  one  that  is  singular.  Sides  1  and  3 
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(hypotenuse  and  base,  respectively)  have  this  corner  in  common.  While  for  Side  2  we  get  14  SD 
at  the  10'6  value,  for  the  other  two  sides  12  for  the  real  part  but  15  for  the  imaginary.  For 
completeness,  we  present  the  rest  of  the  cases  in  Appendix  H. 

Table  6.19:  The  line  integral  of  (6.3)  for  an  OP  above  the  lower  right-hand  comer 

of  the  triangle. 


Casae 

38 

PRESENT  METHOD 

MATLAB  G-K  QUADRATURE 

z' 

SIDE 

RE{INTG 1> 

IM{INTG 1} 

RE{INTG 1> 

IM{INTG 1> 

5.0 

S  1 

-2. 787788771 23972E-01 

1. 142805 17687205E-02 

-2. 787788771 23972E-01 

1. 14280517687205E-02 

S  2 

-1.94184778869090E-01 

1.60631 223094381 E-02 

-1.94184778869090E-01 

1.6063 1223094381  E-02 

S  3 

-1.98613359322592E-01 

4. 1 1400306068688E-03 

-1.98613359322592E-01 

4. 11400306068682E-03 

1.0 

S  1 

8. 1 1487798461838E-01 

-7.94679252600586E-01 

8. 1 148779846 1838E-01 

-7.94679252600586E-01 

S  2 

3.80119469546867E-01 

-5. 36269065701 85 IE-01 

3.80119469546868E-01 

-5. 3626906570 1851  E-01 

S  3 

6.64800195907522E-01 

-5.74698696029768E-01 

6.64800195907522E-01 

-5.74698696029768E-01 

0.01 

S  1 

5.4538328674847 1E+00 

-8.50504527263624E-01 

5.45383286748472E+00 

-8.50504527263625E-01 

S  2 

6. 64757229 191059E-01 

-5. 746947752  31 700  E-01 

6. 6475722919 1059E-01 

-5. 7469477523 1700E-01 

S  3 

5.201 1988933535 1E+00 

-6. 14696002448062E-01 

5.201 19889335352E+00 

-6.14696002448062E-01 

l.E-04 

S  1 

1.00590502389874E+01 

-8. 505 1022 19825 15E-01 

1.00590502389879E+01 

-8. 505 1022 19825 16E-01 

S  2 

6.64800191610610E-01 

-5.74698695637687E-01 

6.64800191610610E-01 

-5.74698695637687E-01 

S  3 

9.80640081789547E+00 

-6. 14700082256520E-01 

9.80640081789576E+00 

-6. 14700082256520E-01 

l.E-06 

S  1 

1.46642204342387E+01 

-8. 505 1022255 1989E-01 

1.4664220434281 2E+01 

-8.50510222551989E-01 

S  2 

6.64800195907092E-01 

-5.74698696029728E-01 

6.64800195907093E-01 

-5.74698696029729E-01 

S  3 

1.441 157101 16020E+01 

-6. 14700082664502E-01 

1. 4411 57 101163 11 E+01 

-6. 14700082664502E-01 
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SUMMARY  AND  CONCLUSIONS 

We  have  considered  the  problem  of  computing  the  inner  integral  of  the  MFIE  to  a  prescribed 
precision.  This  integral  is  given  by  (1.12).  In  evaluating  this  integral,  our  aim  has  been  to  take 
the  greatest  advantage  possible  of  the  results  we  obtained  in  reference  1.  In  Section  1,  we  re¬ 
write  the  integral  as  a  line  integral  over  the  boundary  of  the  integration  triangle  and  a  remaining 
surface  integral  over  the  triangle.  In  Section  3,  we  use  the  same  techniques  as  in  reference  1  to 
evaluate  the  line  integral  to  prescribed  precision.  In  Section  2,  we  show  that  the  remaining 
surface  integral  can  be  written  in  such  a  way  as  to  have  to  evaluate  a  solid-angle  integral  and  one 
that  we  have  already  evaluated  in  reference  1.  In  Section  4,  we  give  an  explicit  formula  for 
measuring  solid  angle.  It  may  be  the  most  compact  fonnula  for  the  purpose.  In  Section  5,  we 
quote  results  from  reference  1  regarding  the  remaining  integral  in  Section  2.  In  Section  6,  we 
provide  a  thorough  validation  of  our  results. 

The  principal  equations  are  the  following.  The  integral  under  consideration  is  given  by  (1.12).  It 
is  re-written  as  two  integrals  in  (1.16).  The  surface  integral  in  (1.16)  is  written  in  terms  of  a 
Maclaurin-series  expansion  in  (2.2).  The  integrals  in  this  expansion  are  given  by  the  iteration 
formulas  (2.6)  -  (2.7),  as  well  as  the  integral  (2.4).  Results  for  the  latter  are  given  in  Section  5,  as 
quoted  from  reference  1.  The  initial  iterate  is  the  solid  angle  integral,  an  explicit  formula  for 
which  is  given  in  Section  4.  The  line  integral  in  (1.16)  is  evaluated  in  Section  3.  The  key  integral 
is  given  by  (3.6).  It  is  expanded  in  a  Maclaurin  series  in  (3.8)  and  evaluated  iteratively. 

The  iterative  formulas  that  we  use  in  this  report  do  not  converge  for  all  values  of  the  OP.  In 
reference  1  we  showed  that  the  two  formulas  in  Section  5,  namely,  (5.4)  -  (5.6)  and  (5.9)  - 
(5.11)  (and  two  other,  similar  sets  further  down)  converge  for  the  OP  being  within  half-a- 
wavelength  from  the  triangle’s  centroid.  Using  the  approach  of  Appendix  B  in  reference  1,  we 
can  show  that  (2.7)  converges  for  observation  points  within  a  wavelength  from  the  centroid.  For 
the  iteration  formulas  (3.37)  and  (3.42),  we  reach  the  same  conclusion.  Thus,  if  the  observation 
point  is  within  half-a  wavelength  from  the  centroid,  our  scheme  converges.  Unfortunately,  time 
did  not  permit  us  to  determine  where  we  should  use  our  method  and  where  we  should  be  using 
quadratures  and  cubatures,  as  we  did  in  Part  II  of  reference  1.  This  part  remains  to  be  done. 
Although  the  task  of  accomplishing  this  is  straightforward,  it  is  a  time-consuming  one.  There 
remains  also  the  outer  integral,  in  the  present  case,  the  first  of  the  iterated  integrals  in  (1.10). 
This  integral  is  over  the  observation  triangle  and  involves  a  basis  function  and  the  inner  integral, 
the  one  we  evaluated  in  this  report.  The  outer  integral  (in  our  approach)  can  only  be  evaluated 
numerically.  It  would  be  interesting  to  conduct  a  study  as  in  Part  II  of  reference  1  to  see  what 
size  cubatures  are  required,  for  a  specified  accuracy. 

In  this  report  and  in  reference  1,  we  considered  the  case  of  perfect  conductors.  When  other 
materials  are  present,  then  a  pair  of  integral  equations  is  required  to  determine  the  unknown 
electric  and  magnetic  current  densities.  These  equations  involve  the  same  kernels  that  appear  in 
the  EFIE  and  the  MFIE.  If  the  materials  are  lossless,  then  the  present  analysis  and  that  in 
reference  1  apply  directly.  If  they  are  lossy,  then  the  wavenumber  becomes  complex  and  results 
not  only  in  an  oscillatory  motion,  as  we  have  encountered  to  this  point,  but  also  in  an 
exponentially  decaying  factor.  The  latter  can  also  be  replaced  by  its  Taylor-series  expansion, 
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multiplying  the  sine  and  cosine  expansions.  The  analysis  from  this  point  on  would  proceed  as  in 
the  present  report  and  reference  1 . 
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ACRONYMS 


EFIE 

Electric  Field  Integral  Equation 

GKQ 

Gauss-Kronrod  Quadrature. 

IM 

Impedance  Matrix 

IT 

Integration  Triangle 

MFIE 

Magnetic  Field  Integral  Equation 

MoM 

Method  of  Moments 

OP 

Observation  Point 

RWG 

Rao-Wilton-Glisson 
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APPENDIX  A 

ANOTHER  TESTING  OF  MFIE 

In  this  Appendix,  we  test  (1.7)  using  h'  x  fm(r')  as  the  testing  function.  We  thus  have 


1  N 

=  J  [A’xf.(r’)]- J“(r’)rfS' 

^  n= 1  j  Tm 

N 

~Yjan\ dS'[n'xim(r')]-h'x\ rfST^r)  x  Vg(r,r')  ,  r'  e  Tm  ,  m  =  1,2,3,...,  N 


n= 1  f 


Using  well-known  vector  identities,  we  can  re-write  this  as 

1  N 

^  n= 1  t  f 

N 

-Xan\dS'Ur')-\dSfn(r)xVg(r,r')  ,  r'eEm,  m  =  1,2,3,...,  A  . 


H=I  J 


The  integral  of  interest  is  the  same  as  in  ( 1 . 1 1 ). 


(A.  1) 


(A.2) 
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APPENDIX  B 

IMAGINARY  PART  OF  SURFACE  INTEGRAL  FOR  OBSERVATION  POINTS 
PROJECTING  ONTO  THE  INTEGRATION  TRIANGLE  AND  IN  CLOSE  PROXIMITY 

TO  IT 

In  (2.2),  we  consider  the  case  where  the  OP  is  above  the  triangle  and  in  close  proximity  to  it.  The 
imaginary  part  of  the  integral  is  given  by 


mfmh-AiJj 


frl. 

-i2nR 

Lj_| 

\d 

sin(2;rR) 

Ii& 

l  *  J 

An  \ 

\dz 

R 

dS. 


(B.l) 


We  perform  the  indicated  differentiation 

,  /  z'  r  2nRcos(2nR)-sin(2nR) 

Im{F(r')|  = - [ - - - { - - - ’-dS 

1  V  n  4 n\  R 

z'  r  2^R[l-cos(2^f?)]  +  [sin(2^f?)-2^R] 

An  ■* 


R 


z'  r  l-cos(2nR)  z'  r  2nR -sin(2^R) 


I 


R2 


dS-—\ 
An 


R 3 


-dS 


(B-2) 


where 


R  =  yj(x-x')2  +(y-y')2  +z'2  =  ^]p2  +z'2  ,  z' ±  0,  p  =  ^(x-x')2  +(y-y’Y  .  (B.3) 

We  expand  the  trigonometric  function  according  to  the  number  of  SDs  we  require  and  write 


lm{F(r-)} = cf  or  w  j R^ds  (-0"  w  , 

1  V  1  2tT  (2/i)!  J  2tT  (2/1  +  1)!  { 


(B.4) 


We  let 

N  =  max  (  A. ,  Ns ) 
and  redefine  (B.4)  as 


(B.5) 


{27r) f^-^z(-ir(2jrrf^ 

1  2tr  (2/1)!  J  2tr  (2/7  +  1)!  J 


J  N-\ 


-yK-rw 

«=1 


1 


1 


(2/i)!  (2/1  +  1)! 


jR2(”-1W5  =  z'X 

t  n= 1 


(2/7  +  1)! 


57 


APPENDIX  B 


NAWCADPAX/TR-20 1 2/200 


/f(-l)'(n  +  l){2x)2'+2 
tZ  (2/1  +  3)! 


We  use  (B.3)  to  re-write  this  as 

:(  —  !)> +  0(2*) 


N-2  i 


2«+2 


Im{F(r')}  =  z'£- 

n=0 


(2«  +  3)! 


-J(p2+z'2)"dS. 


(B.6) 


(B.7) 


We  expand  the  binomial 

Im{F(r')} - h1*" X(" )z«- \ P^dS 


n=0 


(2n  +  3)! 


m= 0  ' 


and  collect  terms  in  powers  of  z' 


2/2+2 


Im{F(r')}  =  z'X- 


n=0 


[2n  +  3)! 


-  J  /?2V5 


7V-2  i 


\  2/2+2 


+z" X  <  *)  ("  +  1)(2”f  j n ) j  2(,-ys 

+  (2«  +  3)! 


tV-2  i 


+z, H)> +1)(2*) 


2/2+2 


n=2 


(2n  +  3)! 


+z'2("-3>«  2  (  1)"("-|-1)(2^)'"'^/  n  \r p*~«»)ds 
„Vi  (2n  +  3)!  alA'-3/} 


+Z 


■2(7V-2)+l  (-1)“(A'-1)(2^) 

(27V  - 1)! 


2/2+2 


R 


Alternatively,  we  can  re-arrange  (B.8)  in  terms  of  equal  powers  of  p  to  obtain 

(-l)"(77  +  l)(2zrz')2 


ImjF(r')}  =  z'(2;r)~ 


\dS 


N-2  | 

z- 


\  2/2 


(2n  +  3)! 


+z'(2^-)" 


J  p2dS 


N-2 

z 


[n +l){2zrz’f" 


+z'(2t r)  JpVS'^]- 

t  n=2 


(2n  +  3)! 
(-l)>  +  l)(2/rz') 


n  r  -2 

n-1  f 


(2n  +  3)! 


n  ,-4 

n-2  f 


(B.8) 


(B.9) 
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+z'(2/r )2 

+zr(2^y 


jp2(N-3)dS 

T 

J+f+iS 


N-2 


y  (  l)  (/7  +  l)(2;rz)  |  n  2(.v  3j 


n=N~3 


(2/1  +  3)!  ^“3 

(-if  (A-1)(2^z')2(7V‘2) 


(2A-1)! 


Shifting  indices,  we  get 


(B.10) 


ImjF(r')}  =  z'(2/r)“ 


JV-2 

Z 

«=0 


(-l)>  +  l)(2/rz')2" 


-z'(2^-)4 


J/rdS 


7V-3 

z 


J  (2/7  +  3)! 

(-!)>  +  2)(n  +  l)(2/rz') 


+z  '(2^)‘Jp4<iS2- 

t  n= 0 


(2/7  +  5)! 

(-1)”  (/z  +  3)(/z  +2)(n  +  l)(2/rz') 


2(2/7  +  7)! 


-z'(2/r)“ 

+z'(2/r)2 


.r 

JV("-2)^ 


N -3 


f  (-1)  {n  +  N-2)(2nz')  (n  +  N.3 
h  (2/7  +  2X-3)! 
(-lf(A-l)(2^)2(W-2) 


(2A-1)! 


Making  the  leading  term  in  each  sum  equal  to  one,  we  get 


5! 

z'3(2/r)6 

7! 


\p2dS 

T 


T 
N- 3 

z 


N-2 

z 

n=0 


(-1)"3!(//  +  1)(2  /rz'f 


(2n  +  3)! 


(-1)"  5!(/7  +  2)(//  +  l)(2zrz')2" 

2  (2/7 +  5)! 

(-1)"  7!(//  +  3)(//  +2)(//  +l)(2zrz') 
31(2/7  +  7)! 


(B.ll) 


z\N-2){2tt)2 


(2N-3)! 


\P 


,2(JV-3) 


^(-1)"  (27V-3)!(/7  +  A-2)(2^z')"/„  +  jV_3 

(  a t  oWt..  t  ir  +i  I  N  ~3 


+z'(2/ry 


\p2(N~2)dS 


(N-2){2n  +  2N-3)\ 
{-\)N  [N -\){2tt)2(N'2) 


(2A-1)! 


(B-12) 
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The  leading  term  in  each  sum  is  equal  to  one.  The  general  term  of  the  first  sum  is  greater  than  the 
general  term  of  the  second  sum,  which  is  greater  than  that  of  the  third  sum,  and  so  on.  Thus,  if 
we  set  the  second  term  of  the  first  sum  to  be  less  than  0.5xl0'M  or  less  than  the  machine  epsilon, 
then  we  only  need  retain  the  first  term  in  each  sum.  For  each  case  we  get 


{Ircz'y 

10 


<0.5-10  M 


,  4  5  - 10  ^ 

z  < - 

2  n 


(B.13) 


and 


(2  nz'X  ,  VlO^ 

- - —  <  £  =>  Z  < - 

10  In 


(B.14) 


From  (reference  15),  we  have  obtained  Table  1  that  presents  the  IEEE  754-2008  standard  and 
gives  the  value  of  machine  epsilon  in  the  last  column.  For  values  of  z'  as  in  (B.13)  or  (B.  14),  we 
can  replace  the  sums  in  (B.12)  by  one  (their  leading  term).  This  corresponds  to  using  only  the 
first  term  in  (B.9) 


N-2 


Im{F(r')}  =  z'£ 


(-i)> +i)(2*r 


n= 0 


(2n  +  3)! 


-J  p2"dS 


(B.15) 


Table  B.  1 :  Machine  epsilon  for  IEEE  754-2008  standard. 


IEEE  754  - 
2008 

Common  Name 

Base 

b 

Fractional  digits 

P 

Machine  Epsilon 

b~p/2 

Value 

binary  16 

half  precision 

2 

10 

2  -11 

4.88e-04 

binarv32 

single  precision 

2 

23 

2 

5.96e-08 

binarv64 

double  precision 

2 

52 

2-53 

l.lle-16 

binarvl28 

quad(ruple) 

precision 

2 

112 

2-m 

9.63e-35 

We  turn  to  the  evaluation  of  the  integrals  in  (B.15).  The  integrands  are  polynomials  in  x  and  y  of 
degree  2 n.  Hence,  we  can  find  an  appropriate  cubature  that  integrates  them  exactly. 
Alternatively,  and  since 


1 


v'(p^>q^)[p'v"+^'vp 

To  A.  ^2 « — 1  .  s-s  2n  ~\  2 n 

1 2np  ■  pp  +2 p  \  =  p 


2  (n  + 1) 
1 


2{n  + 1) 


(B-16) 

(B-17) 
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we  can  use  the  divergence  theorem  in  two  dimensions  to  write 

1P2’JS  =  20TM ) | "  ■  -  2RI # ' ^  (B' 18> 

with  v  the  exterior  unit  normal  vector  to  the  triangle’s  side.  The  line  integrals  can  be  evaluated 
to  machine  precision  using  a  GKQ.  These  integrals  can  also  be  evaluated  analytically.  Using 
Figure  3.2  and  the  definitions  in  that  section,  we  first  write 

|  {v-p)p2nds  =  YJVj-\pP2''ds  (B.19) 

dT  j= 1  Sj 

and  then  evaluate 


■  |  p p“is  =  Vj  ■  j  (a;„  +  </){|a;„  f  +  2 tt  ■  a  Ms  +  s1  \  ds 


“( 

“( 

“( 


J 

+2iJ-*j«s+s 


w 1 


v .  -  a, 


J  J 

)  1 

U 

n 

)Z 


{'M’VVif} 


ds  =  ( 

ri  *=( 


,{(^  +  /.-ay+1)2+|a.+i|  - 


)/{( 

)zfc)(-ua,+1rmr;r%2^ 


(B.20) 

nn 

ds 


m= 0 


/i\(Vay+ 1) 


J  2m  + 1 


2(« 


\2m+l 


(B-21) 
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APPENDIX  C 

ON  THE  SOLID-ANGLE  SUBTENDED  BY  A  PLANAR  FIGURE  AT  A  POINT  IN  SPACE 


We  consider  a  planar  surface  T  and  a  point  in  space.  The  point  is  off  the  surface  and  we  call  it  the 
OP.  From  this  point,  we  draw  a  straight  line  to  a  point  on  the  boundary  of  T  and  we  let  the 
boundary  point  go  once  around  the  boundary.  This  action  generates  a  cone-like  surface  whose 
apex  is  the  OP  and  base  is  the  planar  surface  T.  With  center  the  OP,  we  draw  a  sphere  of  unit 
radius.  We  call  Ts  the  region  of  the  sphere  that  the  cone-like  structure  intercepts.  The  surface 
area  of  Ts  is  the  solid  angle  value  subtended  by  T  at  the  OP  and  is  a  positive  quantity  between  0 
and  4;rsteradians.  The  mathematical  fonnula  for  it  is 


r=1  r2  sin 0d0d(p  =  J  sin ddOdcp  =  ) 

Ts 


(Cl) 


where  A(TS )  stands  for  the  area  of  Ts.  The  second  integral  is  valid  for  any  radius  of  the  sphere 

that  is  less  than  the  minimum  distance  between  the  OP  and  the  boundary.  By  the  divergence 
theorem,  we  can  write  (C.  1)  as 


Q(r';T)  =  Jn-V 

T 


f  i  A 


dS 


(C.2) 


where  R  is  the  distance  from  the  OP  to  a  point  of  T  and  n  is  the  unit  normal  vector  on  T  pointing 
into  the  region  bounded  by  T,  Ts,  and  the  lateral  surface  of  the  cone;  thus,  h  has  a  positive 
component  in  the  direction  of  the  OP.  The  integral  in  (C.2)  can  be  used  for  any  position  of  the 
OP.  If  the  OP  is  reflected  into  the  half-space  away  from  which  the  nonnal  is  pointing,  then  the 
result  will  be  the  negative  of  the  solid  angle.  To  make  things  clear,  if  T  lies  on  the  xv-plane  and 
77  points  in  the  z-direction,  then 


n(r/;r)=Jf.v 


r  \ 

1 


r-r 


dS  =  -\z-W 


f  1  A 


VI1*-1*  ly 


dS,  z’>  0 


where  r’  is  the  image  of  r'  about  the  xv-plane  and  r'  =  x'x  +  y'y  +  z'z 


(C.3) 


63 


APPENDIX  C 


NAWCADPAX/TR-20 1 2/200 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


64 


NAWCADPAX/TR-20 1 2/200 


APPENDIX  D 

COMPARISON  OF  TWO  SOLID-ANGLE  FORMULAS  FOR  POINTS  NEAR  AND 

ON  A  TRIANGLE’S  BOUNDARY 

We  consider  here  the  triangle  of  Figure  6.1  and  the  OPs  of  Table  6.14  when  z'  is  equal  to  zero  or 
very  close  to  it.  At  these  OPs,  we  compute  the  normalized  solid  angle  Fonnula  1  (given  by  (4.1) 
and  (4.2))  and  Formula  2  (given  by  (4.1)  and  (4.3)).  We  display  the  results  in  Table  D.l  for  three 
values  of  z':  10’18,  Iff20  and  0  m.  We  display  both  SP  and  DP  results.  We  discuss  first  the  DP 
results. 

Case  36  represents  a  point  in  the  interior  of  the  triangle  and  we  see  that  both  formulas  provide 
the  correct  result  when  z'  is  equal  to  zero.  The  results  for  z'  equal  to  10"18  and  lO’20  are  also 
correct  to  DP. 

Case  40  represents  a  point  in  the  exterior  of  the  triangle  and  both  formulas  agree  for  all  three 
values  of  z'. 

Case  37  represents  the  left  vertex  of  the  base  of  the  triangle.  In  the  approach  to  it,  both  formulas 
provide  the  correct  answer  to  DP  accuracy.  When  z'  is  equal  to  zero,  however,  Formula  1  treats 
the  point  as  if  it  were  exterior  to  the  triangle,  while  Fonnula  2  sees  it  inside  the  triangle.  The 
correct  answer  is  0.25. 

Case  38  represents  the  right  vertex  of  the  base  of  the  triangle.  In  the  approach  to  it,  both  formulas 
provide  the  correct  answer  to  DP  accuracy.  When  z'  is  equal  to  zero,  however,  both  formulas  see 
a  point  outside  the  triangle.  The  correct  answer  is  0.125. 

Case  39  represents  the  middle  point  of  the  hypotenuse.  For  all  three  values  of  z',  Formula  1  sees 
a  point  exterior  to  the  triangle.  For  the  first  two  values  of  z',  Formula  2  provides  the  correct  value 
(both  numbers  are  close  to  0.5)  but  fails  when  z'  is  equal  to  zero  and  views  the  point  as  an 
exterior  point. 

What  is  happening  in  the  last  three  cases?  The  coordinates  of  the  vertex  in  Case  37  are  (-1/3,  - 
1/3,  0).  When  the  number  1/3  is  rounded  off  in  DP,  it  becomes  a  smaller  number  than  it  actually 
is;  thus,  the  point  (-1/3,  -1/3,  0)  in  Figure  6.1  will  move  to  the  right  and  up,  landing  inside  the 
triangle.  For  reasons,  we  cannot  explain,  both  fonnulas  see  the  desired  position  (-1/3,  -1/3,  0) 
when  z'  is  different  from  zero  but  Formula  1  sees  an  exterior  point  when  z'  is  equal  to  zero,  while 
Formula  2  sees  an  interior  point.  Since  in  the  computer  the  point  is  an  interior  point,  we  may 
claim  that  Formula  2  gives  a  better  answer  than  Formula  1  when  z'  is  zero. 

The  coordinates  of  the  vertex  for  Case  38  are  (2/3,  -1/3,  0).  When  2/3  is  rounded  off,  it  becomes 
larger  while  1/3  becomes  smaller;  thus,  the  point  are  (2/3,  -1/3,  0)  moves  to  the  right  and  up  in 
Figure  6.1,  landing  outside  the  triangle.  This  is  what  both  fonnulas  see  for  z'  equal  to  zero  and 
provide  the  correct  result.  For  z'  different  from  zero,  they  see  Vertex  2  in  Figure  1 . 


65 


APPENDIX  D 


NAWCADPAX/TR-20 1 2/200 


The  coordinates  of  the  point  in  Case  39  are  (1/6,  1/6,  0).  When  1/6  is  rounded  off,  it  becomes 
larger;  thus,  the  midpoint  of  the  hypotenuse  moves  to  the  right  and  up,  landing  outside  the 
triangle.  This  is  what  Formula  1  sees,  while  Formula  2  sees  the  original  point  for  z'  different 
from  zero  and  the  shifted  point  for  z'  equal  to  zero. 

Turning  to  SP,  we  see  that,  in  Cases  36  -  38,  SP  provides  the  same  kind  of  result  as  DP,  with 
Formula  2  performing  slightly  better  than  Formula  1.  In  Case  39,  Formula  1  behaves  as  for  DP, 
while  Formula  2  provides  the  correct  answer  for  points  of  the  triangle  but  only  to  four  SDs.  In 
Case  40,  both  formulas  behave  as  in  DP  but,  for  points  off  the  triangle,  they  both  provide  only 
six  SDs.  In  general,  we  would  advise  avoiding  SP  calculations  for  the  solid  angle. 
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Table  D.l:  Normalized  solid  angle  values  for  the  OPs  of  Table  6.14 


Case  36 

9 

10 

0 

z' 

1.00E-18 

1.00E-20 

0.00E+00 

SP 

Formula  1 

1.000000E+00 

1.000000E+00 

1.000000E+00 

Formula  2 

1.000000E+00 

1.000000E+00 

1.000000E+00 

DP 

Formula  1 

1.00000000000000E+00 

1.00000000000000E+00 

1.00000000000000E+00 

Formula  2 

1.00000000000000E+00 

1.00000000000000E+00 

1.00000000000000E+00 

Case  37 

9 

10 

z' 

1.00E-18 

1.00E-20 

0.00E+00 

SP 

Formula  1 

2.500000E-01 

2.500004E-01 

0.000000E+00 

Formula  2 

2.500000E-01 

2.500000E-01 

1.000000E+00 

DP 

Formula  1 

2.50000000000000E-01 

2.50000000000000E-01 

0.00000000000000E+00 

Formula  2 

2.50000000000000E-01 

2.50000000000000E-01 

1.00000000000000E+00 

Case  38 

9 

10 

z' 

1.00E-18 

1.00E-20 

0.00E+00 

SP 

Formula  1 

1.250000E-01 

1.250003E-01 

0.000000E+00 

Formula  2 

1.250000E-01 

1.250000E-01 

0.000000E+00 

DP 

Formula  1 

1.25000000000000E-01 

1.25000000000000E-01 

0.00000000000000E+00 

Formula  2 

1.25000000000000E-01 

1.25000000000000E-01 

0.00000000000000E+00 

Case  39 

9 

10 

z' 

1.00E-18 

1.00E-20 

0.00E+00 

SP 

Formula  1 

1.000000E+00 

1.000000E+00 

1.000000E+00 

Formula  2 

4.999827E-01 

4.999827E-01 

0.000000E+00 

DP 

Formula  1 

5. 73354097941042 E-03 

5.73416107719460E-05 

0.00000000000000E+00 

Formula  2 

4. 999999987585 19E-01 

4. 999999987585 19E-01 

0.00000000000000E+00 

Case  40 

9 

10 

z' 

1.00E-18 

1.00E-20 

0.00E+00 

SP 

Formula  1 

5.922377E-18 

5.922376E-20 

0.000000E+00 

Formula  2 

5.922380E-18 

5.922379E-20 

0.000000E+00 

DP 

Formula  1 

5. 92237767 136760E-18 

5. 92237767 136760E-20 

0.00000000000000E+00 

Formula  2 

5. 92237767 136760E-18 

5. 92237767 136760E-20 

0.00000000000000E+00 
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APPENDIX  E 

CONVERGENCE  OF  FORMULA  (2.2)  FOR  THE  SURFACE  INTEGRAL 

In  this  Appendix,  we  present  results  as  in  Table  6.16  for  the  remaining  cases  in  Table  6.14.  In 
Table  E.l,  we  display  the  convergence  of  our  algorithm  for  the  surface  integral  for  the  OP  of 
Case  37,  Table  6.14.  Convergence  is  to  15  SD  and  as  a  function  of  the  number  of  terms  in  the 
trigonometric  sums  in  (2.2).  The  entries  in  light  green  are  the  first  ones  for  which  our  formula 
has  stabilized  to  15  SD.  The  OP  is  above  the  right-angle  corner  of  the  IT. 

In  Table  E.2,  we  display  the  convergence  of  our  algorithm  for  the  surface  integral  for  the  OP  of 
Case  38,  Table  6.14.  Convergence  is  to  15  SD  and  as  a  function  of  the  number  of  terms  in  the 
trigonometric  sums  in  (2.2).  The  entries  in  light  green  are  the  first  one  for  which  our  formula  has 
stabilized  to  15  SD.  The  OP  is  above  the  lower  right-hand  corner  of  the  IT. 

In  Table  E.3,  we  display  the  convergence  of  our  algorithm  for  the  surface  integral  for  the  OP  of 
Case  39,  Table  6.14.  Convergence  is  to  15  SD  and  as  a  function  of  the  number  of  terms  in  the 
trigonometric  sums  in  (2.2).  The  entries  in  light  green  are  the  first  one  for  which  our  formula  has 
stabilized  to  15  SD.  The  OP  is  above  the  midpoint  of  the  hypotenuse  of  the  IT. 

In  Table  E.4,  we  display  the  convergence  of  our  algorithm  for  the  surface  integral  for  the  OP  of 
Case  40,  Table  6.14.  Convergence  is  to  15  SD  and  as  a  function  of  the  number  of  terms  in  the 
trigonometric  sums  in  (2.2).  The  entries  in  light  green  are  the  first  one  for  which  our  formula  has 
stabilized  to  15  SD.  The  OP  approaches  a  point  in  the  triangle's  plane  but  slightly  outside  the  IT. 
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Table  E.  1 :  Convergence  of  our  algorithm  for  Case  37. 


z'  = 

5.000 

4.500 

REAL{F(r')} 

IMAG{F(r’)} 

REAL{F(r')} 

IMAG{F(r')} 

15  SD,  Nc=Ns=4 

1.66229 178654614E-03 

4.90079446997064E-03 

2.67743268334130E-04 

5.78612816705256E-03 

15  SD,  Nc=Ns=5 

1.66229 178654613  E-03 

4.90079446997063E-03 

2.67743268334127E-04 

5.78612816705257E-03 

15  SD,  Nc=Ns=6 

1.66229 178654613  E-03 

4.90079446997064E-03 

2.67743268334123E-04 

5.78612816705257E-03 

15  SD,  Nc=Ns=7 

1.66229 178654613  E-03 

4.90079446997064E-03 

2.67743268334125E-04 

5.78612816705257E-03 

15  SD,  Nc=Ns=8 

1.66229 178654613  E-03 

4.90079446997064E-03 

2.67743268334125E-04 

5.78612816705257E-03 

15  SD,  Nc=Ns=9 

1.66229 178654613  E-03 

4.90079446997064E-03 

2.67743268334125E-04 

5.78612816705257E-03 

15  SD,  Nc=Ns=10 

1.66229 178654613  E-03 

4.90079446997064E-03 

2.67743268334125E-04 

5.78612816705257E-03 

15  SD,  Nc=Ns=ll 

1.66229 178654613  E-03 

4.90079446997064E-03 

2.67743268334125E-04 

5.78612816705257E-03 

15  SD,  Nc=Ns=12 

1.66229 178654613  E-03 

4.90079446997064E-03 

2.67743268334125E-04 

5.78612816705257E-03 

z'  = 

1.000 

0.010 

REAL{F(r')} 

IMAG{F(r’)} 

REAL{F(r')} 

IMAG{F(r’)} 

15  SD,  Nc=Ns=4 

-3. 304146503 10833E-02 

3.11999935023589E-03 

-1.23597765705000E-01 

3. 24684202 112266E-05 

15  SD,  Nc=Ns=5 

-3.30414650308276E-02 

3.11999935001453E-03 

-1.23597765854258E-01 

3. 24684605349 130E-05 

15  SD,  Nc=Ns=6 

-3.30414650308277E-02 

3.11999935001460E-03 

-1.235977658541 19  E-01 

3.24684604964339E-05 

15  SD,  Nc=Ns=7 

-3.30414650308277E-02 

3.11999935001460E-03 

-1.235977658541 19  E-01 

3. 246846049646 14E-05 

15  SD,  Nc=Ns=8 

-3.30414650308277E-02 

3.11999935001460E-03 

-1.235977658541 19  E-01 

3. 246846049646 14E-05 

15  SD,  Nc=Ns=9 

-3.30414650308277E-02 

3.11999935001460E-03 

-1.235977658541 19  E-01 

3. 246846049646 14E-05 

15  SD,  Nc=Ns=10 

-3.30414650308277E-02 

3.11999935001460E-03 

-1.235977658541 19  E-01 

3. 246846049646 14E-05 

15  SD,  Nc=Ns=ll 

-3.30414650308277E-02 

3.11999935001460E-03 

-1.235977658541 19  E-01 

3. 246846049646 14E-05 

15  SD,  Nc=Ns=12 

-3.30414650308277E-02 

3.11999935001460E-03 

-1.235977658541 19  E-01 

3. 246846049646 14E-05 

z'  = 

0.0001 

0.000001 

REAL{F(r')} 

IMAG{F(r’)} 

REAL{F(r')} 

IMAG{F(r')} 

15  SD,  Nc=Ns=4 

-1.24986000435264E-01 

3.24638198652974E-07 

-1. 24999859824252  E-01 

3.19908567264832E-09 

15  SD,  Nc=Ns=5 

-1. 2498600061237  IE-01 

3.24685938990200E-07 

-1. 24999859824252  E-01 

3.19908567264832E-09 

15  SD,  Nc=Ns=6 

-1. 24986000612 198E-01 

3.24685891449776E-07 

-1. 24999860001661  E-01 

3.24690652364757E-09 

15  SD,  Nc=Ns=7 

-1. 24986000612 199E-01 

3.24685891484467E-07 

-1.24999860001488E-01 

3.24685888125090E-09 

15  SD,  Nc=Ns=8 

-1. 24986000612 199E-01 

3.24685891484467E-07 

-1.24999860001488E-01 

3. 2468589 1594534E-09 

15  SD,  Nc=Ns=9 

-1. 24986000612 199E-01 

3.24685891484467E-07 

-1.24999860001488E-01 

3. 2468589 1594534E-09 

15  SD,  Nc=Ns=10 

-1. 24986000612 199E-01 

3.24685891484467E-07 

-1.24999860001488E-01 

3. 2468589 1594534E-09 

15  SD,  Nc=Ns=ll 

-1.249860006 12 199E-01 

3.24685891484467E-07 

-1.24999860001488E-01 

3. 2468589 1594534E-09 

15  SD,  Nc=Ns=12 

-1. 249860006 12 199E-01 

3.24685891484467E-07 

-1.24999860001488E-01 

3. 2468589 1594534E-09 
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Table  E.2:  Convergence  of  our  algorithm  for  Case  38. 


z'  = 

5.000 

4.500 

REAL{F(r’)} 

IMAG{F(r')} 

REAL{F(r')} 

IMAG{F(r')} 

15  SD,  Nc=Ns=4 

1. 72890617968269  E-03 

4. 803 16359879893  E-03 

3. 77040226833441 E-04 

5.68301993259905E-03 

15  SD,  Nc=Ns=5 

1.72890617968236E-03 

4. 803 16359879895  E-03 

3.77040226832750E-04 

5. 6830199325989  IE-03 

15  SD,  Nc=Ns=6 

1.72890617968236E-03 

4. 803 16359879895 E-03 

3.77040226832730E-04 

5.68301993259888E-03 

15  SD,  Nc=Ns=7 

1.72890617968236E-03 

4. 803 16359879895  E-03 

3.77040226832735E-04 

5.68301993259889E-03 

15  SD,  Nc=Ns=8 

1.72890617968236E-03 

4. 803 16359879895  E-03 

3.77040226832734E-04 

5.68301993259889E-03 

15  SD,  Nc=Ns=9 

1.72890617968236E-03 

4. 803 16359879895 E-03 

3.77040226832735E-04 

5.68301993259889E-03 

15  SD,  Nc=Ns=10 

1.72890617968236E-03 

4. 803 16359879895  E-03 

3.77040226832735E-04 

5.68301993259889E-03 

15  SD,  Nc=Ns=ll 

1.72890617968236E-03 

4. 803 16359879895  E-03 

3.77040226832735E-04 

5.68301993259889E-03 

15  SD,  Nc=Ns=12 

1.72890617968236E-03 

4. 803 16359879895  E-03 

3.77040226832735E-04 

5.68301993259889E-03 

z'  = 

1.000 

0.010 

REAL{F(r’)} 

IMAG{F(r')} 

REAL{F(r')} 

IMAG{F(r')} 

15  SD,  Nc=Ns=4 

-2. 61083307157073  E-02 

3. 07874207729961 E-03 

-6.20687 139259531  E-02 

3.20426950357757E-05 

15  SD,  Nc=Ns=5 

-2. 61083307 127430E-02 

3.07874207412868E-03 

-6. 20687 168746579  E-02 

3. 204400139 15789E-05 

15  SD,  Nc=Ns=6 

-2. 61083307 127456E-02 

3.07874207413142E-03 

-6.20687168676938E-02 

3.20439982283057E-05 

15  SD,  Nc=Ns=7 

-2.61083307127456E-02 

3.07874207413142E-03 

-6. 20687 168677051  E-02 

3.20439982334936E-05 

15  SD,  Nc=Ns=8 

-2.61083307127456E-02 

3.07874207413142E-03 

-6.20687168677050E-02 

3.20439982334867E-05 

15  SD,  Nc=Ns=9 

-2. 61083307 127456E-02 

3.07874207413142E-03 

-6.20687168677050E-02 

3.20439982334867E-05 

15  SD,  Nc=Ns=10 

-2. 61083307 127456E-02 

3. 07874207413142 E-03 

-6.20687168677050E-02 

3.20439982334867E-05 

15  SD,  Nc=Ns=ll 

-2. 61083307 127456E-02 

3.07874207413142E-03 

-6.20687168677050E-02 

3.20439982334867E-05 

15  SD,  Nc=Ns=12 

-2.61083307127456E-02 

3. 07874207413142  E-03 

-6.20687168677050E-02 

3.20439982334867E-05 

z'  = 

0.0001 

0.000001 

REAL{F(r’)} 

IMAG{F(r')} 

REAL{F(r')} 

IMAG{F(r')} 

15  SD,  Nc=Ns=4 

-6. 2495695825592 1  E-02 

3.18993899829706E-07 

-6.24999537153430E-02 

1.75553548874080E-09 

15  SD,  Nc=Ns=5 

-6. 2495699 1069065  E-02 

3. 20444859397371 E-07 

-6.24999570001650E-02 

3.20801988734267E-09 

15  SD,  Nc=Ns=6 

-6. 24956990989495  E-02 

3. 2044125075905 IE-07 

-6. 2499956992 1975  E-02 

3. 20440649 164474E-09 

15  SD,  Nc=Ns=7 

-6.24956990989627E-02 

3.20441256837172E-07 

-6. 24999569922 107E-02 

3.20441258052063E-09 

15  SD,  Nc=Ns=8 

-6.24956990989627E-02 

3. 20441256830233  E-07 

-6. 24999569922 107  E-02 

3. 20441257358 174E-09 

15  SD,  Nc=Ns=9 

-6.24956990989627E-02 

3.20441256830233  E-07 

-6. 24999569922 107E-02 

3.20441257358174E-09 

15  SD,  Nc=Ns=10 

-6.24956990989627E-02 

3.20441256830233  E-07 

-6. 24999569922 107E-02 

3.20441257358174E-09 

15  SD,  Nc=Ns=ll 

-6.24956990989627E-02 

3.20441256830233  E-07 

-6. 24999569922 107E-02 

3.20441257358174E-09 

15  SD,  Nc=Ns=12 

-6.24956990989627E-02 

3.20441256830233  E-07 

-6. 24999569922 107  E-02 

3. 20441 257358 174E-09 
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Table  E.3:  Convergence  of  our  algorithm  for  Case  39. 


z'  = 

5.000 

4.500 

REAL{F(r’)} 

IMAG{F(r')} 

REAL{F(r')} 

IMAG{F(r')} 

15  SD,  Nc=Ns=4 

1. 62745 189256200E-03 

4.95019906577087E-03 

2. 106867572 12 165E-04 

5.83827756766792E-03 

15  SD,  Nc=Ns=5 

1. 62745 189256203  E-03 

4.95019906577090E-03 

2. 106867572 122 17E-04 

5.83827756766799E-03 

15  SD,  Nc=Ns=6 

1. 62745 189256201 E-03 

4.95019906577089E-03 

2. 106867572 12 195E-04 

5.83827756766797E-03 

15  SD,  Nc=Ns=7 

1.62745 189256202  E-03 

4.95019906577089E-03 

2. 106867572 12200E-04 

5.83827756766797E-03 

15  SD,  Nc=Ns=8 

1.62745 189256202  E-03 

4.95019906577089E-03 

2. 106867572 12 199E-04 

5.83827756766797E-03 

15  SD,  Nc=Ns=9 

1.62745 189256202  E-03 

4.95019906577089E-03 

2. 106867572 12200E-04 

5.83827756766797E-03 

15  SD,  Nc=Ns=10 

1.62745 189256202  E-03 

4.95019906577089E-03 

2. 106867572 12200E-04 

5.83827756766797E-03 

15  SD,  Nc=Ns=ll 

1.62745 189256202  E-03 

4.95019906577089E-03 

2. 106867572 12200E-04 

5.83827756766797E-03 

15  SD,  Nc=Ns=12 

1.62745 189256202  E-03 

4.95019906577089E-03 

2. 106867572 12200E-04 

5.83827756766797E-03 

z'  = 

1.000 

0.010 

REAL{F(r’)} 

IMAG{F(r')} 

REAL{F(r')} 

IMAG{F(r')} 

15  SD,  Nc=Ns=4 

-3.85251 154820082  E-02 

3.14083404792880E-03 

-2. 4576821601 1015  E-01 

3. 26827981 196468E-05 

15  SD,  Nc=Ns=5 

-3.8525 11548 19883  E-02 

3.14083404791340E-03 

-2.4576821601 13 19  E-01 

3. 2682798128 1253E-05 

15  SD,  Nc=Ns=6 

-3.8525 11548 19883  E-02 

3. 1408340479 1340E-03 

-2.457682 1601 13 19  E-01 

3. 2682798128 1800E-05 

15  SD,  Nc=Ns=7 

-3.8525 11548 19883  E-02 

3.14083404791340E-03 

-2. 4576821601 13 19  E-01 

3. 2682798128 1800E-05 

15  SD,  Nc=Ns=8 

-3.8525 1154819883  E-02 

3.14083404791340E-03 

-2. 4576821601 13 19  E-01 

3. 2682798128 1800E-05 

15  SD,  Nc=Ns=9 

-3.8525 11548 19883  E-02 

3. 1408340479 1340E-03 

-2.4576821601 13 19  E-01 

3. 2682798128 1800E-05 

15  SD,  Nc=Ns=10 

-3.8525 11548 19883  E-02 

3. 1408340479 1340E-03 

-2. 457682 1601 13 19  E-01 

3. 2682798128 1800E-05 

15  SD,  Nc=Ns=ll 

-3.8525 11548 19883  E-02 

3. 1408340479 1340E-03 

-2. 457682 1601 13 19  E-01 

3. 2682798128 1800E-05 

15  SD,  Nc=Ns=12 

-3.8525 11548 19883  E-02 

3.14083404791340E-03 

-2. 4576821601 13 19  E-01 

3. 2682798128 1800E-05 

z'  = 

0.0001 

0.000001 

REAL{F(r’)} 

IMAG{F(r')} 

REAL{F(r')} 

IMAG{F(r')} 

15  SD,  Nc=Ns=4 

-2.49957723082499E-01 

3.26829087442671E-07 

-2.49999577236076E-01 

3. 268 104459 13695E-09 

15  SD,  Nc=Ns=5 

-2.49957723083909E-01 

3.26829273902691E-07 

-2. 49999577237499  E-01 

3.26829278735652E-09 

15  SD,  Nc=Ns=6 

-2.49957723083909E-01 

3.26829273854127E-07 

-2. 49999577237499  E-01 

3.26829273878434E-09 

15  SD,  Nc=Ns=7 

-2.49957723083909E-01 

3.268292738541 27E-07 

-2. 49999577237499  E-01 

3.26829273878434E-09 

15  SD,  Nc=Ns=8 

-2.49957723083909E-01 

3.268292738541 27E-07 

-2. 49999577237499  E-01 

3.26829273878434E-09 

15  SD,  Nc=Ns=9 

-2.49957723083909E-01 

3.26829273854127E-07 

-2. 49999577237499  E-01 

3.26829273878434E-09 

15  SD,  Nc=Ns=10 

-2.49957723083909E-01 

3.26829273854127E-07 

-2. 49999577237499  E-01 

3.26829273878434E-09 

15  SD,  Nc=Ns=ll 

-2.49957723083909E-01 

3.26829273854127E-07 

-2. 49999577237499  E-01 

3.26829273878434E-09 

15  SD,  Nc=Ns=12 

-2.49957723083909E-01 

3.268292738541 27E-07 

-2. 49999577237499  E-01 

3.26829273878434E-09 
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Table  E.4:  Convergence  of  our  algorithm  for  Case  40. 


z'  = 

5.000 

4.500 

REAL{F(r’)} 

IMAG{F(r')} 

REAL{F(r')} 

IMAG{F(r')} 

15  SD,  Nc=Ns=4 

1. 63265 177430638E-03 

4.94291929209383  E-03 

2. 19 19745  27066 13E-04 

5.83059487773860E-03 

15  SD,  Nc=Ns=5 

1. 63265 177430639  E-03 

4.94291929209384E-03 

2.19197452706570E-04 

5.83059487773854E-03 

15  SD,  Nc=Ns=6 

1. 63265 177430638E-03 

4.94291929209383  E-03 

2.19197452706587E-04 

5.83059487773856E-03 

15  SD,  Nc=Ns=7 

1. 63265 177430638E-03 

4.94291929209383  E-03 

2.19197452706583E-04 

5.83059487773856E-03 

15  SD,  Nc=Ns=8 

1. 63265 177430638E-03 

4.94291929209383  E-03 

2.19197452706583E-04 

5.83059487773856E-03 

15  SD,  Nc=Ns=9 

1. 63265 177430638E-03 

4.94291929209383  E-03 

2.19197452706583E-04 

5.83059487773856E-03 

15  SD,  Nc=Ns=10 

1. 63265 177430638E-03 

4.94291929209383  E-03 

2.19197452706583E-04 

5.83059487773856E-03 

15  SD,  Nc=Ns=ll 

1. 63265 177430638E-03 

4.94291929209383  E-03 

2.19197452706583E-04 

5.83059487773856E-03 

15  SD,  Nc=Ns=12 

1. 63265 177430638E-03 

4.94291929209383  E-03 

2.19197452706583E-04 

5.83059487773856E-03 

z'  = 

1.000 

0.010 

REAL{F(r’)} 

IMAG{F(r')} 

REAL{F(r')} 

IMAG{F(r')} 

15  SD,  Nc=Ns=4 

-3.75158019087497E-02 

3. 137766491553 15 E-03 

-2.93455187094148E-02 

3. 2651241384150  IE-05 

15  SD,  Nc=Ns=5 

-3.75158019087259E-02 

3. 13776649153445 E-03 

-2. 93455 187090839  E-02 

3.26512413098356E-05 

15  SD,  Nc=Ns=6 

-3.75158019087259E-02 

3.1377 6649 153445  E-03 

-2. 93455 187090841 E-02 

3.26512413098946E-05 

15  SD,  Nc=Ns=7 

-3.75158019087259E-02 

3. 13776649153445 E-03 

-2. 93455 187090841  E-02 

3.26512413098946E-05 

15  SD,  Nc=Ns=8 

-3.75158019087259E-02 

3. 13776649 153445 E-03 

-2. 93455 187090841  E-02 

3.26512413098946E-05 

15  SD,  Nc=Ns=9 

-3.75158019087259E-02 

3.1377 6649 153445  E-03 

-2. 93455 187090841  E-02 

3.26512413098946E-05 

15  SD,  Nc=Ns=10 

-3.75158019087259E-02 

3. 13776649 153445 E-03 

-2. 93455 187090841  E-02 

3.26512413098946E-05 

15  SD,  Nc=Ns=ll 

-3.75158019087259E-02 

3. 13776649 153445 E-03 

-2. 93455 187090841  E-02 

3.26512413098946E-05 

15  SD,  Nc=Ns=12 

-3.75158019087259E-02 

3. 13776649 153445 E-03 

-2. 93455 187090841  E-02 

3.26512413098946E-05 

z'  = 

0.0001 

0.000001 

REAL{F(r’)} 

IMAG{F(r')} 

REAL{F(r')} 

IMAG{F(r')} 

15  SD,  Nc=Ns=4 

-2.98383298898281  E-04 

3. 265 13705511685 E-07 

-2.98383804959292  E-06 

3.26513705640839E-09 

15  SD,  Nc=Ns=5 

-2.98383298895017E-04 

3. 265 1370477571  IE-07 

-2.98383804956028E-06 

3.26513704904855E-09 

15  SD,  Nc=Ns=6 

-2.98383298895019E-04 

3.26513704776294E-07 

-2. 98383804956031  E-06 

3.26513704905455E-09 

15  SD,  Nc=Ns=7 

-2.98383298895019E-04 

3.26513704776294E-07 

-2. 98383804956031  E-06 

3.26513704905455E-09 

15  SD,  Nc=Ns=8 

-2.98383298895019E-04 

3.26513704776294E-07 

-2. 98383804956031  E-06 

3.26513704905455E-09 

15  SD,  Nc=Ns=9 

-2.98383298895019E-04 

3.26513704776294E-07 

-2. 98383804956031  E-06 

3.26513704905455E-09 

15  SD,  Nc=Ns=10 

-2.98383298895019E-04 

3.26513704776294E-07 

-2. 98383804956031  E-06 

3.26513704905455E-09 

15  SD,  Nc=Ns=ll 

-2.98383298895019E-04 

3.26513704776294E-07 

-2. 98383804956031  E-06 

3.26513704905455E-09 

15  SD,  Nc=Ns=12 

-2.98383298895019E-04 

3.26513704776294E-07 

-2. 98383804956031  E-06 

3.26513704905455E-09 
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APPENDIX  F 

CUBATURE  COMPUTATION  OF  DOUBLE-LAYER  POTENTIAL 
We  examine  the  calculation  of  the  double-layer  potential 


1  r  d  f  e 


F(r')  =  ~—  f- 

Arr  J  r)' 7 


-ikR  \ 


\n '  dz  y  R 


dS. 


(F.l) 


The  expression  as  is  presents  problems  due  to  the  presence  of  the  normal  derivative.  We  re-write 
it  in  the  form 


rd 

'  -Is 

dS- 

1 

1  <0 

— >  ^ 

{  R  J 

4  n 

*  dz  \  R  j 

dS. 


(F.2) 


The  second  term  on  the  right  is  related  to  the  solid-angle  integral  and  we  have  calculated  that 
elsewhere.  We  write  for  the  first  term 


1  I 

f5 

d?  ~ _ L 1 

f5 

cos(kR)  -1  -  z'sin(kR) 

4  n  5 

1  dz 

/  R  ) 

uo  —  - 

4  n\ 

1  dz 

R 

dS 


1  I 

r 8 

cos(kR)-l 

dS  +  —  I 

f5 

sin  (kR) 

1  dz 

R 

4  n  * 

\dz 

R 

dS. 


(F.3) 


Performing  the  differentiations,  we  get 


d 

cos  (kR)  - 1 

-kRsin(kR)-  cos(kR)  -  l]  dR  | 

dz 

R 

R2  dz 

z= 0 


? ,  (kR)sin(kR)  +  cos(kR)  - 1]  ^ z'^ 

=k'  UJ 


(F.4) 


and 


d 

sin  (kR) 

(kR) cos  (kR)  -  sin  (kR)  dR  ,  2  (kR)  cos  (kR)  -  sin(kR) , 

fzO 

dz 

R 

R2  dzU=°  '  (kR)2  1 

(F.5) 


In  (F.4)  we  expand  the  trigonometric  functions  in  power  series  and  we  keep  enough  terms  to 
make  the  remainder  smaller  than  a  prescribed  value 


d_ 

dz 


cos (kR)  - 1 
R 


( z'\  U-i 


2«-2 


I  n=0 


(2/i  +  l)! 


(2  n)\ 


(F.6) 
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We  let  Abe  the  greater  of  L  and  M,  and  write 


3  [" cos(A7?)-l]  (z'\2  pA  lV,  ( kR)2"  ^  ( kR)2n 

&L — ^ — J = UJ  (^fj 

= "- V  few-  («r  [  1  -  1  1  -  (-1)"  u J 

UJ  p  [(2/7  +1)!  (2n  +  2)!j  V  (2A-1)! 

=  (Ae  \NY2(_lY  (2n  +  l)(kRf  _  N  (kRf"-l)' 

UJ  ’  (2/7  +  2)!  (  (2A-1)! 


We  repeat  this  process  for  (F.5) 


3  sin  (kR) 

dz  R 


£V  vV  n-ME  y,  ,)■(</ 

R)  W  j  (2,,)!  h(  ’  (2n  +  l)! 


-  2 AVyv  ifAgT!  .2ZVY(  ,r  (” +  ‘)  W:" 

UJ*  Sf  (2n  +  l)!  h  ’  (2n  +  3)! 


We  note  that  this  expression  is  well  behaved  for  small  kR\  in  fact,  if  we  plot  a  normalized  form 
of  (F.5) 


,  .  13  sin  (kR)  (kR)cos(kR)-sin[kR) 

f(kR)= tv  &[++ J = 


we  get  a  rather  slowly  varying  graph,  even  around  the  origin  (see  Figure  F.l).  A  low-order 
cubature  should  then  work  very  well  on  this  function. 

This  is  no  longer  true  for  (F.7)  when  kR  is  small,  as  can  be  seen  in  Figure  F.2  where  we  plot  the 
function  g 


1  3  rcos(£f?)-ll  (kR)sm(kR)  +  cos(W)-l] 

g(tR) = \ =  jj^f 


(F. 10) 


We  return  to  (F.4)  and  write 


3  fcos(M)-ll  _  k2  (  z'\  2  J(M)sin(£fl)  +  [cos(£fl)-l]  l\( z' 


dz  R  2  i  R  j 


2  [R 


P  (  7'\ 

=  ~  ~  +k 2 

2  UJ 


2(kf?)sin(kf?)  +  2  cos(ki?)-l  ~(kR)2  f z' 
2  (kR)2  |U 


(F.l  1) 
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Figure  F.  1 :  The  graph  of  the  function/in  (F.9). 


Figure  F.2:  The  graph  of  the  function  g  in  (F.  10). 
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If  we  employ  (F.7)  in  (F.  1 1),  we  get 


d 

cos(ki?)-l 

k2 

rZ') 

fz1 

k2. 

A2 ( ,r(2»+i mf’  ,  ./.(“A”' 

dz 

R 

2 

UJ 

UJ 

lx  * 

’  (2n  +  2)!  1  ’  (2A-1)! 

k2  fz'/ 


N- 3 


\2«+l 


-  a'  \  y  (-o'  <2”+3)(M)'  +  (_n' 

’  (2/1  +  4)!  '  ’ 


(kR) 


2N-3 


(2N-1)\ 


(F-12) 


which  clearly  shows  that  the  expression  in  braces  is  well  behaved  for  all  kR.  We  can  also 
demonstrate  this  by  looking  at  the  graph  of  the  function 


i  ^  ^  2(ki?)sin(ki?)  +  2[cos(ki?)-l]-(kf?)2 
hkR)  =  2  (kRf 


(F.13) 


in  Figure  F.3. 


Figure  F.3:  The  graph  of  the  function  h  in  (F.13). 
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We  return  to  (F.3)  and  substitute  (F.5)  and  (F.  1 1)  in  it 


1 1 

f5 

d?  ~ _ L  i 

f5 

cos(kR)-l 

dS  +  - —  1 

f5 

sin  (kR) 

4  n  \ 

1  dz 

/  R  ) 

uO  —  - 

4  n\ 

1  dz 

R 

4  n\ 

1  dz 

R 

k2z'  c  1  k\'  r 2(kR)sm(kR)  +  2[cos(kR)-l]-(kRy 

= - — db - - - — ^ - - - —  dS 

8  7T  \r  8  n  \  ( kR ) 

ik3z'  (•  (kR)  cos  (kR)  -  sin  (kR) 

4/r  J  (kR) 


The  last  two  integrals  on  the  right  may  be  computed  using  cubatures.  The  first  integral  is 
identical  to  Q 0  in  (2.4).  It  has  been  thoroughly  studied  in  Section  2  of  reference  1.  The 
appropriate  integral  there  is  given  by  the  second  one  in  (1.21),  p.  9,  with  n  equal  to  zero.  This  is 
transformed  to  three  line  integrals  in  (2.1),  p.  1 1.  These  line  integrals  are  subsequently  evaluated 
starting  with  Subsection  2.2,  p.  17. 
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APPENDIX  G 

CUBATURE  RESULTS  FOR  THE  CASES  OF  TABLE  6.14 

In  this  Appendix,  we  complete  the  presentation  of  the  data  for  the  cubature  calculations  for  the 
points  of  Table  6.14.  The  results  for  Case  36  are  displayed  in  Table  6.17.  Here,  we  use  the  same 
format  to  display  the  results  for  the  rest  of  the  cases.  The  conclusions  are  summarized  in  Figures 
6.4  to  6.7. 
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Table  G.  1  .a:  Comparison  of  our  results  for  the  surface  integral  of  Case  37  with  those 


obtained  through  a  Riemann  sum  and  cubatures. 

z'  = 

5.0 

4.5 

15  SD,  Nc=Ns=8 

MFIE  METHOD 

REAL{F(r')} 

1.66229178654613E-03 

2.67743268334125E-04 

IMAG{F(r')} 

4.90079446997064E-03 

5.78612816705257E-03 

RIEMANN  SUM 

REAL{F(r')} 

1.66229178673289E-03 

2.67743268658897E-04 

N  =  10,000 

IMAG{F(r')} 

4.90079446948171E-03 

5.78612816653607E-03 

CUBATURE  DIM 

SD 

RE{F(r')} 

21 

1.66229178654563E-03 

2. 6774326833307  IE-04 

28 

15 

1.66229178654613E-03 

2.67743268334128E-04 

36 

1.66229178654613E-03 

2.67743268334128E-04 

45 

1.66229178654614E-03 

2.67743268334133E-04 

46 

1.66229178654613E-03 

2.67743268334126E-04 

54 

1.66229178654613E-03 

2.67743268334126E-04 

55 

1.66229178654613E-03 

2.67743268334128E-04 

66 

1.66229178654613E-03 

2.67743268334130E-04 

78 

1.66229178654613E-03 

2.67743268334130E-04 

85 

1.66229178654613E-03 

2.67743268334133E-04 

91 

1.66229178654613E-03 

2.67743268334135E-04 

105 

1.66229178654613E-03 

2.67743268334133E-04 

120 

14 

1.66229178654613E-03 

2.67743268334124E-04 

126 

1.66229178654612E-03 

2.67743268334130E-04 

CUBATURE  DIM 

SD 

IM{F(r')} 

21 

15,  15 

4.90079446997064E-03 

5. 786128 16705257E-03 

28 

4.90079446997064E-03 

5. 786128 16705257E-03 

36 

4.90079446997063E-03 

5.78612816705256E-03 

45 

4.90079446997065E-03 

5. 786128 16705259E-03 

46 

4.90079446997063E-03 

5. 786128 16705256E-03 

54 

4.90079446997063E-03 

5. 786128 16705256E-03 

55 

4.90079446997064E-03 

5. 786128 16705257E-03 

66 

4.90079446997063E-03 

5.78612816705257E-03 

78 

4.90079446997063E-03 

5. 786128 16705257E-03 

85 

4.90079446997063E-03 

5. 786128 16705256E-03 

91 

4.90079446997063E-03 

5.78612816705257E-03 

105 

4.90079446997063E-03 

5. 786128 16705257E-03 

120 

4.90079446997064E-03 

5. 786128 16705257E-03 

126 

4.90079446997063E-03 

5. 786128 16705256E-03 
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Table  G.  1  .b:  Comparison  of  our  results  for  the  surface  integral  of  Case  37  with  those 
obtained  through  a  Riemann  sum  and  cubatures. 


z'  = 

15  SD,  Nc=Ns=8 

MFIE  METHOD 

1.0E+00 

1.0E-02 

REAL{F(r')} 

IMAG{F(r')} 

-3.30414650308277E-02 

3.11999935001460E-03 

-1.23597765854119E-01 

3.24684604964614E-05 

RIEMANN  SUM 

REAL{F(r')} 

-3.30414650272966E-02 

-1.23597765853978E-01 

N  =  10,000 

IMAG{F(r')} 

3. 11999934980801 E-03 

3.24684604943186E-05 

CUBATURE  DIM 

SD 

RE{F(r')} 

21 

-3.30414650438932E-02 

-1.23597766234781E-01 

28 

-3.30414650318885E-02 

-1.23597765889297E-01 

36 

-3. 3041465030853  IE-02 

-1.23597765930566E-01 

45 

-3.30414650308289E-02 

-1.23597765825503E-01 

46 

-3.30414650308252E-02 

-1.23597765902280E-01 

54 

-3.30414650308280E-02 

-1.23597765892955E-01 

55 

-3.30414650308278E-02 

-1.23597765841787E-01 

66 

15 

-3.30414650308277E-02 

-1.23597765854805E-01 

78 

-3.30414650308277E-02 

-1.23597765862033E-01 

85 

-3.30414650308277E-02 

-1.23597765862308E-01 

91 

-3.30414650308277E-02 

-1.23597765848987E-01 

105 

-3.30414650308277E-02 

-1.23597765857442E-01 

120 

-3.30414650308277E-02 

-1.23597765854512E-01 

126 

12 

-3.30414650308277E-02 

-1.23597765854274E-01 

CUBATURE  DIM 

SD 

IM{F(r')} 

21 

3. 11999935001461  E-03 

3.24684604964487E-05 

28 

3. 11999935001461  E-03 

3.24684604964487E-05 

36 

15 

3.11999935001460E-03 

3.24684604964486E-05 

45 

11 

3. 11999935001461  E-03 

3.24684604964489E-05 

46 

3.11999935001460E-03 

3.24684604964486E-05 

54 

3.11999935001460E-03 

3.24684604964487E-05 

55 

3. 11999935001461  E-03 

3.24684604964487E-05 

66 

3.11999935001460E-03 

3.24684604964488E-05 

78 

3.11999935001460E-03 

3.24684604964487E-05 

85 

3.11999935001460E-03 

3.24684604964487E-05 

91 

3. 11999935001461  E-03 

3.24684604964488E-05 

105 

3.11999935001460E-03 

3.24684604964487E-05 

120 

3. 11999935001461  E-03 

3.24684604964486E-05 

126 

3.11999935001460E-03 

3.24684604964487E-05 
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Table  G.l.c:  Comparison  of  our  results  for  the  surface  integral  of  Case  37  with  those 
obtained  through  a  Riemann  sum  and  cubatures. 


z'  = 

15  SD,  Nc=Ns=8 

MFIE  METHOD 

1.0E-04 

1.0E-06 

REAL{F(r')} 

IMAG{F(r')} 

-1.24986000612199E-01 

3. 2468589 1484467E-07 

-1.24999860001488E-01 

3.24685891594534E-09 

RIEMANN  SUM 

REAL{F(r')} 

-1. 2498600061193  IE-01 

-1.24999860007953E-01 

N  =  10,000 

IMAG{F(r')} 

3. 2468589 1454803  E-07 

3.24685891583276E-09 

CUBATURE  DIM 

SD 

RE{F(r')} 

21 

-1.24986000615963E-01 

-1.24999860001526E-01 

28 

-1.24986000612553E-01 

-1.24999860001492E-01 

36 

-1.24986000612965E-01 

-1.24999860001496E-01 

45 

-1.24986000612010E-01 

-1.24999859883931E-01 

46 

-1.24986000612687E-01 

-1.24999860001493E-01 

54 

-1.24986000612612E-01 

-1.24999860001493E-01 

55 

-1. 24986000612 120E-01 

-1.24999859923020E-01 

66 

12 

-1.24986000612220E-01 

-1.24999860001489E-01 

78 

-1.24986000612306E-01 

-1.24999860001490E-01 

85 

-1.24986000612301E-01 

-1.24999860001489E-01 

91 

15 

-1.24986000612160E-01 

-1.24999860001488E-01 

105 

-1.24986000612241E-01 

-1.24999860001489E-01 

120 

-1.249860006122 12E-01 

-1.24999860001489E-01 

126 

-1.249860006122 10E-01 

-1.24999860001489E-01 

CUBATURE  DIM 

SD 

IM{F(r')} 

21 

3. 2468589 1476039  E-07 

3. 2468589 1604692E-09 

28 

3. 2468589 1476039  E-07 

3. 2468589 160469  IE-09 

36 

3. 2468589 1476038E-07 

3. 2468589 1604691E-09 

45 

3. 2468589 1486629  E-07 

3. 24686002 148726E-09 

46 

3. 2468589 1476039  E-07 

3. 2468589 1604692E-09 

54 

10 

3. 2468589 1476040E-07 

3.24685891604690E-09 

55 

11 

3. 24685891483 108E-07 

3.24685965392274E-09 

66 

3. 2468589 1476041  E-07 

3.24685891604691E-09 

78 

3. 2468589 1476041  E-07 

3. 2468589 160469  IE-09 

85 

3. 2468589 1476038E-07 

3. 2468589 160469  IE-09 

91 

3. 2468589 1476038E-07 

3.24685891604691E-09 

105 

3. 2468589 1476040E-07 

3. 2468589 1604691E-09 

120 

3. 2468589 1476040E-07 

3. 2468589 1604690E-09 

126 

3. 2468589 1476040E-07 

3. 2468589 1604690E-09 
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Table  G.2.a:  Comparison  of  our  results  for  the  surface  integral  of  Case  38  with  those 
obtained  through  a  Riemann  sum  and  cubatures. 


z'  = 

15  SD,  Nc=Ns=8 

MFIE  METHOD 

5.0 

4.5 

REAL{F(r')} 

IMAG{F(r')} 

1.72890617968236E-03 

4. 80316359879895  E-03 

3.77040226832734E-04 

5.68301993259889E-03 

RIEMANN  SUM 

REAL{F(r')} 

1. 7289061795 1762  E-03 

3.77040226535453E-04 

N  =  10,000 

IMAG{F(r')} 

4.80316359927620E-03 

5. 683019933 10194E-03 

CUBATURE  DIM 

SD 

RE{F(r')} 

21 

1. 728906179682 19  E-03 

3.77040226832552E-04 

28 

15 

1.72890617968236E-03 

3.77040226832749E-04 

36 

1.72890617968236E-03 

3.77040226832744E-04 

45 

1.72890617968236E-03 

3.77040226832749E-04 

46 

1.72890617968236E-03 

3.77040226832742E-04 

54 

1. 72890617968235  E-03 

3.77040226832744E-04 

55 

1. 72890617968235  E-03 

3.77040226832744E-04 

66 

1.72890617968236E-03 

3.77040226832742E-04 

78 

1.72890617968236E-03 

3.77040226832749E-04 

85 

1.72890617968236E-03 

3.77040226832749E-04 

91 

1.72890617968236E-03 

3.77040226832749E-04 

105 

13 

1. 72890617968235  E-03 

3.77040226832741E-04 

120 

1.72890617968236E-03 

3.77040226832744E-04 

126 

1. 72890617968235  E-03 

3.77040226832742E-04 

CUBATURE  DIM 

SD 

IM{F(r')} 

21 

15, 15 

4. 80316359879895  E-03 

5.68301993259889E-03 

28 

4.803 16359879895  E-03 

5.68301993259889E-03 

36 

4. 80316359879895  E-03 

5.68301993259889E-03 

45 

4.80316359879896E-03 

5.68301993259890E-03 

46 

4.803 16359879895  E-03 

5.68301993259889E-03 

54 

4.803 16359879895  E-03 

5.68301993259888E-03 

55 

4. 80316359879892  E-03 

5.68301993259886E-03 

66 

4. 80316359879895  E-03 

5.68301993259889E-03 

78 

4.803 16359879895  E-03 

5.68301993259889E-03 

85 

4.803 16359879895  E-03 

5.68301993259889E-03 

91 

4. 80316359879895  E-03 

5.68301993259889E-03 

105 

4. 80316359879895  E-03 

5.68301993259889E-03 

120 

4.80316359879896E-03 

5.68301993259890E-03 

126 

4.803 16359879895  E-03 

5.68301993259889E-03 
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Table  G.2.b:  Comparison  of  our  results  for  the  surface  integral  of  Case  38  with  those 
obtained  through  a  Riemann  sum  and  cubatures. 


z'  = 

15  SD,  Nc=Ns=8 

MFIE  METHOD 

1.0E+00 

1.0E-02 

REAL{F(r')} 

IMAG{F(r')} 

-2. 61083307 127456E-02 

3. 07874207413 142  E-03 

-6.20687168677050E-02 

3.20439982334867E-05 

RIEMANN  SUM 

REAL{F(r')} 

-2. 61083307 149989E-02 

-6.20687168677565E-02 

N  =  10,000 

IMAG{F(r')} 

3. 07874207433365  E-03 

3.20439982355682E-05 

CUBATURE  DIM 

SD 

RE{F(r')} 

21 

-2. 610833072 17 141E-02 

-6.20687169957358E-02 

28 

-2. 61083307 132206E-02 

-6.20687168032204E-02 

36 

-2. 61083307 127649E-02 

-6.20687168915915E-02 

45 

-2. 61083307 127448E-02 

-6.20687168431814E-02 

46 

-2. 61083307 127475E-02 

-6.20687168821266E-02 

54 

-2. 61083307 127455E-02 

-6.20687168767932E-02 

55 

-2. 61083307 127457E-02 

-6.20687168541300E-02 

66 

15 

-2. 61083307 127456E-02 

-6.20687168688377E-02 

78 

-2. 61083307 127456E-02 

-6.20687168704573E-02 

85 

-2. 61083307 127456E-02 

-6.20687168707639E-02 

91 

-2. 61083307 127456E-02 

-6.20687168662112E-02 

105 

10 

-2. 61083307 127456E-02 

-6.20687168684284E-02 

120 

-2. 61083307 127456E-02 

-6.20687168685661E-02 

126 

-2. 61083307 127456E-02 

-6.20687168684846E-02 

CUBATURE  DIM 

SD 

IM{F(r')} 

21 

13 

3. 07874207413 143  E-03 

3.20439982334859E-05 

28 

15 

3. 07874207413 142  E-03 

3.20439982334858E-05 

36 

3. 07874207413 142  E-03 

3.20439982334859E-05 

45 

3. 07874207413 143  E-03 

3.20439982334858E-05 

46 

3. 07874207413 142  E-03 

3.20439982334859E-05 

54 

3. 07874207413 142  E-03 

3.20439982334858E-05 

55 

3. 07874207413 141  E-03 

3.20439982334857E-05 

66 

3. 07874207413 142  E-03 

3.20439982334858E-05 

78 

3. 07874207413 142  E-03 

3.20439982334858E-05 

85 

3. 07874207413 142  E-03 

3.20439982334859E-05 

91 

3. 07874207413 142  E-03 

3.20439982334859E-05 

105 

3. 07874207413 142  E-03 

3.20439982334858E-05 

120 

3. 07874207413 143  E-03 

3.20439982334859E-05 

126 

3. 07874207413 142  E-03 

3.20439982334858E-05 
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Table  G.2.c:  Comparison  of  our  results  for  the  surface  integral  of  Case  38  with  those 
obtained  through  a  Riemann  sum  and  cubatures. 


z'  = 

15  SD,  Nc=Ns=8 

MFIE  METHOD 

1.0E-04 

1.0E-06 

REAL{F(r')} 

IMAG{F(r')} 

-6.24956990989627E-02 

3. 20441256830233  E-07 

-6.24999569922107E-02 

3.20441257358174E-09 

RIEMANN  SUM 

REAL{F(r')} 

-6.24956990988243E-02 

-6.24999569919056E-02 

N  =  10,000 

IMAG{F(r')} 

3. 20441256846239  E-07 

3.20441256973727E-09 

CUBATURE  DIM 

SD 

RE{F(r')} 

21 

-6. 24956991001 147E-02 

-6.24999569922222E-02 

28 

-6.24956990982312E-02 

-6.24999569922033E-02 

36 

-6.24956990991425E-02 

-6.24999569922125E-02 

45 

-6.24956990988784E-02 

-6.24999568746547E-02 

46 

-6.24956990990698E-02 

-6.24999569922117E-02 

54 

-6.24956990990267E-02 

-6.24999569922113E-02 

55 

-6.24956990988958E-02 

-6.24999569235308E-02 

66 

-6.24956990989744E-02 

-6.24999569922108E-02 

78 

-6. 2495699098982  IE-02 

-6.24999569922109E-02 

85 

-6. 2495699098982  IE-02 

-6.24999569922109E-02 

91 

13, 15 

-6. 24956990989635  E-02 

-6.24999569922107E-02 

105 

-6.24956990989641E-02 

-6.24999569922107E-02 

120 

-6.24956990989685E-02 

-6.24999569922107E-02 

126 

-6.24956990989677E-02 

-6.24999569922107E-02 

CUBATURE  DIM 

SD 

IM{F(r')} 

21 

3. 20441256825502  E-07 

3.20441256952951E-09 

28 

3. 20441256825502  E-07 

3.20441256952952E-09 

36 

3. 20441256825502  E-07 

3.20441256952951E-09 

45 

3. 20441256836092  E-07 

3.20441367496987E-09 

46 

3. 20441256825502  E-07 

3. 2044125695295  IE-09 

54 

10 

3. 20441256825503  E-07 

3.20441256952952E-09 

55 

3. 2044125683 1688E-07 

3. 2044132 1536074E-09 

66 

3. 20441256825502  E-07 

3.20441256952952E-09 

78 

3. 20441256825503  E-07 

3.20441256952952E-09 

85 

3. 20441256825503  E-07 

3.20441256952951E-09 

91 

8 

3. 20441256825502  E-07 

3.20441256952953E-09 

105 

3. 20441256825502  E-07 

3.20441256952952E-09 

120 

3. 20441256825503  E-07 

3.20441256952952E-09 

126 

3. 20441256825503  E-07 

3.20441256952951E-09 
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Table  G.3.a:  Comparison  of  our  results  for  the  surface  integral  of  Case  39  with  those 
obtained  through  a  Riemann  sum  and  cubatures. 


z'  = 

15  SD,  Nc=Ns=8 

MFIE  METHOD 

5.0 

4.5 

REAL{F(r')} 

IMAG{F(r')} 

1.62745 189256202  E-03 

4. 95019906577089  E-03 

2. 106867572 12 199E-04 

5.83827756766797E-03 

RIEMANN  SUM 

REAL{F(r')} 

1.62745 189236729  E-03 

2.10686756873228E-04 

N  =  10,000 

IMAG{F(r')} 

4. 95019906626393  E-03 

5.83827756818758E-03 

CUBATURE  DIM 

SD 

RE{F(r')} 

21 

1.62745189256148E-03 

2.10686757211024E-04 

28 

15 

1.62745189256202  E-03 

2. 106867572 12207E-04 

36 

1.62745189256202  E-03 

2. 106867572 12207E-04 

45 

1.62745189256203  E-03 

2. 106867572 122 16E-04 

46 

1.62745189256202  E-03 

2. 106867572 12206E-04 

54 

12 

1.62745189256202  E-03 

2.10686757212201E-04 

55 

1.62745189256202  E-03 

2. 106867572 12206E-04 

66 

1.62745 189256202  E-03 

2. 106867572 12209E-04 

78 

1.62745189256202  E-03 

2. 106867572 122 13E-04 

85 

1.62745189256203  E-03 

2. 106867572 12207E-04 

91 

1.62745189256203  E-03 

2. 106867572 122 13E-04 

105 

1.62745189256202  E-03 

2. 106867572 12211E-04 

120 

1.62745 189256202  E-03 

2. 106867572 12211E-04 

126 

1.62745189256202  E-03 

2. 106867572 122 13E-04 

CUBATURE  DIM 

SD 

IM{F(r')} 

21 

15 

4. 95019906577089  E-03 

5.83827756766798E-03 

28 

15 

4. 95019906577089  E-03 

5.83827756766797E-03 

36 

4. 95019906577089  E-03 

5.83827756766797E-03 

45 

4.95019906577090E-03 

5.83827756766798E-03 

46 

4. 95019906577089  E-03 

5.83827756766797E-03 

54 

4. 95019906577089  E-03 

5.83827756766797E-03 

55 

4. 95019906577089  E-03 

5.83827756766797E-03 

66 

4. 95019906577089  E-03 

5.83827756766797E-03 

78 

4. 95019906577089  E-03 

5.83827756766797E-03 

85 

4. 95019906577089  E-03 

5.83827756766797E-03 

91 

4. 95019906577089  E-03 

5.83827756766798E-03 

105 

4. 95019906577089  E-03 

5.83827756766797E-03 

120 

4.95019906577088E-03 

5.83827756766796E-03 

126 

4. 95019906577089  E-03 

5.83827756766797E-03 
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Table  G.3.b:  Comparison  of  our  results  for  the  surface  integral  of  Case  39  with  those 


obtained  through  a  Riemann  sum  and  cubatures. 

z'  = 

1.0E+00 

1.0E-02 

15  SD,  Nc=Ns=8 

MFIE  METHOD 

REAL{F(r')} 

IMAG{F(r')} 

-3.85251154819883E-02 

3. 1408340479 1340E-03 

-2.45768216011319E-01 

3.26827981281800E-05 

RIEMANN  SUM 

REAL{F(r')} 

-3.85251154867750E-02 

-2. 457682 16011701E-01 

N  =  10,000 

IMAG{F(r')} 

3. 14083404812 180E-03 

3.26827981303243E-05 

CUBATURE  DIM 

SD 

RE{F(r')} 

21 

-3. 8525115557 1443E-02 

-2.45768238119501E-01 

28 

-3.85251154844081E-02 

-2.45768223087241E-01 

36 

-3.85251154822801E-02 

-2.45768222768720E-01 

45 

-3.85251154819970E-02 

-2. 457682 18058558E-01 

46 

-3.85251154819630E-02 

-2. 457682 11272347E-01 

54 

-3.85251154819900E-02 

-2. 457682 16765815E-01 

55 

10 

-3.85251154819879E-02 

-2. 457682 16054808E-01 

66 

15 

-3.85251154819883E-02 

-2.45768214392943E-01 

78 

-3.85251154819883E-02 

-2. 457682 17743930E-01 

85 

-3.85251154819883E-02 

-2. 457682 16711269E-01 

91 

-3.85251154819883E-02 

-2. 457682 14753590E-01 

105 

-3.85251154819883E-02 

-2. 457682 15655422E-01 

120 

-3.85251154819883E-02 

-2. 457682 17248556E-01 

126 

-3.85251154819883E-02 

-2. 457682 16240576E-01 

CUBATURE  DIM 

SD 

IM{F(r')} 

21 

15, 13 

3.14083404791340E-03 

3.26827981281858E-05 

28 

3. 1408340479 1340  E-03 

3.26827981281880E-05 

36 

3.14083404791340E-03 

3.26827981281867E-05 

45 

3. 1408340479 1341  E-03 

3.26827981281870E-05 

46 

3. 1408340479 1340  E-03 

3.26827981281865E-05 

54 

3. 1408340479 1340  E-03 

3.26827981281866E-05 

55 

3.14083404791340E-03 

3.26827981281868E-05 

66 

3.14083404791340E-03 

3.26827981281865E-05 

78 

3. 1408340479 1340  E-03 

3.26827981281874E-05 

85 

3. 1408340479 1340  E-03 

3.26827981281861E-05 

91 

3.14083404791340E-03 

3.26827981281866E-05 

105 

3. 1408340479 1340E-03 

3.26827981281858E-05 

120 

3. 1408340479 1340E-03 

3.26827981281866E-05 

126 

3. 1408340479 1340E-03 

3.26827981281862E-05 
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Table  G.3.c:  Comparison  of  our  results  for  the  surface  integral  of  Case  39  with  those 
obtained  through  a  Riemann  sum  and  cubatures. 


z'  = 

15  SD,  Nc=Ns=8 

MFIE  METHOD 

1.0E-04 

1.0E-06 

REAL{F(r')} 

IMAG{F(r')} 

-2.49957723083909E-01 

3.26829273854127E-07 

-2.49999577237499E-01 

3.26829273878434E-09 

RIEMANN  SUM 

REAL{F(r')} 

-2.49957723083307E-01 

-2.49999577256479E-01 

N  =  10,000 

IMAG{F(r')} 

3.26829273871477E-07 

3.26829274000744E-09 

CUBATURE  DIM 

SD 

RE{F(r')} 

21 

-2.49957723309702E-01 

-2.49999577239756E-01 

28 

-2.49957723158600E-01 

-2.49999577238245E-01 

36 

-2.49957723152704E-01 

-2.49999577238186E-01 

45 

-2.49957723106961E-01 

-2.49999577237729E-01 

46 

-2.49957723035387E-01 

-2.49999577237013E-01 

54 

-2.49957723091232E-01 

-2.49999577237572E-01 

55 

11, 14 

-2. 49957723084043  E-01 

-2.49999577237500E-01 

66 

-2.49957723066964E-01 

-2.49999577237329E-01 

78 

-2.49957723102660E-01 

-2.49999577237686E-01 

85 

-2.49957723091748E-01 

-2.49999577237577E-01 

91 

-2.49957723070645E-01 

-2.49999577237366E-01 

105 

-2.49957723081466E-01 

-2.49999577237474E-01 

120 

-2.49957723097360E-01 

-2.49999577237633E-01 

126 

-2.49957723086500E-01 

-2.49999577237524E-01 

CUBATURE  DIM 

SD 

IM{F(r')} 

21 

3.26829273850200E-07 

3.26829273979439E-09 

28 

11 

3.26829273850206E-07 

3.26829273979443E-09 

36 

3.26829273850184E-07 

3.26829273979447E-09 

45 

3.26829273850196E-07 

3.26829273979437E-09 

46 

3.26829273850186E-07 

3.26829273979442E-09 

54 

3.26829273850187E-07 

3.26829273979442E-09 

55 

3. 26829273850185  E-07 

3.26829273979444E-09 

66 

3.26829273850186E-07 

3.26829273979443E-09 

78 

3. 26829273850185  E-07 

3.26829273979442E-09 

85 

3. 26829273850265  E-07 

3. 268292739797 17E-09 

91 

3.26829273850184E-07 

3.26829273979442E-09 

105 

3.26829273850156E-07 

3.26829273979450E-09 

120 

9 

3. 26829273850185  E-07 

3. 26829273979435  E-09 

126 

3.26829273850184E-07 

3.26829273979440E-09 
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Table  G.4.a:  Comparison  of  our  results  for  the  surface  integral  of  Case  40  with  those 
obtained  through  a  Riemann  sum  and  cubatures. 


z'  = 

15  SD,  Nc=Ns=8 

MFIE  METHOD 

5.0 

4.5 

REAL{F(r')} 

IMAG{F(r')} 

1.63265 177430638E-03 

4. 94291929209383  E-03 

2.19197452706583E-04 

5.83059487773856E-03 

RIEMANN  SUM 

REAL{F(r')} 

1.63265 177408649  E-03 

2.19197452322990E-04 

N  =  10,000 

IMAG{F(r')} 

4. 94291929265363  E-03 

5.83059487832946E-03 

CUBATURE  DIM 

SD 

RE{F(r')} 

21 

1. 63265 177430584E-03 

2.19197452705390E-04 

28 

15 

1.63265177430639  E-03 

2.19197452706573E-04 

36 

1.63265 177430639  E-03 

2.19197452706573E-04 

45 

1. 63265 177430640E-03 

2.19197452706578E-04 

46 

1.63265 177430639  E-03 

2. 1919745270657  IE-04 

54 

1.63265 177430639  E-03 

2.19197452706573E-04 

55 

1.63265 177430639  E-03 

2.19197452706573E-04 

66 

1.63265 177430639  E-03 

2.19197452706569E-04 

78 

1.63265 177430639  E-03 

2.19197452706575E-04 

85 

1. 63265 177430640E-03 

2.19197452706578E-04 

91 

14 

1. 63265 177430640E-03 

2.19197452706580E-04 

105 

1.63265 177430639  E-03 

2.19197452706569E-04 

120 

1.63265 177430639  E-03 

2.19197452706573E-04 

126 

1.63265 177430639  E-03 

2.19197452706561E-04 

CUBATURE  DIM 

SD 

IM{F(r')} 

21 

15 

4.94291929209384E-03 

5.83059487773856E-03 

28 

15 

4. 9429 1929209383  E-03 

5.83059487773856E-03 

36 

4. 94291929209383  E-03 

5.83059487773856E-03 

45 

4.94291929209384E-03 

5.83059487773857E-03 

46 

4. 9429 1929209383  E-03 

5.83059487773856E-03 

54 

4. 9429 1929209383  E-03 

5.83059487773856E-03 

55 

4. 94291929209383  E-03 

5.83059487773856E-03 

66 

4. 94291929209383  E-03 

5.83059487773856E-03 

78 

4. 9429 1929209383  E-03 

5.83059487773856E-03 

85 

4. 9429 1929209383  E-03 

5.83059487773856E-03 

91 

4. 94291929209383  E-03 

5.83059487773856E-03 

105 

4. 94291929209383  E-03 

5.83059487773856E-03 

120 

4. 9429 1929209382  E-03 

5.83059487773854E-03 

126 

4. 9429 1929209383  E-03 

5.83059487773856E-03 
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Table  G.4.b:  Comparison  of  our  results  for  the  surface  integral  of  Case  40  with  those 
obtained  through  a  Riemann  sum  and  cubatures. 


z'  = 

15  SD,  Nc=Ns=8 

MFIE  METHOD 

1.0E+00 

1.0E-02 

REAL{F(r')} 

IMAG{F(r')} 

-3.75158019087259E-02 

3. 13776649 153445  E-03 

-2.93455187090841E-02 

3.26512413098946E-05 

RIEMANN  SUM 

REAL{F(r')} 

-3. 75158019 140408E-02 

-2.93455187120835E-02 

N  =  10,000 

IMAG{F(r')} 

3.13776649177036E-03 

3.26512413123283E-05 

CUBATURE  DIM 

SD 

RE{F(r')} 

21 

-3.75158019812679E-02 

-2.93455303684519E-02 

28 

-3. 7515801911065  IE-02 

-2.93455223333366E-02 

36 

-3.75158019090099E-02 

-2. 934552 18856173E-02 

45 

-3.75158019087327E-02 

-2.93455194872424E-02 

46 

-3.75158019087016E-02 

-2.93455165345567E-02 

54 

-3.75158019087273E-02 

-2.93455186197702E-02 

55 

-3.75158019087256E-02 

-2.93455184475968E-02 

66 

15 

-3.75158019087259E-02 

-2.93455180058965E-02 

78 

-3.75158019087259E-02 

-2.93455 191161826E-02 

85 

-3.75158019087259E-02 

-2.93455188876310E-02 

91 

-3.75158019087259E-02 

-2.93455182707200E-02 

105 

-3.75158019087259E-02 

-2.93455186095185E-02 

120 

-3.75158019087259E-02 

-2.93455189809726E-02 

126 

12 

-3.75158019087259E-02 

-2.93455187119994E-02 

CUBATURE  DIM 

SD 

IM{F(r')} 

21 

3. 13776649 153446E-03 

3.26512413098932E-05 

28 

3. 13776649 153446E-03 

3.26512413098933E-05 

36 

15, 13 

3. 13776649153445  E-03 

3.26512413098934E-05 

45 

3.13776649153446E-03 

3.26512413098934E-05 

46 

3. 13776649 153445  E-03 

3.26512413098932E-05 

54 

3. 13776649 153446E-03 

3.26512413098933E-05 

55 

3.13776649153446E-03 

3.26512413098933E-05 

66 

3.13776649153446E-03 

3.26512413098934E-05 

78 

3. 13776649 153446E-03 

3.26512413098934E-05 

85 

3. 13776649 153446E-03 

3.26512413098933E-05 

91 

3. 13776649 153446E-03 

3.26512413098933E-05 

105 

3. 13776649153445  E-03 

3.26512413098933E-05 

120 

3. 13776649 153445  E-03 

3.26512413098932E-05 

126 

3. 13776649 153445  E-03 

3.26512413098933E-05 
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Table  G.4.c:  Comparison  of  our  results  for  the  surface  integral  of  Case  40  with  those 


obtained  through  a  Riemann  sum  and  cubatures. 

z'  = 

1.0E-04 

1.0E-06 

15  SD,  Nc=Ns=8 

MFIE  METHOD 

REAL{F(r')} 

IMAG{F(r')} 

-2.98383298895019E-04 

3.26513704776294E-07 

-2.98383804956031E-06 

3.26513704905455E-09 

RIEMANN  SUM 

REAL{F(r')} 

-2. 983832989262 11E-04 

-2.98383804987222E-06 

N  =  10,000 

IMAG{F(r')} 

3.26513704800458E-07 

3.26513704929694E-09 

CUBATURE  DIM 

SD 

RE{F(r')} 

21 

-2.98383416437008E-04 

-2.98383922498115E-06 

28 

-2.98383335647856E-04 

-2.98383841708918E-06 

36 

-2.98383331021954E-04 

-2.98383837083002E-06 

45 

-2.98383306868899E-04 

-2.98383812929930E-06 

46 

-2.98383276881133E-04 

-2.98383782942117E-06 

54 

-2.98383297940533E-04 

-2.98383804001537E-06 

55 

-2.98383296208764E-04 

-2.98383802269768E-06 

66 

-2. 9838329 1728056E-04 

-2.98383797789053E-06 

78 

-2.98383303054725E-04 

-2.98383809115745E-06 

85 

-2.98383300707697E-04 

-2.98383806768710E-06 

91 

-2.98383294414452E-04 

-2.98383800475453E-06 

105 

-2.98383297902541E-04 

-2.98383803963552E-06 

120 

-2.98383301680386E-04 

-2.98383807741404E-06 

126 

9, 10 

-2.98383298924706E-04 

-2.98383804985717E-06 

CUBATURE  DIM 

SD 

IM{F(r')} 

21 

13 

3.265 13704776282  E-07 

3.26513704905445E-09 

28 

13 

3.26513704776278E-07 

3.26513704905448E-09 

36 

3.26513704776280E-07 

3.26513704905446E-09 

45 

3.26513704776280E-07 

3.26513704905448E-09 

46 

3. 265 13704776279  E-07 

3.26513704905447E-09 

54 

3.26513704776278E-07 

3.26513704905447E-09 

55 

3. 265 13704776279  E-07 

3.26513704905446E-09 

66 

3. 265 13704776279  E-07 

3.26513704905446E-09 

78 

3. 265 13704776279  E-07 

3.26513704905447E-09 

85 

3. 265 13704776279  E-07 

3.26513704905446E-09 

91 

3. 265 13704776279  E-07 

3.26513704905447E-09 

105 

3.26513704776278E-07 

3.26513704905447E-09 

120 

3.26513704776278E-07 

3.26513704905446E-09 

126 

3.265 13704776279  E-07 

3.26513704905446E-09 
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APPENDIX  H 

LINE  INTEGRAL  RESULTS  FOR  REST  OF  CASES  OF  TABLE  6.14 

In  this  Appendix,  we  present  the  rest  of  the  line  integral  results  for  the  points  of  Table  6.14.  The 
first  result,  Case  38,  was  presented  in  Table  6.19. 

Table  H.l :  The  line  integral  of  (6.3)  for  an  OP  above  the  interior  point  (-7/30,  -7/30) 
of  the  triangle.  The  minimum  distance  of  this  point  to  a  side  is  10  mA.  Despite  the  small 
distance,  the  G-K  quadrature  remains  stable  in  the  approach  to  this  point  and  agrees 

with  our  method  to  at  least  14  SD. 


Case 

36 

PRESENT  METHOD 

MATLAB  G-K  QUADRATURE 

z' 

Side 

RE{INTG 1> 

IM{INTG„1> 

RE{INTG 1> 

IM{INTG 1} 

5.0 

S  1 

-2.79991492725852E-01 

8.51159900317967E-03 

-2.79991492725852E-01 

8. 51 1599003 17923E-03 

S  2 

-1.98953683726784E-01 

3. 1 37393530 10368E-03 

-1.98953683726784E-01 

3. 137393530 1 0371 E-03 

S  3 

-1.98953683726784E-01 

3. 1 37393530 10368E-03 

-1.98953683726784E-01 

3. 137393530 10371 E-03 

1.0 

S  1 

8. 404962 177542 53 E-01 

-8.04197056952424E-01 

8. 4049621 77542 54E-01 

-8.04197056952425E-01 

S  2 

6.93958697338809E-01 

-5. 778 167902 12706E-01 

6.93958697338809E-01 

-5. 778167902 12707E-01 

S  3 

6.93958697338809E-01 

-5. 778 167902 12706E-01 

6.93958697338809E-01 

-5. 778167902 12707E-01 

0.01 

S  1 

1.90555978762984E+00 

-8. 604 14763093 154E-01 

1.90555978762985E+00 

-8. 60414763093 154E-01 

S  2 

3.68186409929430E+00 

-6.17941085308805E-01 

3. 68 18640992943 1E+00 

-6. 17941085308805E-01 

S  3 

3.68186409929430E+00 

-6.17941085308805E-01 

3. 68 18640992943 1E+00 

-6. 17941085308805E-01 

l.E-04 

S  1 

1.90582347488829E+00 

-8. 604204976 17676E-01 

1.90582347488830E+00 

-8.60420497617676E-01 

S  2 

3.6903582200468 1E+00 

-6.17945177996947E-01 

3.69035822004681E+00 

-6. 17945177996947E-01 

S  3 

3.6903582200468 1E+00 

-6.17945177996947E-01 

3. 69035822004681 E+00 

-6. 17945177996947E-01 

l.E-06 

S  1 

1. 90582350 126170E+00 

-8.60420498191 129E-01 

1. 90582350 126 170E+00 

-8.60420498191 130E-01 

S  2 

3. 6903590741 3992E+00 

-6. 17945 1784062 17E-01 

3. 690359074 13992E+00 

-6. 17945 1784062 17E-01 

S  3 

3. 6903590741 3992E+00 

-6. 17945 1784062 17E-01 

3. 690359074 13992E+00 

-6. 17945 1784062 17E-01 
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Table  H.2:  The  line  integral  of  (6.3)  for  an  OP  above  the  90-deg  corner  of  the  triangle. 
As  the  OP  approaches  this  point,  the  G-K  quadrature  begins  to  lose  SD  over  Sides  2  and  3, 
the  two  sides  that  have  this  comer  in  common.  As  expected,  the  loss  occurs  in  the  real  part 

of  the  integral  and  its  maximum  is  3  SD. 


Case 

37 

PRESENT  METHOD 

MATLAB  G-K  QUADRATURE 

z' 

Side 

RE{II\ITG 1} 

IM{INTG 1} 

RE{INTG 1> 

IM{INTGJ} 

5.0 

S  1 

-2.78895727915674E-01 

1.16009605809 162E-02 

-2. 7889572791 5675E-01 

1.16009605809165E-02 

S  2 

-1.98613359322592E-01 

4. 11400306068684E-03 

-1.98613359322592E-01 

4. 11400306068682E-03 

S  3 

-1.98613359322592E-01 

4. 11400306068684E-03 

-1.98613359322592E-01 

4. 11400306068682E-03 

1.0 

S  1 

7.57961656479037E-01 

-7.943 13693 10 1739E-01 

7.57961656479037E-01 

-7.9431 3693 101740E-01 

S  2 

6.64800195907522E-01 

-5.74698696029768E-01 

6.64800195907522E-01 

-5.74698696029768E-01 

S  3 

6.64800195907522E-01 

-5.74698696029768E-01 

6.64800195907522E-01 

-5.74698696029768E-01 

0.01 

S  1 

1.54106739353266E+00 

-8. 50 1284262764 17E-01 

1.54106739353267E+00 

-8. 50 1284262764 17E-01 

S  2 

5. 201 1988933535 1E+00 

-6.14696002448062E-01 

5.20119889335352E+00 

-6.14696002448062E-01 

S  3 

5. 201 1988933535 1E+00 

-6. 14696002448062E-01 

5.20119889335352E+00 

-6. 14696002448062E-01 

l.E-04 

S  1 

1.54122507689600E+00 

-8.50134119929503E-01 

1.54122507689600E+00 

-8.50134119929504E-01 

S  2 

9.80640081789547E+00 

-6.14700082256520E-01 

9. 8064008 1789577E+00 

-6. 1470008225652  IE-01 

S  3 

9.80640081789547E+00 

-6.14700082256520E-01 

9. 8064008 1789577E+00 

-6. 1470008225652 IE-01 

l.E-06 

S  1 

1. 54 1225092666 18E+00 

-8.50134120498870E-01 

1. 541225092666 18E+00 

-8.50134120498870E-01 

S  2 

1. 44 1157101160 19E+01 

-6. 1470008266450  IE-01 

1.44115710116314E  +  01 

-6.14700082664502E-01 

S  3 

1.44 1157 101160 19E+01 

-6. 1470008266450  IE-01 

1.44115710116314E  +  01 

-6.14700082664502E-01 
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Table  H.3:  The  line  integral  of  (6.3)  for  an  OP  above  the  midpoint  of  the  triangle's  hypotenuse. 
As  the  OP  approaches  this  point,  the  G-K  quadrature  begins  to  lose  SD  over  Side  1,  the  side 
that  contains  this  point.  As  expected,  the  loss  occurs  in  the  real  part  of  the  integral  and  its 

maximum  is  3  SD. 


Case 

39 

PRESENT  METHOD 

MATLAB  G-K  QUADRATURE 

z' 

Side 

RE{INTG„1> 

IM{INTG 1> 

RE{INTG„1> 

IM{INTG 1> 

5.0 

S  1 

-2. 8 1880767906 103E-01 

2. 935372943875 15E-03 

-2.81880767906104E-01 

2.93537294387527E-03 

S  2 

- 1.98635 108609686E-01 

4. 144991 67474095 E-0 3 

- 1. 98635 108609687E-01 

4. 14499 1674741 29  E-03 

S  3 

- 1.98635 108609686E-01 

4. 144991 67474095 E-0  3 

-1. 98635 108609687E-01 

4. 14499 1674741 29  E-03 

1.0 

S  1 

1. 027494 14262907E+00 

-8.2 1944754365 135E-01 

1.027494 14262907E+00 

-8. 2 1944754365 135E-01 

S  2 

6. 49008525848461 E-01 

-5.74633334885979E-01 

6.49008525848462E-01 

-5.74633334885979E-01 

S  3 

6. 49008525848461 E-01 

-5.74633334885979E-01 

6.49008525848462E-01 

-5.74633334885979E-01 

0.01 

S  1 

9.80559256320577E+00 

-8.78883957009733E-01 

9.80559256320577E+00 

-8.78883957009734E-01 

S  2 

1. 65043 165737456E+00 

-6. 1462876 1354497E-01 

1. 65043 165737457E+00 

-6. 1462876 1354498E-01 

S  3 

1. 65043 165737456E+00 

-6. 1462876 1354497E-01 

1. 65043 165737457E+00 

-6. 1462876 1354498E-01 

l.E-04 

S  1 

1.901 59400618 196E+01 

-8. 788897647 10227E-01 

1.90 1594006 18 196E+01 

-8. 788897647 10228E-01 

S  2 

1.65073 12407526 1E+00 

-6. 14632840972885E-01 

1. 650731 24075262E+00 

-6. 14632840972886E-01 

S  3 

1.65073 12407526 1E+00 

-6.14632840972885E-01 

1. 650731 24075262E+00 

-6. 14632840972886E-01 

l.E-06 

S  1 

2.82262804435970E+01 

-8.78889765290998E-01 

2. 8226280573558 1E+01 

-8.78889765290999E-01 

S  2 

1.65073 12707 18 16E+00 

-6. 1463284 1380848E-01 

1.650731 270718 17E+00 

-6. 14632841380848E-01 

S  3 

1.65073 12707 18 16E+00 

-6. 1463284 1380848E-01 

1.650731 270718 17E+00 

-6. 14632841380848E-01 
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Table  H.4:  The  line  integral  of  (6.3)  for  an  OP  above  the  exterior  point  1/5,  1/5)  of  the  triangle. 
The  minimum  distance  of  this  point  to  the  hypotenuse  is  4.71  mA.  Despite  the  small  distance, 
the  G-K  quadrature  remains  stable  in  the  approach  to  this  point  and  agrees  with  our  method 

to  15  SD. 


Case 

40 

PRESENT  METHOD 

MATLAB  G-K  QUADRATURE 

z' 

Side 

RE{TNTG 1} 

IM{INTG 1} 

RE{INTG 1> 

IM{INTG 1> 

5.0 

S  1 

-2.81867912873243E-01 

2.9744697 1612503E-03 

-2.81867912873244E-01 

2.9744697 16 12546E-03 

S  2 

-1.98486089338058E-01 

4.58230465634543E-03 

-1.98486089338059E-01 

4.58230465634484E-03 

S  3 

-1.98486089338058E-01 

4.58230465634543E-03 

-1.98486089338059E-01 

4.58230465634484E-03 

1.0 

S  1 

1. 0259565 1253827E+00 

-8.21820720358887E-01 

1. 0259565 1253827E+00 

-8.21820720358887E-01 

S  2 

6. 35056002 113453E-01 

-5.73240619963098E-01 

6.350560021 13454E-01 

-5.73240619963099E-01 

S  3 

6. 35056002 113453E-01 

-5.73240619963098E-01 

6.350560021 13454E-01 

-5.73240619963099E-01 

0.01 

S  1 

6.66104112899468E+00 

-8. 787548900 18984E-01 

6.661041 12899468E+00 

-8. 787548900 18985E-01 

S  2 

1.55348859766489E+00 

-6.13179408371834E-01 

1.55348859766489E+00 

-6.13179408371834E-01 

S  3 

1.55348859766489E+00 

-6. 13179408371834E-01 

1.55348859766489E+00 

-6. 13179408371834E-01 

l.E-04 

S  1 

6. 7050203048359 1E+00 

-8.78760697209028E-01 

6. 70502030483591 E+00 

-8.78760697209028E-01 

S  2 

1.55374448005436E+00 

-6. 13183482245833E-01 

1.55374448005436E+00 

-6.13183482245833E-01 

S  3 

1.55374448005436E+00 

-6. 13183482245833E-01 

1.55374448005436E+00 

-6. 13183482245833E-01 

l.E-06 

S  1 

6. 70502480 107593E+00 

-8.78760697789748E-01 

6. 70502480 107593E+00 

-8.78760697789748E-01 

S  2 

1.55374450564806E+00 

-6. 13 183482653221 E-01 

1.55374450564806E+00 

-6. 13 183482653221  E-01 

S  3 

1.55374450564806E+00 

-6. 13 183482653221 E-01 

1.55374450564806E+00 

-6. 13 183482653221 E-01 
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